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FOREWORD 


The  extension  and  application  of  the  ZEUS  code  to  high,  hypersonic  Mach  numbers  by  incorporat¬ 
ing  high-temperature  equilibrium  air  effects  is  described.  The  ZEUS  code  is  a  space-marching  Euler 
solver.  It  uses  an  automatic  multiple  zone  grid  generation  technique  and  a  finite  volume,  second-order 
Godunov  scheme  for  steady  supersonic  flow.  High-temperature  effects  are  achieved  by  using  curve 
fits  of  the  thermodynamic  properties  of  equilibrium  air.  Sample  calculations  for  a  sharp  cone  and  a 
blunt  nosed  cylinder  are  provided. 


This  work  was  supported  by  the  Naval  Weapons  Center  at  China  Lake. 


CARL  W.  LARSON,  Head 
Physics  &  Technology  Division 


NAVSWC  TR  91-86 


CONTENTS 

Chapter  Page 

1  INTRODUCTION  . 1 

2  HIGH-TEMPERATURE  EQUILIBRIUM  AIR . 3 

2.1  COMPUTED  THERMODYNAMIC  PROPERTIES  . 3 

2.2  EFFECTIVE  GAMMA . 4 

3  NUMERICAL  PROCEDURE  . 5 

3.1  GOVERNING  EQUATIONS . 5 

3.2  DECODING . 6 

3.3  NEW  DENSITY  . 7 

3.4  OBLIQUE  SHOCKS . 9 

3.5  EXPANSION  WAVES . 10 

3.6  SHOCK  ANGLE .  10 

4  ZEUS  INPUT  FILE  (TAPE  5) .  13 

4.1  INTEGRATION  CONTROL .  13 

4.2  BOUNDARY  CONDITIONS .  14 

4.3  OUTPUT  CONTROL .  15 

4.4  AERODYNAMIC  DATA .  16 

4.5  SEPARATION  MODELING  .  16 

4.6  GAS  TYPE .  17 

4.7  UNITS . 17 

4.8  INITIAL  PLANE  DATA  .  17 

4.9  MESH  CLUSTERING .  18 

5  PRE-PROCESSING  PROGRAMS  . .  19 

5.1  BLUNT2 .  19 

iii 


NAVSWC  TR  91-86 


CONTENTS  (Cont.) 

Chapter  Page 


6  SAMPLE  CASES . 21 

6.1  HYPERSONIC  SHARP  CONE . 21 

6.1.1  Geometry . 21 

6.1.2  ZEUS  Input  File  —  I  . 23 

6.1.3  ZEUS  Input  File  —  II . 27 

6.1.4  Results . 31 

6.2  HYPERSONIC  BLUNT  NOSE  CYLINDER . 31 

6.2.1  Geometry . 31 

6.2.2  START  File . 33 

6.2.3  ZEUS  Input  File . 34 

6.2.4  Results . 38 

6.3  HYPERSONIC  AIRCRAFT  . 39 

6.3.1  Geometry  —  Nose  Section . 39 

6.3.2  Zeus  Input  File  —  Nose  Section . 41 

6.3.3  Geometry  —  Wing  Section  . 45 

6.3.4  CONVERT  Input  File . 49 

6.3.5  Zeus  Input  File  —  Wing  Section . 50 

6.3.6  Results . 54 

6.4  CONCLUDING  REMARKS . 54 

REFERENCES . 65 

APPENDIX 

A  ZEUS  INPUT  FILE  . A-1 

B  SAMPLE  CASE  OUTPUT . B-1 

B.1  HYPERSONIC  SHARP  CONE . B-1 

B.1.1  ZEUS  Output  —  II . B-1 

B.2  HYPERSONIC  BLUNT  NOSE  CYLINDER . B-6 

B.2.1  BLUNT2  Output . B-6 

B.2.2  ZEUS  Output . B-11 


IV 


NAVSWC  TR  91-86 


ILLUSTRATIONS 

Page 


CARTESIAN  AND  CYLINDRICAL  COORDINATE  SYSTEM  FOR  ZEUS,  A 
SPACE-MARCHING  EULER  SOLVER . 55 

CONTROL  VOLUME  NOMENCLATURE . 56 

FLOW  DISCONTINUITIES  IN  THE  CASE  ^  >  1  (A.)  OBLIQUE  SHOCK  WAVE 
(B.)  EXPANSION  FAN . 57 

SURFACE  (A.)  PRESSURES  (B.)  DENSITIES  (C.)  TEMPERATURES  ON 
CONES  AT  AN  ALTITUDE  OF  100,000  FT.  A  COMPARISON  BETWEEN 
CALORICALLY  PERFECT  AND  EQUILIBRIUM  AIR  RESULTS  FOR  ZEUS  AND 
HUDGINS’®  COMPUTATIONS . 58 

COMPUTED  PRESSURE  CONTOURS  FOR  A  SHARP  CONE  WITH  A  25® 
HALF-ANGLE  AT  MACH  20,  AN  INCIDENCE  OF  20®,  AND  AN  ALTITUDE  OF 
100,000  FT.  RESULTS  ARE  COMPARED  BETWEEN  (A.)  EQUILIBRIUM  AIR 
AND  (B.)  CALORICALLY  PERFECT . 59 

DENSITY  CONTOURS  FOR  A  HYPERSONIC  SPHERICALLY  BLUNT  NOSED 
CYLINDER  AT  MACH  20  AND  AN  INCIDENCE  OF  0®  FOR  (A.)  CALORICALLY 
PERFECT  AND  (B.)  EQUILIBRIUM  CHEMICALLY  REACTING  AIR.  RESULTS 
ARE  COMPARED  WITH  COMPUTATIONS  OF  ZEUS  AND  PALMER’® . 60 


TEMPERATURE  CONTOURS  FOR  A  HYPERSONIC  SPHERICALLY  BLUNT 
NOSED  CYLINDER  AT  MACH  20  AND  AN  INCIDENCE  OF  0®  FOR  (A.) 
CALORICALLY  PERFECT  AND  (B.)  EQUILIBRIUM  CHEMICALLY  REACTING 
AIR.  RESULTS  ARE  COMPARED  WITH  COMPUTATIONS  OF  ZEUS  AND 
PALMER’®  . 61 

COMPUTED  DENSITY  CONTOURS  FOR  A  HYPERSONIC  SPHERICALLY 
BLUNT  NOSED  CYLINDER  AT  MACH  20  AND  AN  INCIDENCE  OF  20®  FOR 
EQUILIBRIUM  AIR  AT  AN  AXIAL  STATION  9.5  DIAMETERS  DOWNSTREAM  62 


NAVSWC  TR  91-86 


ILLUSTRATIONS  (CONT.) 


Figure  Page 

9  COMPARISON  OF  COMPUTED  PRESSURE  PROFILES  FOR  EQUILIBRIUM  AIR 
ON  THE  WINDWARD  SYMMETRY  PLANE  AT  MACH  24.5,  AN  INCIDENCE  OF 
1°  AND  Z  =  108.67  FT  FOR  ZEUS  AND  THE  PNS  SOLUTIONS  OF  KORTE^^  63 

10  COMPARISON  OF- COMPUTED  PRESSURE  CONTOURS  FOR  EQUILIBRIUM 

AIR  AT  MACH  24.5,  AN  INCIDENCE  OF  1°  AND  Z  =  108.67  FT  FOR  ZEUS 
AND  THE  PNS  SOLUTIONS  OF  KORTE’^ . 64 


VI 


NAVSWC  TR  91-86 


CHAPTER  1 
INTRODUCTION 


This  report  describes  the  extension  and  application  of  the  ZEUS  code  to  high,  hypersonic  Mach 
numbers  by  incorporating  high-temperature  equilibrium  air  or  “real  gas”  effects.  This  approach 
is  important  for  accurate  aerodynamic  predictions  in  support  of  atmospheric  reentry  vehicles  with 
high  maneuverability  and  future  hypersonic  weapons.  The  program  employs  an  automatic  multiple 
zone  grid  generation  technique  which  is  particularly  suited  to  describing  missile-type  geometries  in 
supersonic  flight.  The  ZEUS  code,'“^  whose  numerical  algorithm  is  well  described  in  References  1-2 
and  applications  in  References  2-9,  is  a  supersonic/hypersonic  space-marching  Euler  solver  which 
is  capable  of  computing  internal  as  well  as  external  flows.  The  present  report  provides  recent 
modifications  to  the  program  (described  in  Reference  2)  for  real  gas  effects  and  user  instructions 
for  its  application. 

The  ZEUS  code  uses  a  zonal,  finite  volume,  second  order  Godunov  scheme  to  integrate  the  Euler 
equations.  The  computation  proceeds  for  an  initial  data  plane  on  which  the  flow  field  is  specified,  as 
shown  in  Figure  1,  and  marches  in  a  streamwise  direction  until  the  end  of  the  region  of  interest  is 
encountered.  In  the  case  of  a  missile,  the  initial  data  plane  is  taken  near  the  nose  of  the  missile,  and 
the  flow  field  on  this  plane  can  be  defined  using  conical  flow,  a  blunt  body  solution,  or  uniform  flow. 
The  calculation  is  marched  down  the  axis  of  the  missile  until  the  base  of  the  missile  is  encountered. 
Throughout  the  computation,  the  flow  field  must  remain  supersonic  everywhere. 

At  high  flight  Mach  numbers  and  altitudes,  computation  of  the  flow  field  in  local  chemical  and 
thermodynamic  equilibrium  is  necessary  to  account  for  the  high-temperature  or  real  gas  effects.  These 
high  temperature  effects  are  taken  into  account  by  using  the  curve  fits  of  Scrinivansan  et  al.'®  to 
compute  the  thermodynamic  properties  of  equilibrium  air  as  well  as  modifying  parts  of  the  ZEUS 
code. 
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CHAPTER  2 

HIGH-TEMPERATURE  EQUILIBRIUM  AIR 


The  computation  of  missiles  or  atmospheric  reentry  bodies  at  high  hypersonic  Mach  numbers 
and  high  altitudes  requires  the  flow  field  to  be  calculated  in  local  thermodynamic  and  chemical 
equilibrium.  To  solve  an  equilibrium  flow  field,  two  thermodynamic  state  variables  (such  as  T 
and  p)  need  to  be  expressed  in  terms  of  two  other  state  variables  (such  as  p  and  h).''  The  high- 
temperature  thermodynamic  properties  of  equilibrium  chemically  reacting  air  are  obtained  from 
statistical  thermodynamics  as  discussed  in  Reference  11.  To  simplify  matters,  tables  of  thermodynamic 
properties  of  high-temperature  air  have  been  constructed  (see  Hilsenrath  and  Klein^^)  which  can  be  used 
in  numerical  flow  field  calculations  via  a  “table  look-up”  procedure  which  interpolates  between  discrete 
entries  from  the  tables.  Furthermore,  polynomial  curve  fits  of  the  tabulated  thermodynamic  data  for 
equilibrium  air  as  discussed  in  Reference  10  have  been  correlated  for  computational  convenience  and 
efficiency.  This  latter  method  is  the  approach  implemented  in  the  ZEUS  code. 


2.1  COMPUTED  THERMODYNAMIC  PROPERTIES 

The  tabulated  thermodynamic  properties  of  equilibrium  air  are  curve  fit  as  discussed  in  Refcnce 
10  and  are  implemeted  in  ZEUS  through  the  equilibrium  air  subroutines  used  in  NASA  Langley’s 
CFL3DE  program,'^  a  time-dependent  Navier-Stokes  program.  Each  subroutine  correlates  different 
properties.  For  instance,  subroutine; 


•  TGASl  —  computes  p,a,T  as  a  function  of  e  and  p. 

•  TGAS2  —  computes  s  as  a  function  of  e  and  p. 

•  TGAS4  —  computes  as  a  function  of  p  and  p. 

•  1GAS5  —  computes  p  as  a  function  of  p  and  s. 

•  TGAS7  —  computes  a  as  a  function  of  p  and  s. 


where,  p  =  pressure,  p  =  density,  T  =  temperature,  e  =  internal  energy,  h  =  enthalpy,  5  =  entropy, 
and  a  =  sound  speed. 
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2.2  EFFECTIVE  GAMMA 

It  is  convenient  to  define  an  “effective  gamma”  as  follows: 

This  is  used  in  decoding  the  flow  variables  as  discussed  in  Section  3.2. 
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CHAPTER  3 

NUMERICAL  PROCEDURE 


V"  =  ^  (5) 

Here,  U  is  the  flux  in  the  ^z-direction  w-hich  passes  through  the  shaded  cell  ends  while  the  F’s  are 
the  fluxes  associated  with  the  remaining  cell  edges.  The  >4’s  are  the  cell  edge  areas,  h  is  the  vector 
normal  to  the  cell,  and  u.  v.  w  are  the  cartesian  velocity  components.  Equations  (2)  -  (5)  are  closed 
using  the  constant  total  enthalpy  condition 

Ho  =  h  +  -{ir  +  v‘  +  u'‘)  (6) 

where  h=h(p,p). 
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3.2  DECODING 


After  the  corrector  step,  the  flow  variables  p,  p,  u,  v,  w  must  be  determined  or  decoded  from  the 
computed  conservative  variables  f/j  ,  U2  ,0^  ,  [/4,  as  displayed  in  Equation  (3).  This  is  done  using  an 
iterative  procedure  similar  to  that  used  by  Solomon  et.  al.*^ 

By  definition,  we  have 

U  =  Uzllh  (7) 


r  =  lUlUi 


(8) 


P  —  (t'2  —  Uiw)jA 


(9) 


p  =  Ihlvo/A 

Substitute  Equations  (7)  -  (10)  into  Equation  (6)  which  yields; 


where: 


H'  =  2Ho-~ 


I'l  +  t'l 


Substitue  Equations  (9-11)  into  Equation  (1)  to  get  a  quadratic  for  w. 
corresponding  to  iv~  >  Tp/p  is 

(2[r  + 


IV  = 


^iir  +  i] 


where; 


^)  =  (r2-i) 


(10) 

(11) 

(12) 

The  root  of  this  equation 

(13) 

(14) 


In  the  case  of  a  perfect  gas,  Y  =  -),  which  is  known,  and  Equations  (7)  -  (10),  and  (13)  provide 
the  desired  decoding  formulas.  For  a  more  general  gas  (such  as  equilibrium  air  )  where  h=h(p,p)  is 
curve  fit  to  thermodynamic  data,  F  is  not  known  a-priori,  and  the  decoding  can  not  be  done  in  closed 
form.  In  this  case,  the  decoding  is  done  by  solving  the  nonlinear  equation 

h{p.p)-^-[H'  -rv^]=0  (15) 

for  w.  In  the  code.  Equation  (15)  is  solved  iteratively  using  the  Secant  Method.  This  is  done  in  the 
following  manner; 


(a.)  Guess  a  value  of  F 1  from  p,  p,  h  computed  in  the  predictor  step  and  substitute  F 1  into  Equation 
(13)  to  get  w. 
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(b.)  Decode  values  of  u,  v  and  compute  new  values  of  p,  p  with  the  value  of  w  computed  in  Step  (a.) 

(c.)  Solve  Equation  (15)  by  using  the  equilibrium  air  subroutine,  TGAS4,  for  h(p,p).  This  gives  the 
first  error  estimate,  ffPi). 

(d.)  Use  the  values  of  p,  p,  h  from  Step  (b.)  to  compute  second  guess  value  of  r2  "  *ng  Equation  (1). 

(e.)  Use  r2  to  compute  new  values  of  w  from  Equation  (13)  and  p,  p  from  Equations  (9)  and  (10) 
which  gives  a  new  value  of  h(p,p)  using  subroutine  TGAS4. 

(f.)  Solve  Equation  (15)  using  the  new  values  of  w,  p,  p,  h  to  get  second  error  estimate,  f(r2). 


(g.)  Compute  new  value  of  Fs  by 


TifiT2)-T2f{T,) 

/(r2)-/(ri) 


(16) 


(h.)  Use  Fs  to  compute  nev/  values  of  w,p,  /),  A  as  in  Step  (e.)  and  solve  Equation  (15)  to  get  f(F3). 

(i.)  If  <  10”^  is  true,  then  the  solution  for  w  is  converged.  If  it  is  false,  then  go  to  Step 

(e.)  and  continue  until  converged. 


3.3  NEW  DENSITY 

The  predictor  step  uses  Euler’s  equations  in  non-conservation  form  to  predict  the  properties  of  u. 
V,  w,  and  p.  For  a  perfect  gas  computation,  the  known  total  enthalpy  constraint  is  used  to  determine  p. 
For  equilibrium  air,  the  density  must  be  obtained  by  numerically  solving  the  differential  equation  for  p. 

To  construct  the  differential  equation,  dp,  the  equation  for  total  enthalpy  is  ^^Titten  as 

2 

Ho  =  h{p,p)+^j  (17) 

where  q  is  the  velocity.  Rearranging  Equation  (17)  for  h(p,p)  gives 

2 

Mp.p)  =  i/o-y  (18) 

Casting  Equation  (18)  in  terms  of  differentials  which  gives 
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and  rearranging  Equation  (19)  for  dp  yields 


dp  = 


-9^9-  (If)  ^ 


dp 


(t), 


(20) 


Using  the  following  finite  difference  expressions  for  the  differentials.  Equation  (20)  can  be  computed 
in  closed  form: 


p 


h{p,p)  -  h{p,p  +  Ap) 

Ap 

(21) 

h{p,p)  -  /i(p+  Ap,p) 

Ap 

(22) 

dp  =  p2~  p\ 

(23) 

dp  =  P2-  Pi 

(24) 

d9  =  92-  9i 

(25) 

Ap  =  .OOlp] 

(26) 

Ap  =  .OOlpi 

(27) 

Equations  (21)  -  (27)  are  substituted  into  Equation  (20)  and  solved  for  the  new  density,  p2  with 
Pi?  Pi?  <li?  P2?  92?  h  given.  Here  the  values  of  h  are  determined  using  the  equilibrium  air  subroutine, 
TGAS4.  This  yields  the  following  expression  for  p2 


71(72  -  71)  + 

001pi,pi)l 

(P2  -  Pi) 

L  .OOlp, 

j 

r*(pi.pi)-*(P).iooipi)] 

[ 

■OOlp,  J 

(28) 


where 


71 


92 


1  t 

+  rr  +  trr 


+  r?  +  u 


2 

'•7 


(29) 
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3.4  OBLIQUE  SHOCKS 

The  inviscid  boundary  conditions  at  a  surface  require  that  the  flow  be  tangent  to  the  surface.  ZEUS 
requires  that  the  velocity  vector  be  turned  through  either  a  shock  or  an  expansion  wave  to  align  the 
velocity  with  the  surface  as  illustrated  in  Figure  3.  To  solve  for  flow  properties  across  an  oblique 
shock  wave  in  equilibrium  air,  a  closed  form  solution  does  not  exist  Refering  to  Figure  3(a),  the 
following  energy  equation  is  solved  iteratively; 

lii  +  2/i_(p_,p-)  =  -H  2/i+(p+,/)+)  (30) 

where  for  an  oblique  shock, 

=  VI  sin  /3 


ujf  =  V-  cos  ^ 


(1  —  tan  0cot  0) 
(tan  0  +  cot  0) 


p+=p-u.!u+  ^  ^ 

p+  =  p-u-{it-  ~  u+)  +  p- 
e  =  cos"^  [(n+  •n_)/|n^||7T_|] 

and  0  is  the  shock  wave  angle  to  be  determined.  Here  the  —  and  +  subscripts  denote  the  upstream  and 
downstream  conditions,  respectively.  V_  is  the  upstream  velocity  component  tangent  to  the  surface, 
w_  and  are  the  velocity  components  normal  to  the  shock,  n_  ,  n+  are  the  surface  normal  vectors, 
p_ ,  are  the  pressures,  and  p— ,  /?+  are  the  densities,  and  0  is  the  deflection  angle. 

To  solve  Equations  (30)  -  (31),  the  Secant  Method  iterative  procedure  is  implemented  in  a  similar 
fashion  as  in  Section  3.2  for  decoding.  Computations  proceed  by  initializing 

sin/?  >  max  [sin  sin(l/A/_)]  (32) 

where  Af_  is  the  upstream  Mach  number.  The  error  is  estimated  using  Equation  (30)  as  follows: 

fis\u0)  =  +  2h.(p^.p_)  -  w".  -  2h  +  {p+.pj,)  (33) 

Tlie  value  of  sm0  is  increased  until  the  relative  error  is  less  than  10“^  or 
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/(sin  ,3) 
^+(P+</’+) 


10“* 


(34) 


This  results  in  the  final  properties  behind  the  shock,  p+,  p+,  for  which  the  velocity  is  tangent  to 
the  surface. 


3.5  EXPANSION  WAVES 


To  satisfy  the  inviscid  boundary  condition  at  a  surface,  often  the  velocity  vector  must  be  turned 
through  a  Prandtl-Meyer  expansion  to  align  the  flow  tangent  to  the  surface.  As  in  the  case  for 
oblique  shock  turns,  a  closed  form  solution  does  not  exists  for  equilibrium  air  conditions;  therefore, 
the  expansion  must  be  computed  iteratively. 


Referring  to  Figure  3(b),  the  equation  for  a  centered  expansion  can  be  written  in  differential  form; 


—  1 


for  0  <  (p  <  6 


(35) 


where 


=  2Hq  -  2h{p.p) 

S-  =  5(p_,/)_) 

p  =  p(5_,p)  (36) 

=  a‘{p-p) 

6  =  cos~^  [(n^.  ■  n_)/|u+j|rT_j]  =  turn  angle 

The  —  and  +  subscripts  represent  upstream  and  downstream  conditions,  respectively.  The  unsubscripted 
variables  are  integrated  values.  Computations  proceed  by  integrating  Equations  (35)  -  (36)  on 
the  interval  0  <  0  <  6  subject  to  the  initial  conditions  p  =  and  q  =  V_  at  (p=0  using  the 
second-order  Improved  Euler  method.  Subdivisions  are  every  .377°.  The  final  results  at  0=6  give 

P+  P(i>=9^  P-\-  P<pz=$'r  1^  q<i=$- 

3.6  SHOCK  ANGLE 


The  computation  of  the  bow  shock  angle  is  determined  using  information  along  the  positive 
characteristic  at  the  cell  adjacent  to  the  bow  shock.  The  shock  angle  is 


3+  =  tan  ' 


0  4- 


(37) 
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or  in  terms  of  velocities 


=  tan 


-1 


' 

W  IVi+w^ 


-  1 


Here,  the  variables  are  Roe  averaged  in  the  following  manner: 

(4v-;+  +  /.!v„-) 


V...  = 


War  = 


pI  +  pI 


P'l+Pl 


cL  =  (Tav  -  1) 


-  \{yL  + 


(38) 


(39) 


Ho 


ot 


(^pIHo+  +  plHo-^ 

pI  + pi 


T 


av 


(4r+  + 

4  +  ^’- 


The  +  and  -  subscripts  represent  propenies  in  the  free  stream  and  at  the  cell  adjacent  to  the  shock, 
respectively.  Vi,  is  the  velocity  compxjnent  normal  to  the  cell  edge;  w  is  the  velocity  component  in  the 
streamwise  direction;  a  is  the  speed  of  sound;  Ho  is  the  total  enthalpy;  and  T  is  the  effective  gamma. 
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CHAPTER  4 

ZEUS  INPUT  FILE  (TAPE  5) 


The  ZEUS  input  file,  as  described  in  Reference  2,  allows  the  user  to: 

1.  Control  the  manner  in  which  ZEUS  integrates  the  Euler  equations; 

2.  Select  boundary  condition  options; 

3.  Specify  the  amount  of  output  generated; 

4.  Choose  the  separation  model; 

5.  Select  the  type  of  initial  data  to  be  used; 

6.  Specify  mesh  size  and  clustering. 

The  set  up  commands  and  options  for  the  current  version  of  ZEUS  include  the  items  outlined  above 
from  Reference  2.  However,  two  more  options  are  included  which  are  inserted  between  items  4  and 
5  above  in  the  input  file.  They  are: 

•  Selection  of  the  gas  type; 

•  Specification  of  the  dimensional  units. 

Moreover,  two  input  parameters  are  added  to  item  1  and  one  parameter  to  item  3  above.  The  complete 
input  file  is  listed  in  Appendix  A.  The  following  sections  give  a  description  of  the  input  variables,  in 
order  of  their  occurrance  in  file  zeuse.da. 


4.1  INTEGRATION  CONTROL 


ZETAEND  —  Terminate  the  calculation  when  ZETA  >  ZETAEND.  The  last  step  size  is  adjusted 
to  place  the  final  integration  plane  at  ZETAEND. 

KEND  —  Terminate  the  calculation  after  executing  KEND  steps.  If  the  initial  step  number  is 
STEPI,  the  termination  step  number  if  STEPI+KEND. 

FCFL  —  Step  size  safety  factor.  The  step  size  is  determined  by  multiplying  the  allowed  step 
size,  based  on  the  CFL  condition,  by  FCFL.  For  values  of  FCFL  >  1,  the  calculation  should  be 
unstable.  A  value  of  .9  is  typically  used.  Program  abort  can  sometimes  be  rectified  by  using  a 
smaller  value.  Also,  oscillations  in  surface  properties  may  be  reduced  by  decreasing  FCFL. 
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XKI  —  Limiting  constant  for  cells  not  adjacent  to  a  surface.  The  recommend  value  for  this 
quantity  is  1,  however,  XKI  may  be  assigned  any  value  between  2  and  0.  Increasing  XKI  promotes 
increased  accuracy  and  sharper  shock  capturing  but  may  lead  to  oscillations  in  the  solution.  If 
XKI  =  0,  the  scheme  reduces  to  the  first  order  Godunov  method. 

lAPR  —  The  Davis  approximate  Riemann  problem  is  applied  if  lAPR  =  0,  while  the  complete 
Riemann  problem  is  solved  if  lAPR  =  1 .  On  most  computing  systems  lAPR  =  0  is  faster,  more 
robust  and  gives  nearly  the  same  results  as  lAPR  =  1.  However,  lAPR  =  0  is  more  dissipative, 
and  differences  between  these  solvers  can  be  noted  in  solutions  featuring  strong  crossflow  shocks. 
For  equilibrium  air  (IGAS  =  1),  lAPR  must  be  set  to  0. 

ISHFLUX  —  The  flux  at  the  cell  edge  next  to  the  bow  shock  is  determined  using  the  full  Riemann 
problem  if  ISHFLUX  =  0  or  the  free  stream  conditions  if  ISHFLUX  =  1.  For  most  cases,  a  value 
of  0  can  be  used;  however,  use  1  if  wiggles  occur  at  the  shock. 

ISHANG  —  The  bow  shock  angle  is  determined  using  information  from  the  full  Riemann  problem 
if  ISHANG  =  0  or  determined  from  Roe  averaged  variables  as  described  in  Section  3.6  if  ISHANG 
=  1.  In  general,  a  value  of  0  can  be  used;  however,  for  equilibrium  air  a  value  of  1  is  recommended. 


4.2  BOUNDARY  CONDITIONS 


NXKE  —  Number  of  edges  at  which  the  default  limiter  setting  will  not  be  applied.  Cells  adjacent 
to  a  surface  use  different  limiter  values  than  the  interior  point  setting  of  XKI.  Increasing  the 
limiter  values  improve  the  solution  accuracy  in  smooth  flow  regions,  but  can  lead  to  surface 
pressure  oscillations  near  shocks  and  fin  edges.  Different  limiter  values  can  be  selected  for  each 
edge  of  each  zone.  If  NXKE  =  0,  default  limiter  values  will  be  applied  at  all  edges.  These  are  2 
on  edges  1  and  3  and  0  on  edges  2  and  4.  The  presumption  here  is  that  a  smooth  flow,  such  as  that 
generated  by  a  tangent  ogive  body,  will  occur  near  edges  1  or  3.  However,  shocks  or  expansions 
may  occur  near  edges  2  and  4  as  a  consequence  of  surface  slope  discontinuities  and  fin  edges.  If 
these  default  values  do  not  fit  the  case  under  consideration,  set  NXKE  to  the  number  of  edges  at 
which  default  values  are  to  be  changed.  List  the  zone  number,  edge  number,  and  limiter  value  for 
each  edge  which  is  to  be  changed  at  the  indicated  position  in  the  input  file. 

NSUR  —  Number  of  edges  which  do  not  feamre  the  default  surface  type.  The  user  has  the  latitude 
of  specifying  the  type  of  surface  which  occurs  on  each  zone  edge.  The  possibilities  are: 

0  —  interior;  do  not  apply  boundary  conditions; 

1  —  solid  surface;  apply  tangent  flow  boundary  conditions; 

2  —  free  stream  boundary;  cell  value  fixed  at  free  stream  conditions; 

3  —  fitted  shock  or  expansion. 

If  NSUR  =  0,  edge  1  will  be  type  1  and  edge  3  will  be  type  1  unless  a  shock  is  being  fitted  (i.e., 
ISHOCK  =1).  In  this  case  edge  3  will  automatically  be  set  to  3.  If  edge  3  is  defined  outside 
the  bow  shock  in  the  free  stream,  set  to  2.  Edges  2  and  4  will  be  set  to  0,  however,  a  surface 
will  be  assumed  to  occur  for  u;  <  s  <  ft.  To  change  these  options,  set  NSUR  to  the  number  of 
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edges  to  be  altered  and  for  each  edge  list  zone  number,  edge  number  and  option  number.  Insert 
this  list  at  the  indicated  position  in  the  input  file. 

DFAC  —  At  locations  where  local  Mach  number  is  too  small  to  allow  the  flow  to  be  turned  parallel 
to  the  wall,  the  turn  angle  is  diminished.  This  is  accomplished  by  computing  the  maximum  turn 
angle  associated  with  this  Mach  number.  The  turn  angle  to  be  used  in  the  calculation  is  calculated 
by  multiplying  this  maximum  angle  by  DFAC.  DFAC  should  be  set  between  0  and  1  and  is 
typically  taken  to  be  .9. 


4.3  OUTPUT  CONTROL 


Tape  6  (output  file)  The  first  and  last  crossflow  planes  are  always  written  to  tape  6. 


BPRINT  —  Print  the  crossflow  plane  if  the  step  number  is  evenly  divisible  by  IPRINT.  If  IPRINT 
=  100,  the  crossflow  plane  will  be  printed  at  step  100,  200,  etc.  The  crossflow  plane  at  the  first 
and  last  step  of  the  calculation  will  also  be  printed. 

NSKIP  —  The  crossflow  plane  output  only  contains  those  N  planes  which  are  evenly  divisible 
by  NSKIP.  Set  NSKIP  to  1  if  all  n  planes  are  to  be  printed. 

MSKIP  —  The  crossflow  plane  output  only  contains  those  M  planes  which  are  evenly  divisible 
by  MSKIP.  Set  MSKIP  to  1  if  all  m  planes  are  to  be  printed. 

ISKIP  —  Print  the  step  size  if  the  step  number  is  evenly  divisible  by  ISKJP.  Set  ISKIP  to  1  if 
the  step  size  is  to  be  printed  after  every  step. 


PLOTZA  (Force  and  Surface  Pressure  Summary) 


JSPPR  —  Following  the  completion  of  the  final  computational  step,  PLOTZA,  which  contains 
a  list  of  forces  and  pressures  at  zone  edges  for  each  step,  is  read,  and  surface  pressures  along 
each  edge  may  be  printed.  The  default  option  is  to  list  pressures  only  along  edge  1  of  each  zone. 
JSPPR  is  the  number  of  edges  at  which  this  output  convention  is  to  be  changed.  For  each  altered 
edge,  list  the  zone  number,  edge  number  and  the  print  option;  1  print,  0  don’t  print  this  edge. 
Insert  this  list  at  the  indicated  postion  in  the  input  file. 

JSPDIS  —  Print  the  location  of  each  grid  point  on  the  edges  specified  by  JSPPR  in  addition  to 
the  surface  pressures  if  JSPDIS  is  1  (does  not  appear  on  PLOTZA).  Do  not  print  if  JSPDIS  is  0. 
KSKIP  —  Write  the  surface  pressures,  forces  and  moments  to  PLOTZA  at  steps  with  numbers 
which  are  evenly  divisible  by  KSKIP. 

NPRT  —  On  edges  2  and  4,  only  list  the  surface  pressures  for  n  planes  less  than  NPRT. 
DELZA  —  Write  surface  pressures,forces  and  moments  to  PLOTZA  at  step  zeta  intervals  of 
DELZA. 


PLOTZe  (Plot  tape)  Crossflow  plane  at  first  and  last  step  is  always  written  to  PLOTZC 
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IPLOT  —  Write  the  crossflow  plane  to  PLOTZC  at  step  numbers  which  are  evenly  divisible  by 
IPLOT. 

DELZC  —  Write  the  crossflow  plane  to  DELZC  at  zeta  intervals  of  DELZC. 

BPLOTN  —  Number  of  target  zeta  stations  at  which  crossflow  plane  will  be  written  to  PLOTZC. 
If  IPLOTN  >  0,  a  list  of  these  stations,  in  ascending  order  should  be  included  as  indicated  in 
the  input  file. 


4.4  AERODYNAMIC  DATA 


AREF  —  Reference  area  to  be  used  in  the  force  and  moment  calculations. 

XLREF —  Reference  length  to  be  used  in  the  force  and  moment  calculations. 
lEFORCE  —  Controls  the  output  of  forces  on  individual  edges.  Application  of  this  option  allows 
forces  on  individual  fins  to  be  determined.  Set  lEFORCE  =  1  to  invoke  this  option.  Otherwise 
let  lEFORCE  =  0. 


4.5  SEPARATION  MODELING 

rVIS  —  Type  of  separation  model: 

0  —  no  separation  modeling; 

1  —  clipping; 

2  —  forced  separation. 

If  options  0  or  1  are  selected,  the  remainder  of  the  information  in  this  section  is  unnecessary.  Note 
that  option  2  can  only  be  used  with  cylindrical  coordinates.  If  option  2  is  used,  the  information  listed 
below  is  needed: 

1  NSEP  —  The  number  of  separation  lines.  Use  of  more  than  one  separation  line  is  necessary  if 
separation  on  a  body  without  pitch  plane  symmetry  is  to  be  modeled.  Additional  separation  lines 
can  also  be  used  to  model  secondary  separation. 

2.  For  each  separation  line  specify: 

•  ISSIDE  —  the  side  of  the  body  on  which  separation  occurs;  0-0°<c>j<180°,  1-  180°  < 
<ps  <  360°,  where  <pj  is  the  location  of  the  separation  line; 

•  ISEP  —  number  of  points  used  to  define  the  separation  line; 

•  ZSSEP  —  Zeta  value  at  which  separation  is  staned; 

•  ZESEP  —  Zeta  value  at  which  separation  is  terminated; 

•  PHICD,  PHIAD,  BETACDjBETAAD  —  Flow  direction  in  degrees.  Typical  values  are  20°, 
20°,  20°,  5°,  respectively. 

•  A  list  of  ISEP  pairs  of  separation  line  coordinates  (ZEPZ,  ZEPP),  where  ZEPZ  and  ZEPP 
are  the  z  and  c>  coordinates  for  each  point. 
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4.6  GAS  TYPE 

IGAS  —  Type  of  gas 

0  —  Air  is  treated  as  a  perfect  gas; 

1  —  Air  is  treated  in  chemical  equilibrium.  Thermodynamic  properties  are  determined  using 
the  subroutines  for  equilibrium  air.  lAPR  must  equal  0.  Use  of  CXDNES  is  not  recommended, 
and  the  flowfield  must  be  initialized  from  uniform  conditions  (IMOD  =  0)  or  using  the  blunt 
body  program,  BLUNT2. 


4.7  UNITS 


lUNITS  —  Dimensional  units  of  problem.  All  input  units  must  be  consistant.  Output  units  are 
the  same. 

0  —  English  units  are  used  for  input  and  output. 

(p  =  Ibf  /  ft**2  ;  p  =  slug  /  ft**3  ;  T  =  °R 
1  —  Metric  units  are  used  for  input  and  output. 

(p  =  N  /  m**2  ;  p  =  kg  /  m**3  ;  T  =  °K  ; 


4.8  INITIAL  PLANE  DATA 

IMOD 

0  —  Use  a  uniform  initial  flow  field  to  stan  the  calculation. 

1  —  Read  the  initial  flow  field  from  file  START.  This  flow  field  can  be  generated  by: 

•  The  RESTART  file  from  another  ZEUS  run; 

•  Program  CONVERT  which  rezones  RESTART  files; 

•  Program  CONES  which  generates  an  approximate  conical  flow  field; 

•  Program  NOSETIP  or  BLUNT2  which  generates  the  starting  data  plane  from  a  blunt  body 
solution  for  flow  over  a  sphere. 


;  x,y,z  =  ft); 
x,y,z  =  m). 


If  IMOD  =1,  the  remainder  of  the  data  in  this  section  is  unnecessary.  Instructions  for  applying 
CONVERT,  CONES  and  NOSETIP  are  given  in  Reference  2.  For  IMOD  =  0,  the  initial  flow  field 
and  mesh  size  are  defined  using  the  following  information: 


IZN  —  Number  of  zones; 

NA  —  Number  of  cells  in  the  direction; 
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MA  —  Number  of  cells  in  the  7/  direction  summed  over  all  zones.  This  information  is  followed 
by  a  list  of  the  number  of  cells  in  the  rj  direction  in  each  zone; 

ZETA  —  z  coordinate  of  initial  data  plane; 

ALPHA  —  Angle  of  attack  in  degrees  (see  Figure  1); 

BETA  —  Yaw  angle  in  degrees  (see  Figure  1); 

XMINF  —  Free  stream  Mach  number; 

PINF  —  Free  stream  pressure.  The  magnitude  of  the  free  stream  pressure  scales  out  of  solution. 
For  example,  if  the  free  stream  pressure  is  doubled,  all  solution  pressures  are  doubled.  Setting 
PINF  =  1  insures  that  solution  pressures  are  scaled  by  PINF.  However,  for  IGAS  =  1,  PINF  must 
be  dimensional; 

DINF  —  Free  stream  density.  The  magnitude  of  the  free  stream  density  also  scales  out  of  the 
solution.  However,  for  IGAS  =  1,  DINF  must  be  dimensional; 

ICORD  —  0,1,2  for  cartesian,  cylindrical  or  elliptical  coordinates. 

lASYM  —  0  -  flow  is  not  symmetric  about  pitch  plane;  1  -  flow  is  symmetric  about  pitch  plane 
symmetry; 

ISHOCK  —  Set  to  1  if  edge  3  is  to  be  fitted  as  a  shock  or  expansion;  otherwise  set  to  zero. 


4.9  MESH  CLUSTERING 

IMESHF  —  0  if  the  mesh  is  to  be  uniform  in  the  ^  direction,  1  for  a  clustered  mesh.  If  clustering 
in  this  direction  is  used,  a  list  of  the  f  clustering  function  values  at  cell  centers  must  be  inserted 
into  the  input  file  as  indicated.  Stan  the  list  at  the  cell  center  adjacent  to  edge  1. 

IMESHG  —  0  if  a  the  mesh  is  to  be  uniform  in  the  t)  direction,  1  for  a  clustered  mesh.  If 
clustering  in  this  direction  is  used,  a  different  clustering  function  may  be  used  in  each  zone. 
Hence  it  is  necessary  to  specify  g  values  at  cell  centers  in  each  zone.  Insen  cell  center  g  values 
into  the  input  file  as  indicated,  starting  at  zone  1  and  ending  at  the  final  zone.  Start  the  list  at 
the  cell  center  adjacent  to  edge  4  . 
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CHAPTER  5 

PRE-PROCESSING  PROGRAMS 


This  section  describes  the  -programs  used  to  generate  the  START  file  for  a  ZEUS.  Programs 
CXDNES,  NOSETIP,  GEOTEST,  and  CONVERT  are  described  in  Reference  2  and  will  not  be  further 
discussed  here. 


5.1  BLUNT2 


The  BLUNT2  program  generates  the  START  tape  for  flow  over  a  sphere  for  air  as  both  a  perfect 
gas  and  in  chemical  equilibrium.  The  program  in  turn  creates  the  ZEUS  START  file.  BLUNT2  can 
be  applied  for  Mach  numbers  greater  than  1.5;  however,  the  freestream  Mach  number  must  be  large 
enough  to  insure  that  the  axial  component  of  the  flow  is  supersonic  at  the  sphere- afterbody  junction 
so  that  the  ZEUS  code  is  marchable. 

The  BLUNT2  INPUT  file  is  prepared  by  editing  the  following  file: 


The  following  file  describes  the  required  input  for  BLUNT2 .  To  execute  BLUNT2, 
replace  the  ()  symbols  in  the  value  column  with  the  desired  variable  value. 

Add  information  supporting  specific  options  between  the  lines  containing 
the  >>>>>  and  <<<<<  symbols.  Dc  not  eliminate  any  lines  from  this  file. 


VALUE  VARIABLE  DESCRIPTION 


Read  in  the  title  of  comiputaton  using  A-format  with  a  maximrum  of  72  characters. 
>>>>>start  title 
<<<<<er,d  title 


() 

0 

0 

o 


*****FREESTREAM  CONDITIONS***** 

AMINE  Mach  number 

PINE  Pressure  used  if  IATMP=0;  otherwise,  set  PINE=1.  For  English 

units  (IUMTS=0)  use  ibf/ft**2,  and  for  m.etric  units  {rJNITS=l) 
use  N /m.*  *2  . 

TINE  Temperature  used  if  IATM?=0;  otherwise  set  TINE=1 .  For  English 
units  (lUNITS-C)  use  Rankine,  and  for  m^etric  units  (IUNITS=1) 
use  Kelvin. 

lAT.MP  C-use  input  values  of  FINE  and  TINE;  59-use  1959  ARDC  Tables; 

or  62-use  1  962  C .  S .  Standard  Tables  to  com.pute  PINE  and  TINE. 
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i'l  Altitude  used  it  IATM?  =  S9  cr  52  to  cortpute  FIN?  and  TIN'; 

otherwise,  if  IATNiP=C,  set  ALT=C.  For  English,  units  {IT'MTS=0) 
use  feet,  and  for  metric  units  (IL!MTS=1)  use  meters. 


*****GAS  AND  UNITS  OPTIONS***** 

O  IRG  0,  i  -  air  is  a  perfect  gas  or  in  chemical  equilibriu.' 

()  lUNITS  0,  1  -  dimensions  are  in  English  or  metric  units. 


»****SPHElMi  DATA***** 

()  RN  Nose  radius  of  the  sphere  in  ft  for  IUNITS=0  or  m  for  IUNITS=1. 

0  YMAX  Maximum  angle  (in  degrees)  over  the  sphere  for  which  the  flow 

field  will  be  computed.  The  0  deg  line  correspo.nds  to  the 
stagnation  stream, line.  (Typically  12C.) 

0  NMAX  Number  of  mesh  points  between  tre  body  and  the  shoc’<. 

(Typically  18) . 

0  MMAX  Number  of  mesh  points  in  the  streamiwise  direction.  (Typically  18) 

0  IMAX  Numbier  of  tire  steps.  Generally  at  least  500  steps  is  needed 

to  co.n verge  -  ’r.e  solution. 


•****ZEUS  DATA*'*** 

0  ZSTART  Axial  distance  from  the  nosetip  to  the  starting  plane  in 

ft  for  IUNITS=C  or  m  for  IUNITS=:. 

()  ZO  Axial  distance  fro.--'  the  nosetip  to  the  center  of  the  sphere  in 

ft  for  IuNITS=0  or  m  for  IUNITS=I 

0  lASYM  0,  1  -  nor.  symmetric  cr  symmetric  ficwfield  computation.  If 

BETA  not  equal  to  0,  then  set  IASYM=0. 

0  NA  Number  of  m.esh  points  between  the  body  and  shock  for  the 

ZEUS  starting  crossflow  plane. 

0  MA  Number  of  mesh  points  in  the  circumiferential  direction  for  the 

ZEUS  starting  crossflow  plane. 

()  ALPHA  Angle  of  attack  (in  degrees). 

()  BETA  Angle  of  yaw  (in  degrees). 
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CHAPTER  6 
SAMPLE  CASES 


This  section  discusses  the  computation  of  sample  cases  using  ZEUS  for  equilibrium  chemically 
reacting  air. 


6.1  HYPERSONIC  SHARP  CONE 


This  case  compares  perfect  gas  and  equilibrium  air  computations  for  a  sharp  cone  at  0°  incidence 
and  an  altitude  of  100,000  ft.  Solutions  were  computed  for  three  different  hypersonic  similarity 
parameters,  Kc  =  A/ocsin^c  =  5.18,  8.45,  14.79.  Calculations  were  performed  using  the  geometry 
described  in  cylindrical  coordinates  on  a  one  zone  18  x  18  mesh.  The  solution  was  initiated  using 
uniform  flow  conditions  at  ZETA  =  .  1  and  terminated  at  ZETA  =  50.  In  addition,  results  are  provided, 
as  well  as  output  found  in  Appendix  B,  for  the  sharp  cone  at  the  same  altitude  for  Moc  =  20  and  20° 
incidence  on  a  18  x  36  single  zone  mesh. 


6.1.1  Geometry 


A  description  of  the  edges  of  the  zone  are  as  follows: 


EDGE  1:  b  =  ( tan  6c.  bz  =  tan  6c  ,  =  0. 

EDGE  2:  =  tt,  =  tv  =  0. 

EDGE  3:  not  needed  (shock  fitted) 

EDGE  4:  CT  =  0.,  ctz  =  ar  =  0. 


The  required  geometry  routines  are  as  follows: 


SUBROUTINE  BCONST (ZETA, BOFZ, BZOFZ) 

C . THIS  SUBROUTINE  CALCULATES  THE  MAJOR  AXIS  OF  AN  ELLIPSE 

C  (BOFZ)  AND  ITS  Z  DERIVATIVE  (BZOFZ) 

RETURN 
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END 

SUBROUTINE  EDGEl (IN, ZETA, YY, BT, BZT, BYT) 

C . THIS  SUBROUTINE  DEFINES  EDGEl 

C . USER  DEFINES  BT,BZT  AND  BYT  AS  A  FUNCTION  OF  ZETA  (ZETA) 

C  Y  OR  PHI  (YY)  AND  ZONE  NUMBER  (IN) . 

C  INPUT:  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  BT  -  X  OR  R  LOCATION  OF  EDGE  1 

C  BZT  -  Z  DERIVATIVE  OF  BT . 

c  by:  derivative  of  bt. 

C . CHOOSE  CONE  ANGLE,  THETAC=1 5 , 2 0 , 25, ETC . 


THETAC=25. 

PI=4 . *ATAN (1 . ) 

RAD=PI/180 . 

THETAR=THETAC*RAD 
BT=ZETA*TAN (THETAR) 

BZT=IAN (THETAR) 

BYT=0 . 

RETURN 
END 

SUBROUTINE  EDGE2 (IN, ZETA, XX, THT, THZT, TKXT) 

THIS  SUBROUTINE  DEFINES  EDGE2 

USER  DEFINES  THT, THZT  AND  THXT  AS  A  FUNCTION  OF  ZETA  (ZETA) 
C  X  OR  R  (XX)  AND  ZONE  NUMBER  (IN) . 

C  INPUT:  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  THT  -  Y  OR  PHI  LOCATION  OF  F'' ' 

C  THZT  -  Z  DERIVATIVE  OF  THT 

C  THXT  -  X  DERII'Ari”"  -F  THT. 

C 

THT=4.*ATAN(1.) 

THZT=0. 

THXT=0 . 

RETURN 

END 

SUBROUTINE  EDGES (IN, ZETA, PHI , CT, CZT,CXT) 


C . USER  DEFINES  CT,CZT,CYT  IF  BOW  SHOCK  IS  NOT  BEING  TRACKED. 

C  Y  OR  PHI  (YY)  AND  ZONE  NUMBER  (IN) . 

C  INPUT:  IN  -  ZONE  HUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  CT  -  X  OR  R  LOCATION  OF  EDGE  3 

C  CZT  -  Z  DERIVATIVE  OF  CT. 

C  CYT  -  Y  DERIVATIVE  OF  CT . 

C 

RETURN 


END 

SUBROUTINE  EDGE4 (IN, ZETA, XX, SGT, SGZT, SGXT) 


C . THIS  SUBROUTINE  DEFINES  EDGE4 

C . THE  USER  MUST  PROVIDE  A  DEFINITION  OF  SGT, SGZT, SGXT,  AS  A  FUNCTION 

C  OF  ZETA  (ZETA),  X  OR  R  (YY)  AND  ZONE  NUMBER  (IN). 

C  INPUT:  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  SGT  -  Y  OR  PHI  LOCATION  OF  EDGE  4 
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C  SGZT  -  Z  DERIVATIVE  OF  SGT . 

C  SGXT  -  X  DERIVATIVE  OF  SGT. 

C 

SGT=0 . 

SGZT=0. 

SGXT-0. 

RETURN 

END 


SUBROUTINE  PEDGE (IN, II, ZETA, POUTL, PINL) 

C  THIS  SUBROUTINE  CALCULATES  THE  FIN  EDGE  LOCATION 

C  REQUIRED  INPUT 

C  IN  -  ZONE  NUMBER 

C  II  -  1  FOR  EDGE  4,  2  FOR  EDGE  2 

C  ZETA  -  ZETA 

C 

POUTL=0 . 

PINL=0. 

C . USER  PROVIDES  DEFINITION  OF  POUTL  AND  PINL 

RETURN 

END 


6.1.2  ZEUS  Input  File  —  I 


The  following  ZEUS  input  file  is  used  for  the  sharp  cone  with  a  25°  half-angle  at  Mach  35  and 
0°  incidence  (note  that  the  flowfield  is  initiated  with  uniform  conditions  and  the  START  file  is  not 
needed): 

The  following  file  describes  the  required  input  for  ZEUS.  To  execute  ."EUS, 
replace  the  ()  syrtbols  in  the  value  colur-n  with  the  desired  variable  value. 

Add  information  supporting  specific  optio.ns  betwee.n  the  lines  containing 
the  >>>>>  and  <<<<<  symbols.  Variables  enclosed  by  [j  should  be  places  by 
themselves  on  a  lire.  Do  not  eliminate  any  lines  from  this  fi^e. 


VALUE  VARIABLE  DESCRIPTION 


50  . 
5000  . 

.9 
1 . 

0 

0 


.....integration  control***** 

ZETAEND  Terminate  calculation  if  ZETA  >  ZET.^END. 

KEND  Terminate  calculation  if  STEP  NUMBER  >  KENT. 

FCFL  Step  size  safety  factor  (1.  >  FCFL  >  C.)  Typically  use  .9  . 

XKI  Interior  point  limiting  ccnsta.nt  (2.  >  XKI  >  C.)  Typica^Ay  use  .. 

lAPR  0,1  for  apWoxim.ate  or  complete  Riemanr,  Problem.  Typi-ally  use  C. 
Must  use  C  for  equilibrium,  air. 

ISHFLUX  0,1  for  using  complete  Riemar.n  Problem  or  freestream  properties 
to  compute  the  flux  at  the  cell  edge  next  tc  the  bow  shock. 
Typically  use  0. 

ISHANG  0,1  for  using  com.plete  Riem.ann  Problem  or  Roe  averaged  variables 
to  compute  si-ock  angle.  Typically  use  C.  Per  equilibrium  air,  a 
value  cf  1  is  preferred. 


■BOUNDARY  CCNDITICNS - 
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0  NXKE  Number  of  edges  at  which  default  limiter  setting  won't  be  used. 

Default  settings  are  2  on  edges  1  or  3  and  0  on  edges  2  or  4  . 
For  each  edge  at  which  the  default  limiter  is  non-standard,  add 
a  line  containing  the  following  information: 
zone  number  (KEZ),  edge  number  (KEE)  and  lim.iter  value  (XKE) 
>>>>>start  list  of  [KEZ, KEE, XKE] 

<<<<<end  list 

0  NSUR  Number  of  edges  not  featuring  default  surface  types.  Default 

values  are  1  on  edges  1  and  3  and  0  on  edge  2  and  4.  Program 
will  automatically  account  a  fitted  shock  on  edge  3  or  fin 
surfaces  on  edges  2  or  4 .  For  each  edge  at  which  the  surface 
type  is  not  of  default,  add  a  line  that  specifies  the  zone 
number (KSURZ) ,  edge  number (KSURE)  and  surface  type  (KSUR) , 
»>>>start  list  of  [KSURZ,  KSURE,  KSUR] 

<<<<<ena  list 

.9  DFAC  In  case  of  a  local  wall  Mach  number  which  is  to  small  to  allow 

the  flow  to  be  turned  parallel  to  the  wall,  the  turn  angle  is 
multiplied  by  the  constant  DFAC  (  1.  >  DFAC  >  0.) . 


250 

IPRINT 

3 

NSKIP 

6 

MSKIP 

10 

ISKIP 

0 

JSPPR 

0 

JSPDIS 

50 

KSKIP 

100 

NPRT 

10000. 

DELZA 

10000  IPLOT 
10000.  DELZC 
0  IPLOTN 


1 .  AREF 

1.  XLREF 


OUTPUT  CONTROL***** 


- Printer - 

Print  crossflow  plane  if  step  number  is  evenly  divided  by  IPRINT 

Print  n  planes  which  are  evenly  divisible  by  NSKIP 

Print  r  planes  which  are  evenly  divisible  by  MSKIP 

Print  step  size  if  step  number  is  evenly  divided  by  ISKIP 

- PLOTZA/Force-Pressure  Summ.ary - 

In  default  mode,  only  edge  1  of  each  zone  will  be  written  to 
PLOTZA.  JSPPR  is  the  number  of  additional  edges  which  are  not 
to  be  written  in  the  default  m.anner.  For  each  edge  added  or 
deleted  from,  this  print  list, include  a  line  whimu  specifies 
zone  number  (JS?Z),  edge  number  (JSPE)  and  print  code  (JS?)  . 

JSP  =  1  will  write  this  zone  edge  to  PLOTZA,  0  will  not. 

>>>>>st;rt  list  of  [JSPZ , JSPE, JSP ; 

<<<<<end  list 

1  -  print  grid  point  locations  for  the  edges  specified  by  CSPPR. 
C  -  do  not  print. 

Write  surface  properties  on  selected  edges  to  PLOTZA  if 
if  STEP  NUMBER  is  evenly  divisible  by  KSKl? 

On  summary  sheet,  print  surface  pressures  of  edges  2  or  if 
n  <  NPST . 

Write  to  PLOTZA  if  (ZETA  -  ZETA  at  last  write)  >  DELZA 
- PLOTZC  (Plot  file) 

Write  to  PLOTZC  if  STEP  NUMBER  is  evenly  divisible  by  IPLOT. 
Write  to  PLOTZC  if  (ZETA  -  ZETA  at  last  write)  >  DELZC 
Number  of  target  Z  stations  (ZTARGET)  at  which  PLClzC  will 
be  written.  If  IPLOTN  >  0,  Include  a  list  of  these 
stations  on  the  next  line  (m.aximum  of  20.  Stations  miUst  be 
listed  in  ascending  order. 

>>>>>start  ZTARGET  list 
<<<<<end  list 


****AERODYNA.MIC  DATA***** 

Reference  area  used  in  calculating  force  and  moment  coefficients. 
Reference  length  used  in  calculating  m.omen.t  coefficients  and 
center-of-pressure. 
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0  lEFORCE  0,  1  don't,  do  print  force  and  moments  for  individual  edges. 


.....separation  modeling***** 

0  IVIS  0,1,2  for  no  modeling,  clipping,  and  forced  separation.  Clipping 

and  forced  separation  can  only  be  used  in  conjunction  with 
cylindrical  coordinates.  If  IVIS  =  2,  specify  the  number  of 
separation  lines  (NSEP). 

>>»>  NSEP 
<<<<< 

For  each  separation  line  specify: 

ISSIDE  -  0  separation  line  located  between  0  and  180  degrees, 

1  separation  line  located  between  180  and  360  degrees. 
ISEP  -  number  of  prints  usee  to  define  the  separation  line. 

ZSSEP  -  Zeta  v.alue  at  which  separation  is  started 
ZESEP  -  Zeta  value  at  which  separation  is  termdnated 
PHICD  -  separation  angle  phiC  (degrees) 

PHIAD  -  separation  angle  phiA  (degrees) 

BETACD  -  sparstion  andgle  BetaC  (degrees) 

BETAAD  -  separation  angle  BetaA  (degrees) 

>>>>>  [ISSIDE, ISEP, ZSSEP, ZESEP. PHICD, PHIAD, BETACD, BETAAD] 

<<<<< 

A  list  of  ISEP  pairs  of  points  (ZEP2,ZE?P)  describing  the  zeta 
and  phi  coordinates  of  the  separatio.n  line.  This  list  should 
start  at  the  most  forward  section  and  move  aft.  Start  the  list 
for  each  separation  line  on  a  separate  line  and  list  in  same 
order  as  used  in  the  above  data. 

>>>>>  start  list  of  ZEPZ,ZEPP 
<<<<< 


*****GAS  TYPE***** 

1  IGAS  0-perfect  gas;  !-equiIibrium  air 


1 


0 


1 

18 

18 


18 


*****D.MTS  ***** 

lUNITS  0-english  (p=lbf /f t* *2,  rho=siug/ft* *3,  t=rankine,  xyz=ft); 

1-metric  (p=n/m**2,  rho=kg/m**3,  t  =  kelvin,  xyz=m,)  . 

All  input  units  must  be  consistant.  Output  units  are  the  same. 


*****in:tial  data***** 

IMOD  0,1  for  new  start  or  restart  respectively.  For  a  restart,  the 
initial  flow  field  is  read  from  TAPES.  Additional  information 
is  not  needed  ana  the  remaining  entries  in  t.hio  section  should 
be  disregarded.  A  new  start  assumes  a  uniform,  flow  field  which 
is  described  in  the  following  section. 

- Mesh  Size - 

IZN  Number  of  zones. 

NA  Number  of  cells  in  xi  direction  (i.e.,  between  edges  i  and  3) 

MA  Number  of  cells  in  eta  direction  (i.e.,  between  edge  4  of  zone  1 

and  edge  2  of  zone  IZN) 

List  the  number  of  M  planes  in  each  zone,  starting  from  zone  i 
and  ending  at  zone  MA. 

>>>>>Start  list  of  MA  for  each  zone 

<<<<<end  list 
- Flow  Properties - 

ZETA  Zeta  ccorainate  of  initial  data  pi,an.e 
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0 .  ALPHA 

0.  BETA 

35.  XMINF 

1185.5  PIN? 
.017861  DINF 


Angle  of  attack  (degrees) 
Angle  of  yaw  (degrees) 
Mach  njmber 
Pressure 
Density 


- Problem  Definition - 

PROBLEM  DEFINITION 

ICORD  0,1,2  for  cartesian,  cylindrical  or  elliptic  coordinates 

lASYM  0  -  no  pitch  plane  symmetry;  1  -  pitch  plane  symmetry. 

ISHOCK  Set  to  1  if  edge  3  is  to  be  fitted  by  the  calculation 
(either  as  a  shock  or  sonic  line);  ot.herwise  set  to  0. 

ISHOCK  =  1  is  only  valid  for  ICORD  =1  or  2.  If  edge  3 
is  to  be  fitted,  the  r  and  phi  (cylindrical  coordinates) 
or  s,  tau  (elliptic  coordinates)  of  the  shock  location  on  each 
M  plane,  (SS,  SSZ)  must  be  defined,  starting  from  M  =  1. 
etc,  . 

>>>>>Start  list  of  (SS.SSZ] 


>>>>>End  list 


**»**MESH  DEFINITION’**** 


IMESHF  0  for  uniform  mesh  in  xi  direction.  I  for  clustered  mesh.  If 
a  clustered  mesh  is  to  be  used,  t.he  clustering  function  (FN)  ; 
the  xi  direction  at  the  NA  cell  centers  must  be  defined 
starting  at  cell  1,  locates  adjacent  to  edge  1  anc  miov.ng  tc 
edge  3  (usually  body  to  shock). 

>>>>>start  list  of  FN 
>>>>>end  list 

IMESHG  0  for  uniform  mesh  in  eta  direction.  1  for  clustered  mesh.  ] 
a  clustered  mesh  is  tc  be  used,  the  clustering  function  (GN)  ; 
the  eta  direction  at  the  MA(i)  cel.^  centers  in  each  zone  mrust 
be  defined,  starting  with  zone  1  and  ending  with  zone  IZN.  Sti 
with  cell  1,  locatec  adjacent  tc  edge  4, and  end  at  the  cell 
adjacent  to  edge  2.  The  list  for  each  zone  should  be  started 
or.  a  new  card 
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6.1.3  ZEUS  Input  File  —  D 

The  following  ZEUS  input  file  is  used  for  the  sharp  cone  with  a  25°  half-angle  at  Mach  20  and 
20°  incidence  (note  that  the  flowfield  is  initiated  with  uniform  conditions  and  the  START  file  is  not 
needed): 


The  following  file  describes  the  required  input  for  ZEUS.  To  execute  ZEUS, 
replace  the  ()  symbols  in  the-  value  column  with  the  desired  variable  value 
Add  information  supporting  specific  options  between  the  lines  containing 
the  >>>>>  and  <<<<<  symbols.  Vaiiabj.es  enclosed  by  []  should  be  placed  by 
themselves  on  a  line.  Do  ret  eliminate  any  lines  from  this  file. 


VALUE  VARIABLE  DESCRIPTION 


•“••INTEGRATION  CONTROL*^*^^ 

ZETAEND  Term, irate  calculation  if  ZETA  >  ZETAEND. 

KEND  Term.inate  calculation  if  STEP  NUMBER  >  KEND . 

FCFL  Step  size  safety  factor  (1.  >  FCFL  >  C.)  Typically  use  .9  . 

XKI  Interior  po.ir.t  lim.itinc  constant  {2.  >  XKI  >  C  . )  Typically  use  1. 

lAPR  0,1  for  approximate  or  com.plete  Riemann  Problem..  Typically  use  C. 
Must  use  0  for  equilibrium,  air. 

ISHFLUX  0, 1  for  using  complete  Riemann  Problem,  or  freestream  properties 
to  compute  the  flux  at  the  cell  edge  next  to  the  bow  shock. 
Typically  use  0. 

ISHANG  0,1  for  using  complete  Riemann  Problem  or  Roe  averaged  variables 
to  com.pute  shock  angle.  Typically  use  0.  For  equilibrium,  air,  a 
value  of  1  is  preferred. 


•••••BOUNDARY  CONDIT IONS • • • • • 

NXKE  Number  of  edges  at  which  default  lim.iter  setting  won't  be  used. 
Default  settings  are  2  on  edges  1  or  3  and  0  on  edges  2  or  <  . 
For  each  edge  at  which  the  default  lim.iter  is  non-standard,  add 
a  line  containing  the  following  information: 
zone  num.ber  (KEZ),  edge  num.ber  (KEE)  and  lim.iter  value  (XKE) 
>>>>>start  list  of  IKE2, KEE,XKE] 

0 

<<<<<end  list 

NSUR  Number  of  edges  not  featuring  default  surface  types.  Default 
values  are  1  on  edges  1  and  3  and  C  on  edge  2  and  <  .  Frograr- 
wiil  autom.at i ca  1 1  y  account  a  fitted  shock  c.n  edge  3  or  fin 
surfaces  on  edges  2  or  <  .  For  each  edge  at  which,  the  surface 
type  is  not  of  default,  add  a  line  that  specifies  the  zone 
number ( KSURZ ) ,  edge  number (KSURE)  and  surface  type  (KSUR). 
>>>>>start  list  of  I  KSURZ , KSURE , KSUR ; 

<<<<<end  list 

DFAC  In  case  of  a  local  wall  Mach  nu.mber  which  is  to  small  to  allow 
the  flow  to  be  turned  parallel  to  the  wall,  the  turn  angle  is 
multiplied  by  the  constant  DFAC  <  1.  >  DFAC  >  C.). 


•••••OUTPUT  CONTRCl.  •  •  *  •  • 

- Printer - 

2hC  IFRINT  Print  crossflow  plane  if  step  number  is  evenly  divided  by  'PRINT 

3  NSKI?  Print  n  planes  which  are  ever.ly  divisible  by  .NS.^1? 
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.9 
1 . 

0 

0 
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Print  m  planes  which  are  evenly  civisirle  by  .VSKIP 
Print  step  size  if  step  nufier  is  ever.-y  diVidea  by  ISKIP 

- PLOTZA/Force-Pressure  Surttary - 

In  default  mode,  only  edge  1  of  each  zone  wil^  be  written  to 
PLOT2A.  JSPPR  is  the  number  of  additronal  edges  which  are  not 
to  be  written  in  the  default  manner.  For  each  edge  added  or 
deleted  from  this  print  list, include  a  line  which  specifies 
zone  number  (JSPZ)  ,  edge  .numbier  <JSFE)  and  print  code  (JSP). 

JSP  =  1  will  write  this  zone  edge  to  PLCTZA,  0  will  not. 

>>>>>start  list  of  ( JSPZ, JSPE, JSP; 

<<<<<end  list 

1  -  print  grid  point  locations  for  the  edges  specified  by  JSPPR. 
0  -  do  not  print. 

Write  surface  properties  on  selected  edges  to  PLCTZA  if 
if  STEP  NUMBER  is  evenly  divisible  by  KSKI? 

On  summary  sheet,  print  surface  pressures  of  edges  2  or  4  if 
n  <  NPRT. 

Wri’.o  uc  i-LClZA  if  (ZETA  -  ZETA  at  last  write)  >  DELZA 
- PLOTZC  (Plot  file) 

Write  to  PLOTZC  if  STEP  NUMBER  is  evenly  divisible  by  IPLOT. 
Write  to  PLOTZC  if  (ZETA  -  ZETA  at  last  write)  >  DELZC 
Number  of  target  Z  stations  (ZTARGET)  at  wnich  PLOTZC  will 
be  written.  If  IPLOTN  >  C,  Include  a  list  of  these 
stations  on  the  next  line  (maximum  of  20).  Stations  must  be 
listed  in  ascending  order. 

>>>>>start  ZTARGET  list 
<<<<<end  list 


•**• 'AERODYNAMIC  DATA*"" 

1.  AREF  Reference  area  used  in  calculating  force  and  moment  coefficients. 

1.  XLREF  Reference  length  used  in  calculating  m.cment  coefficients  and 

center-of-pressure . 

0  lEFO.RCE  0,  1  don't,  do  print  force  and  moments  for  individual  edges. 


""'SEPARATION  MODELING'"" 

C  IVIS  C,I,2  for  no  m.odeling,  clipping,  and  forced  separation.  Clipping 

and  forced  separation  nan  only  be  used  in  conjunction  with 
cylindrical  coordinates.  If  IVIS  =  2,  specify  the  number  of 
separation  lines  (NSEP). 

>>>>s  NSEP 
<<<<< 

For  each  separation  line  specify: 

ISSIDE  -  0  separation,  line  located  cetween  C  and  18C  degrees, 

1  separation  line  located  between  181.  and  36C  degrees, 
ISEP  -  number  of  points  useo  to  define  tne  separation  line. 

ZSSEP  -  Zeta  value  at  which  separation  is  started 
ZESEP  -  Zeta  value  at  which  separation,  is  terminateo 
PHICD  -  separation  angle  phiC  (degrees) 

PHIAD  -  separation,  angle  phiA  (degrees) 

BETACD  -  sparation  an.dgle  BetaC  (degrees) 

BETAAD  -  separation  angle  BetaA  (degrees) 

>>>>>  [ISSIDE,  ISEP, ZSSEP, ZESu?,F.nICD,P.n -  AD, BrTACD,E  r  T  AAD  . 

<<<<< 

A  list  of  ISEP  pairs  of  points  (ZE?Z,ZE??)  describing  the  zeta 
and  phi  coordinates  of  the  separation  line.  Tn.is  list  shculo 
start  at  the  most  fcrwarc  section  and  move  aft.  Start  the  list 


6 

MSKIP 

1C 

ISKIP 

0 

JSPPR 

0 

JSPDIS 

50 

KSKIP 

100 

NPRT 

10000. 

DELZA 

10000 

IPLOT 

10000. 

DELZC 

0 

IPLOTN 
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for  each  separation  line  on  a  separate  line  and  list  in  same 
order  as  used  in  the  above  data. 

>>>>>  start  list  of  ZE?Z,ZEPP 
<<<<< 


.....gas  type***** 

1  IGAS  0-perfect  gas;  l-equilibrium.  air 


»»***UNITS  *•*•* 

1  lUNITS  0-english  (p=lbf/ft**2,  rho=slug/f t • *3,  t=rankine,  xyz=ft); 

1-metric  (p=n/m**2,  rho=kg/m**3,  t=kelvin,  xyz=m) . 

All  input  units  must  be  consistant.  Output  units  are  the  same. 


»«***INITIAL  DATA***** 

0  IMOD  0,1  for  new  start  or  restart  respectively.  For  a  restart,  the 

Initial  flow  field  is  read  from  TAPE3 .  Additional  information 
is  not  needed  and  the  remaining  entries  in  this  section  should 
be  disregarded.  A  new  start  assumes  a  uniform  flow  field  which 
is  described  in  t.he  following  section. 


1  IZN 

18  NA 

36  MA 


36 


.1  ZETA 

20 .  ALPHA 

0.  BETA 

20.  XyiNF 

1185.5  PINF 
.017861  DINF 


1  ICORD 

1  lASYM 

1  ISHOCK 


.30  2.5 
.30  7.5 
.30  12.5 
.30  17.5 
.30  22.5 
.30  27.5 


- Mesh  Size - 

Number  of  zones. 

Num.ber  of  cells  in  xi  direction  (i.e.,  between  edges  1  and  3) 
Number  of  cells  in  eta  direction  (i.e.,  between  edge  4  of  zone  1 
and  edge  2  of  zone  IZN) 

List  the  number  of  M  planes  in  each  zone,  starting  from,  zone  1 
and  ending  at  zone  MA. 

>>>>>Start  list  of  MA  for  each  zone 

<<<<<end  list 

- Flow  Properties - 

Zeta  coordinate  of  initial  data  plane 
Angle  of  attack  (degrees) 

Angle  of  yaw  (degrees) 

Mach  number 
Pressure 
Density 


- Problem  Defir, itio.n - 

PROBLEM  DEFINITION 

0,1,2  fcr  cartesian,  cylindrrcal  or  elliptic  coordinates 
0  -  no  pitch  plane  syiroretry;  1  -  pitch  plane  sym.metry. 

Set  to  1  if  edge  3  is  to  be  fitted  by  the  calculation 
(either  as  a  shocK  or  sonic  line);  otherwise  set  to  C. 

ISHOCK  =  1  is  only  valid  for  ICORD  =1  or  2.  If  edge  3 
is  tu  be  fitted,  the  r  and  phi  (cylindrical  coordinates) 
or  s,  tau  ^(elliptic  coordinates)  of  t.he  shock  location  on  each 
M  plane,  (SS,  SSZ )  must  be  defined,  starting  from  ,M  =  i, 
etc,  . 

>>>>>Start  list  of  ;SS,SSZ; 
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30 

32.1' 

30 

37.5 

30 

ri.b 

30 

4'' .  5 

3C 

52.5 

30 

57.5 

30 

62.5 

30 

67.5 

30 

72.5 

30 

77 . 5 

30 

82.5 

30 

87.5 

30 

92.5 

30 

97.5 

30 

102.5 

30 

107.5 

30 

112.5 

30 

117.5 

30 

122.5 

30 

127 . 5 

30 

132.5 

30 

137.5 

30 

142.5 

30 

147 . 5 

30 

152.5 

30 

157.5 

30 

162.5 

30 

167.5 

30 

172.5 

30 

177.5 

>>>>>End  list 


***»*MESH  DEFINITION****’ 

0  IMESKF  0  for  uniform  mesh  in  xi  direction.  1  for  clustered  mesh.  I 

a  clustered  mesh  is  to  be  used,  the  clustering  functio.*-.  (FN> 
the  xi  direction  at  the  NA  ce*l  ce.nters  must  be  defi.ned 
starting  at  cell  1,  located  adjacent  to  edge  1  a.td  moving  to 
edge  3  (usually  body  to  shock). 

>>>>>start  list  of  FN 
>>>> >end  list 

0  IKESKG  0  for  uniform,  mesh  in  eta  direction.  1  for  clustered  mesh. 

a  clustered  mcsn  -s  to  be  used,  the  cluster ino  function  (GN) 


'.e  eta  direction  at  the  XA  ( i )  cell  centers  in  eaci 


zone  must 


be  defined,  starting  with  zcne  1  and  ending  with  zone  IZN.  Start 


with  ceil 


located  adjacent 


to  ecge  ,and 


adjacent  to  edge  2.  The 
on  a  new  card 

start  list  GN 


i  St  for  each  zone 


ena 

shot 


at  tne  cex- 
id  be  startec 


end  list<<<<<<< 
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6.1.4  Results 


Conical  flow  solutions  have  been  computed  for  a  right  circular  sharp  cone  with  6c  =  25°  and 
for  Mach  20  at  0°  incidence.  Results  are  compared  with  tabulated  cone  properties  by  Hudgins^®  at 
an  altitude  of  100,000  ft.  Figure  4  displays  surface  property  comparisons  between  equilibrium  air 
and  calorically  perfect  results  as  a  function  of  the  hypersonic  similarity  parameter,  IQ.  Computed 
results  are  in  good  agreement  with  tabulated  results  for  values  of  IQ  =  5.18,  8.45,  and  14.79.  The 
characteristics  of  equilibrium  chemically  reacting  air  are  well  predicted  which  show  that  for  chemically 
reacting  air  (1)  the  pressure  is  strongly  dependent  on  the  mechanical  aspects  of  the  flow  and  not  the 
chemical  reactions,  (2)  the  density  ratio  is  higher  implying  a  smaller  shock  layer  thickness,  and  (3) 
the  equilibrium  temperatures  are  lower  than  calorically  perfect  results,  a  result  of  the  dissociated  air. 
Figure  5  illustrates  pressure  contours  for  Mach  20  at  20°  incidence.  Equilibrium  air  results  demonstate 
a  smaller  shock  layer  thickness  as  compared  to  calorically  perfect  results. 


6.2  HYPERSONIC  BLUNT  NOSE  CYLINDER 


This  case  consists  of  a  spherically  blunt  nosed  cylinder  at  Mach  20,  incidences  of  0°  and  20°, 
and  an  altitude  of  20,000  meters  for  equilibrium  air.  The  output  for  20°  is  provided  in  Appendix  B. 
The  calculation,  which  uses  a  cylindrical  coordinate  description  of  the  geometry,  is  divided  into  two 
sections.  The  first  is  the  spherically  blunt  nose  calculation  using  the  program  BLUNT2  on  a  18  x 
25  mesh  which  generates  a  ZEUS  staning  solution  at  ZETA  =  .5.  The  second  computes  the  cylinder 
body  using  a  36  x  36  single  zone  mesh  and  is  terminated  at  ZETA  =  20. 


6.2.1  Geometry 


A  description  of  the  edges  of  the  zone  are  as  follows: 


EDGE  1:  b  ==  .5,  b,  =  ^  0. 

EDGE  2:  xl'  =  t:,  Vz  =  i'r  -  0. 
EDGE  3:  not  needed  (shock  fitted) 
EDGE  4:  a  =  0.,  =  (Tr  -  0. 


The  required  geometry  routines  are  as  follows: 


SUBROUTINE  BCONST ( ZETA , BCFZ , BZCFZ ) 

C . THIS  SUBROUTINE  CALCULATES  THE  MAJOR  AXIS  OF  AN  ELLIPSE 

C  (BOFZ)  AND  ITS  Z  DERIVATIVE  (BZOFZ) 

RETURN 
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END 

SUBROUTINE  EDGEl (IN, ZETA, YY , BT, BZT, BYT) 

C . THIS  SUBROUTINE  DEFINES  EDGEl 

C . USER  DEFINES  BT,BZT  AND  BYT  AS  A  FUNCTION  OF  ZETA  (ZETA) 

C  Y  OR  PHI  (YY)  AND  ZONE  NUMBER  (IN) . 

C  INPUT:  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  BT  -  X  OR  R  LOCATION  OP  EDGE  1 

C  BZT  -  Z  DERIVATIVE  OF  BT. 

C  BYT  -  Y  DERIVATIVE  OF  BT. 


BT-.  5 
BZT-0 . 

BYT-0 . 

RETURN 

END 

SUBROUTINE  EDGE2 (IN, ZETA, XX, THT, THZT, THXT) 


C . THIS  SUBROUTINE  DEFINES  EDGE2 

C . USER  DEFINES  THT, THZT  AND  THXT  AS  A  FUNCTION  OF  ZETA  (ZETA) 

C  X  OR  R  (XX)  AND  ZONE  NUMBER  (IN) . 

C  INPUT:  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  THT  -  Y  OR  PHI  LOCATION  OF  EDGE  2 

C  THZT  -  Z  DE^TVATIVE  OF  THT. 

C  THXT  -  X  ''ERIVATIVE  OF  THT. 


C 

THT=4.*ATAN(1. ) 

THZT«=0. 

THXT-0. 

RETURN 

END 

SUBROUTINE  EDGES (IN, ZETA, PHI , CT, CZT, CXT) 


C . USER  DEFINES  CT,CZT,CYT  IF  BOW  SHOCK  IS  NOT  BEING  TRACKED. 

C  Y  OR  PHI  (YY)  AND  ZONE  NUMBER  (IN) . 

C  INPUT:  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT;  CT  -  X  OR  R  LOCATION  OF  EDGE  3 

C  CZT  -  Z  DERIVATIVE  OF  CT. 

C  CYT  -  Y  DERIVATIVE  OF  CT . 


C 

RETURN 

END 

SUBROUTINE  EDGE4 (IN, ZETA, XX, SGT, SGZT, SGXT) 


C . THIS  SUBROUTINE  DEFINES  EDGE4 

C . THE  USER  MUST  PROVIDE  A  DEFINITION  OF  SGT, SGZT, SGXT ,  AS  A  FUNCTION 

C  OF  ZETA  (ZETA),  X  OR  R  (YY)  AND  ZONE  NUMBER  (IN). 

C  INPUT;  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  SGT  -  Y  OR  PHI  LOCATION  OF  EDGE  4 

C  SGZT  -  Z  DERIVATIVE  OF  SGT. 

C  SGXT  -  X  DERIVATIVE  OF  SGT. 

C 

SGT=0. 

SGZT=0 . 
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SGXT-O. 

RETURN 

END 

SUBROUTINE  PEDGE (IN, II, ZETA, POUTL, PINL) 

C  THIS  SUBROUTINE  CALCULATES  THE  FIN  EDGE  LOCATION 

C  REQUIRED  INPUT 

C  IN  -  ZONE  NUMBER 

C  II  -  1  FOR  EDGE  4,  2  FOR  EDGE  2 

C  ZETA  -  ZETA 

C 

POUTL-0 . 

PINL-0. 

C . USER  PROVIDES  DEFINITION  OF  POUTL  AND  PINL 

RETURN 

END 


6.2.2  START  File 


The  start  file  for  20°  incidence  is  generated  using  the  BLUNT2  program  using  the  following  input: 


The  following  file  describes  the  required  input  for  BLUNT2 .  To  execute  BL'JNTZ, 
replace  the  ()  symbols  in  the  value  column  with  the  desired  variable  value. 

Add  information  supporting  specific  options  between  the  lines  containing 
the  >>>>>  and  <<<<<  symbols.  Do  not  eliminate  any  lines  from  this  file. 


VALUE  VARIABLE  DESCRIPTION 


.. » * •title* • *  * • 

Bead  in  the  title  of  computaton  using  A-format  with  a  m.axim.um  of  72  characters. 
>>>>>start  title 

Hypersonic  Blunt  Nose  Cylinder,  M=20,  Alpha=0,  Alt=2C,000  m 
<<<<<end  title 


20.  AMINE 

1.  PINE 

1.  TINE 

59  lATMP 

2000C.  ALT 


•••••EREESTREAM  CONDITIONS***** 

Mach  number 

Pressure  used  if  IATM?=C;  otherwise,  set  PINE=1.  Eor  English 
units  (IUNITS=0)  use  ibf/ft**2,  and  for  metric  u.nits  {IUNITS=1) 
use  N/m*  *2 . 

Temperature  used  if  IATMP=0;  otherwise  set  TINE=1.  Eor  English 
units  (IUNITS=C)  use  Rankine,  and  for  metric  units  (:UNITS=1) 
use  Kelvin. 

0-use  input  values  of  PINE  and  TINE;  59-use  1959  ARDC  Tables; 
or  62-use  1962  U.S.  Standard  Tables  to  com.pute  PINE  and  TINE. 
Altitude  used  if  IATMP=59  or  62  to  compute  PINE  and  TINE; 
otherwise,  if  IATMP=0,  set  ALT=C ,  Eor  English  units  {IUNITS=C) 
use  feet,  and  for  metric  units  (IUNITS=1)  use  meters. 


*****GAS  AND  UNITS  OPTIONS***** 

1  IRC  0,  1  -  air  is  a  perfect  gas  or  in  chemical  equilibrium.. 

1  lUNITS  C,  1  -  dime.nsions  are  in  English  or  .'"etric  units. 
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.5 

120. 

RN 

YMAX 

18. 

NMAX 

25. 

200C 

MMAX 

I  MAX 

.5 

Z START 

.5 

ZO 

1 

lASYM 

36 

NA 

36 

MA 

20. 

ALPHA 

0. 

BETA 

.....sphere  data***** 

Nose  radius  of  the  sphere  in  ft  for  IUNITS=0  or  tr.  for  IUNITS=1  . 
Maximum  angle  (in  degrees)  over  the  sphere  for  which  the  flow 
field  will  be  com.puted.  The  0  deg  line  corresponds  to  the 
stagnation  streamline.  (Typically  120.) 

Number  of  mesh  points  between  the  body  and  the  shoc)«  . 

(Typically  18)  . 

Number  of  mesh  points  in  the  streamwise  direction.  (Typically  18) 
Number  of  time  steps.  Generally  at  least  500  steps  is  needed 
to  converge  the  solution. 


•****ZE'JS  DATA***** 

Axial  distance  from  the  nosetip  to  the  starting  plane  in 
ft  for  IUNITS=0  or  m  for  IUNITS=i. 

Axial  distance  from  the  nosetip  to  the  center  of  the  sphere  in 
ft  for  IUNITS=C  or  m  for  lUNITS=i. 

0,  1  -  non  symmetric  or  sym.metric  flowfield  computation.  If 
BETA  not  equal  to  C,  then  set  IASYM=0. 

Number  of  mesh  points  between  the  body  and  shoc)^  for  the 
ZEUS  starting  crossflow  plane. 

Number  of  mesh  points  in  the  circumferential  direction  for  the 
ZEUS  starting  crossflow  plane. 

Angie  of  attac)c  (in  degrees). 

Angle  of  yaw  (in  degrees). 


6.2.3  ZEUS  Input  File 


The  following  ZEUS  input  file  is  used  for  the  afterbody: 


The  following  file  describes  the  required  input  for  ZEUS.  To  execute  ZEUS, 
replace  the  ()  symbols  in  the  value  column  with  the  desired  variable  value 
Add  informiation  supporting  specific  options  between  the  lines  containing 
the  >>>>>  and  <<<<<  sym.bols.  Variables  enclosed  by  []  should  be  placed  by 
them.selves  on  a  line.  Do  not  eliminate  any  lines  from  this  file. 


VALUE  VARIABLE  DESCRIPTION 


*****integrat:on  control***** 


10  . 

ZETAEND 

Terminate  calculation  i 

f  ZETA  >  ZETAEND. 

5000  . 

KEND 

Term.inate  calculation  i 

f  STEP  NUMBER  >  KEND. 

.9 

FCFL 

Step  size  safety  factor 

( : .  >  FCFL  >  0 . )  Typ 

ica  1 

1 y  u  s  e  .9  . 

1  . 

XKI 

Interior  point  lim,iting 

constant  (2.  >  XKI  > 

C.) 

Typically  use  1 

C  lAPR  C,1  for  approxim.ate  or  complete  Riem.anr.  Froclem,.  Typically  use  C. 

Must  use  0  for  equilibrium  air. 

1  ISHFLUX  C,  1  for  using  complete  Riem.ann  Problem  or  freestrea"  properties 

to  compute  the  fl„x  at  tn.e  cell  eoge  next  to  the  bow  shock. 
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Typically  use  0. 

1  ISHANG  0,1  for  using  complete  Riemann  Problem  or  Roe  averaged  variables 

to  compute  shock  angle.  Typically  use  0.  For  equilibrium  air,  a 
value  of  1  is  preferred. 


••‘••BOUNDARY  CONDIT IONS • • • • • 

1  NXKE  Number  of  edges  at  which  default  limiter  setting  won't  be  used. 

Default  settings  are  2  on  edges  1  or  3  and  0  on  edges  2  or  4 . 
For  each  edge  at  wr.icr.  the  default  limiter  is  non-standard,  add 
a  line  containing  the  following  information; 

zone  number  (KEZl,  edge  nu.mber  (KEE)  and  limiter  value  (XKE) 
»>>>start  list  of  lKE2,  KEE,  XKE] 

1  1  1. 

<<<<<end  list 

0  NSUR  Number  of  edges  not  featuring  default  surface  types.  Default 

values  are  1  on  edges  1  and  3  and  C  on  edge  2  and  4 .  Program 
will  automatically  account  a  fitted  shock  on  edge  3  or  fin 
surfaces  on  edges  2  or  ^ .  For  each  edge  at  which  the  surface 
type  is  not  of  default,  add  a  line  that  specifies  the  zone 
number (KSURZ) ,  edge  number (KSURE)  and  surface  type  (KSUR) . 
>>>>>start  list  of  [KSURZ, KSURE, KSUR] 

<<<<<end  list 

.9  DFAC  In  case  of  a  local  wall  Kach  number  which  is  to  sm.all  to  allow 

the  flow  to  be  turned  parallel  to  the  wall,  the  turn  angle  is 
multiplied  by  the  constant  DFAC  (  1.  >  DFAC  >  C.). 


•••••OUTPUT  CONTROL***** 


- Pri 

nte  r - 

25C 

IPRINT 

Print 

crossflow  plane  if 

6 

NSKIP 

Print 

n  planes  which  are 

6 

MSKIP 

Print 

m  planes  which  are 

10 

ISKIP 

Print 

step  size  if  step 

step  number  is  eve.nly  aivided  by  IPRINT 
evenly  divisible  by  NSKI? 
evenly  divisible  by  MSKIP 
number  is  evenly  divided  by  ISKIF 


- PLOTZA/Force-Pressure  Summary - 

0  JSPPR  In  de fault]  mode ,  only  edge  1  of  each  zone  will  be  wri.tten  to 

PLOTZA.  JSPPR  is  the  number  of  additional  edaes  wnich  are  .not 
to  be  written  in  tne  default  manner.  For  each  edge  added  or 
deleted  from;  this  print  list,  include  a  ^ine  which,  specifies 
zone  number  (JSPZ),  edge  number  (JSPE)  anc  print  code  (JS?)  . 
JS?  =  1  will  write  this  zone  edge  to  PLCTZA,  C  will  not. 

>>>>>start  list  of  i JSPZ , JSPE , JS? ; 

<<<<<er;d  list 
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0  JSPDIS 

40  KSKIP 

100  NPRT 

10000.  DELZA 


1  -  print,  grid  poi.nt  locatio.ns  for  t.he  edges  specified  by  JSPPR. 
0  -  do  not  print. 

Write  s'jrface  properties  on  selected  edges  to  PLOTZA  if 
if  STEP  NUMBER  is  evenly  divisible  by  KSKIP 

On  summary  sheet,  print  surface  pressures  of  edges  2  or  4  if 
n  <  NPRT. 

Write  to  PLOTZA  if  (ZETA  -  ZETA  at  last  write)  >  DELZA 


10000  IPLOT 
10000.  DELZC 
2  IPLOTN 


- PLOTZC  (Plot  file) 

Write  to  PLOTZC  if  STEP  NUMBER  is  evenly  divisible  by  IPLOT. 
Write  to  PLOTZC  if  (ZETA  -  ZETA  at  last  write)  >  DELZC 
Number  of  target  Z  stations  (ZTARGET)  at  which  PLOTZC  will 
be  written.  If  IPLOTN  >  C,  Include  a  list  of  these 
stations  on  the  next  line  (maximum  of  20.  Stations  must  be 
listed  in  ascending  order. 

>>>>>start  ZTARGET  list 


.5  10. 


<<<<<end  list 


***'*ae.rodykam:c  data*'*** 

1.  AREF  Reference  area  used  in  calculating  force  and  moment  coefficients. 

1.  XLREF  Reference  length  used  in  calculating  m.oment  coefficients  and 

cent er-of -pres su re . 

0  lEFORCE  C,  1  don't,  do  print  force  and  moments  for  individual  edges. 


•'•••SEPARATION  MODELING* ••• • 

0  IVIS  0,1,2  for  no  mode.ing,  clipping,  and  forced  separation.  Clipping 

and  forced  separation  can  only  be  used  in  conjunction  with, 
cylindrical  coordinates.^  If  IVIS  =  2,  specify  the  num.ber  of 
separation  lines  (NSE?)  . 

>>>>>  NSEP 
<<<<< 

For  each  separation  line  specify: 

ISSIDE  -  C  separation  line  locatec  between.  C  ano  18T  degrees, 

1  separation  line  located  between  18C  and  360  degrees. 
ISEP  -  number  of  points  used  to  define  the  separation  line. 

ZSSE?  -  Zeta  value  at  which  separation  is  startec 
ZESE?  -  Zeta  value  at  which  separation  is  terminated 
PHICD  -  separation  angle  phiC  (degrees) 

FHIAD  -  separation  angle  phiA  (aegrees) 


BETAC3 

“  sparatior. 

andcle 

BetaC 

(aegrees) 

BETAAD 

-  separation 

angle 

BetaA 

(degrees) 

>>>>> 

;iss:de, :se? 

, ZSSEP, 

,  ZESEP, 

?r.ZZ2,  ?Hi 
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<<«< 

A  list  of  ISEP  pairs  of  points  (2EPZ,ZEPP)  describing  the  zeta 
and  phi  coordinates  of  the  separation  line.  This  list  should 
start  at  the  most  forward  section  and  move  aft.  Start  the  list 
for  each  separation  line  on  a  separate  line  and  list  in  same 
order  as  used  in  the  above  data. 

>>>>>  start  list  of  ZEPZ.ZEPP 
<<<<< 


*****GAS  TYPE***** 

1  IGAS  0-perfect  gas;  1 -equi librium  air 


*****UNITS  *»*•* 

1  lUNITS  0-english  {p=lbf / f t *  *  2 ,  rho=slug/ft* *3,  t=ran)cine,  xyz  =  ft); 

1-metric  (p=n/m**2,  rho=kg/m**3,  t=kelvin,  xyz=m). 

All  input  units  must  be  consistant.  Output  units  are  the  same. 


*****INIT1AL  DATA***** 

1  IMOD  0,1  for  new  start  or  restart  respectively.  For  a  restart,  the 

initial  flow  field  is  read  from  TAPE3 .  Additional  information 
is  not  needed  and  the  remaining  entries  in  this  section  should 
be  disregarded.  A  new  start  assumes  a  uniform  flow  field  which 
is  described  in  the  following  section. 


- Mesh 

Size  — 

IZN 

Number  of 

zones . 

NA 

Number  of 

cells 

in 

xi  direction  (i.e.. 

between  edges  1 

and  3 ) 

MA 

Number  of 

cells 

in 

eta  direction  (i.e. 

,  between  edge  4 

of  zf.ne  1 

and  edge  2  of  zone  IZN) 


List  the  number  of  M  planes  in  each  zone,  starting  from,  zone  1 
and  ending  at  zone  MA. 

>>>>>Start  list  of  MA  for  each  zone 
<<<<<end  list 


- Flow  Properties - 


() 

ZETA 

Zett  '■ocrdinate  of  initial 

0 

ALP.HA 

Angle  of  attack  (degrees) 

{) 

BETA 

Angle  of  yaw  (degrees) 

0 

XMINF 

Mach  number 

0 

PINF 

Pressure 

0 

DINF 

Density 
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0 

0 

0 


- Problem  Definition - 

PROBLEM  DEFINITION 

ICORD  0,1,2  for  cartesian,  cylindrical  or  elliptic  coordinates 

lASYM  0  -  no  pitch  plane  symmetry;  1  -  pitch  plane  symmetry. 

ISHOCK  Set  to  1  if  edge  3  is  to  be  fitted  by  the  calculation 
(either  as  a  shock  or  sonic  line);  otherwise  set  to  0. 

ISHOCK  ==  1  is  only  v/alid  for  ICORD  =1  or  2.  If  edge  3 
is  to  be  fitted,  the  r  and  phi  (cylindrical  coordinates) 
or  s,  tau  (elliptic  coordinates)  of  the  shock  location  on  each 
M  plane,  (SS,  SSZ)  must  be  defined,  starting  from  M  =  1, 
etc,  . 

»»>Sf.art  list  of  [SS,SSZ] 

>>>>>Eri(i  list 


***»*MESK  DEFINITION***** 

0  IMESHF  0  for  uniform  m.esh  in  xi  direction.  1  for  clustered  mesh.  If 

a  clustered  mesh  is  to  be  used,  the  clustering  function  (FN)  in 
the  xi  direction  at  the  NA  cell  centers  m.ust  be  defined 
starting  at  cell  1,  located  adjacent  to  edge  1  and  moving  to 
edge  3  (usually  body  to  shock) . 

>>>>>start  list  of  FN 
>>>>>end  list 

0  IMESHG  0  for  uniform  mesh  in  eta  direction.  1  for  clustered  mesh.  If 

a  clustered  mesh  is  to  be  used,  the  clustering  function  (GN)  in 
the  eta  direction  at  the  MA(i)  cell  centers  in  each  zone  must 
be  defined,  starting  with  zone  1  and  ending  with  zone  IZN.  Start 
with  cell  1,  located  adjacent  to  edge  <,and  end  at  the  cell 
adjace.nt  to  edge  2.  The  list  for  each  zone  should  be  started 
on  a  new  card 

>>>>>>>>start  list  GN 

ero  1 '-St.<-<<<<<< 


6.2.4  Results 


Computations  were  performed  for  a  spherically  blunt  nosed  cylinder  at  an  altitude  of  20  km  for 
Mach  20  at  0°  and  20°  incidences.  ZEUS  computations  are  compared  with  computed  results  of 
Palmer*^  for  0°  incidence.  Density  and  temperature  contour  plots  in  Figures  6  and  7  compare  well  for 
both  calorically  perfect  and  equilibrium  air.  ZEUS  predicts  the  exhibited  higher  densities  and  thinner 
shock  layer  and  lower  temperatures  associated  with  the  chemically  reacting  air.  Figure  8  illustrates 
density  contours  for  the  20°  incidence  case  at  an  axial  station  9.5  diameters  downstream. 
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6.3  HYPERSONIC  AIRCRAFT 


The  hypersonic  aircraft  configuration'^  of  Figure  9  was  computed  at  Mach  24.5  at  1°  angle  of 
incidence.  A  uniform  36  x  18  single  zone  mesh  was  implemented  over  the  nose  cone.  For  the  wing 
section,  the  flow  field  was  rezoned  into  two  zones  with  each  zone  containing  a  72  x  36  uniform  mesh. 
Here,  the  outer  boundary,  edge  3,  was  defined  outside  the  bow  shock  in  the  free  stream. 


6.3.1  Geometry  —  Nose  Section 


A  description  of  the  edges  of  the  zone  are  as  follows: 


EDGE  1:  6  =  =  tan^c  ’  =  0- 

EDGE  2:  i/i  =  tt,  =  V’r  =  0. 

EDGE  3:  not  needed  (shock  fitted) 

EDGE  4:  (7  =  0.,  at  =  Or  =  0. 


The  required  geometry  routines  are  as  follows: 


SUBROUTINE  BCONST (ZETA, BOFZ, BZOFZ) 

C . THIS  SUBROUTINE  CALCULATES  THE  MAJOR  AXIS  OF  AN  ELLIPSE 

C  (BOFZ)  AND  ITS  Z  DERIVATIVE  (BZOFZ) 

RETURN 

END 

SUBROUTINE  EDGEl (IN, ZETA, YY, BT, BZT, BYT) 

C . THIS  SUBROUTINE  DEFINES  EDGEl 

C . USER  DEFINES  BT,BZT  AND  BYT  AS  A  FUNCTION  OF  ZETA  (ZETA) 

C  Y  OR  PHI  (YY)  AND  ZONE  NUMBER  (IN). 

C  INPUT:  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  BT  -  X  OR  R  LOCATION  OF  EDGE  1 

C  BZT  -  Z  DERIVATIVE  OF  BT. 

C  BYT  -  Y  DERIVATIVE  OF  BT . 


BYT=0 . 

BZT=0 . 

BT-5. 08333 

IF ( ZETA . LE . 63 . 0 ) THEN 
THETAC=4 . 6 
PI=4.*ATAN(1.) 
RAD-PI/180 . 
THETAR=THETAC*RAD 
BT-ZETA*TAN (THETAR) 
BZT=TAN(THETAR) 

BYT-0 . 

ENDIF 
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RETURN 

END 

SUBROUTINE  EDGE2 (IN, ZETA,XX, THT, THZT, THXT) 
C  ....THIS  SUBROUTINE  DEFINES  EDGE2 


C . USER  DEFINES  THT, THZT  AND  THXT  AS  A  FUNCTION  OF  ZETA  (ZETA) 

C  X  OR  R  (XX)  AND  ZONE  NUMBER  (IN) . 

C  INPUT:  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  THT  -  Y  OR  PHI  LOCATION  OF  EDGE  2 

C  THZT  -  Z  DERIVATIVE  OF  THT. 

C  THXT  -  X  DERIVATIVE  OF  THT. 


C 

THT-=4.*ATAN(1.) 

THXT-0. 

THZT-0. 

RETURN 

END 

SUBROUTINE  EDGES (IN, ZETA, PHI , CT , CZT, CXT) 


C . USER  DEFINES  CT,CZT,CYT  IF  BOW  SHOCK  IS  NOT  BEING  TRACKED. 

C  Y  OR  PHI  (YY)  AND  ZONE  NUMBER  (IN) . 

C  INPUT:  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  CT  -  X  OR  R  LOCATION  OF  EDGE  3 

C  CZT  -’Z  DERIVATIVE  OF  CT. 

C  CYT  -  Y  DERIVATIVE  OF  CT. 


C 


RETURN 

END 

SUBROUTINE  EDGE4 (IN, ZETA, XX, SGT, SGZT, SGXT) 

C . THIS  SUBROUTINE  DEFINES  EDGE4 

C . THE  USER  MUST  PROVIDE  A  DEFINITION  OF  SGT, SGZT, SGXT,  AS  A  FUNCTION 

C  OF  ZETA  (ZETA),  X  OR  R  (YY)  AND  ZONE  NUMBER  (IN). 

C  INPUT:  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  SGT  -  Y  OR  PHI  LOCATION  OF  EDGE  4 

C  SGZT  -  Z  DERIVATIVE  OF  SGT. 

C  SGXT  -  X  DERIVATIVE  OF  SGT. 

C 

SGT=0 . 

SGXT=0 . 

SGZT=0. 

RETURN 


END 

SUBROUTINE  PEDGE (IN, II , ZETA, POUTL, PINL) 


C  THIS  SUBROUTINE  CALCULATES  THE  FIN  EDGE  LOCATION 

C  REQUIRED  INPUT 

C  IN  -  ZONE  NUMBER 

C  II  -  1  FOR  EDGE  4,  2  FOR  EDGE  2 

C  ZETA  -  ZETA 

C 

C . USER  PROVIDES  DEFINITION  OF  POUTL  AND  PINL 

POUTL=0 . 

PINL=0. 
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RETURN 

END 


6.3.2  Zeus  Input  File  -  -  Nose  Section 

The  following  input  file  is  used  for  the  hypersonic  aircraft  nose  section  run  at  Mach  24.5  and  1° 
incidence.  Computations  were  initiated  from  uniform  free  stream  conditions. 

The  following  file  describes  the  required  input  for  ZEUS.  To  execute  ZEUS, 
replace  the  0  symbols  in  the  value  colujrin  with  the  desired  variable  value. 

Add  information  supporting  specific  options  between  the  lines  containing 
the  »»>  and  ««<  symbols.  Variable  enclosed  by  []  should  be  placed  by 
themselves  on  a  line.  Do  not  eliminate  any  lines  from  this  file. 


VALUE  VARIABLE  DESCRIPTION 

****»INTEGRATION  CONTROL***** 

ZETAEND  Terminate  calculation  if  ZETA  >  ZETAEND. 

KEND  Terminate  calculation  if  STEP  NUMBER  >  KEND. 

FCFL  Step  size  safety  factor  (1.  >  FCFL  >  0.)  Typically  use  .9  . 

XKI  Interior  point  limiting  constant  (2.  >  XKI  >  0.)  Typically  use  1. 

lAPR  0,1  for  approximate  or  complete  Riemann  Problem.  Typically  use  0. 
Must  use  0  for  equilibrium  air. 

ISHFLUX  0,1  for  using  complete  Riemann  Problem  or  freestream  properties 
to  compute  the  flux  at  the  cell  edge  next  to  the  bow  shock. 
Typically  use  0. 

ISHANG  0,1  for  using  complete  Riemann  Problem  or  Roe  averaged  variables 
to  compute  shock  angle.  Typically  use  0.  For  equilibrium  air,  a 
value  of  1  is  preferred. 

*****B0UNDARY  CONDITIONS***** 

0  NXKE  Number  of  edges  at  which  default  limter  setting  won't  be  used. 

Default  settings  are  2  on  edges  1  or  3  and  0  on  edges  2  or  4 . 

For  each  edge  at  which  the  default  limiter  is  non-standard,  add 
a  line  containing  the  following  information: 

zone  number  (KEZ) ,  edge  number  (KEE)  and  limiter  value  (XKE) 
»»>start  list  of  [KEZ,  KEE,  XKE] 

«<«end  list 

0  NSUR  Number  o^  edges  not  featuring  default  surface  types.  Default 

values  are  1  on  edges  1  and  3  and  0  on  edge  2  and  4 .  Program 

will  automatically  account  a  fitted  shock  on  edge  3  or  fin 
surfaces  on  edges  2  or  4 .  For  each  edge  at  which  the  surface 
type  is  not  of  default,  add  a  line  that  specifies  the  zone 
number (KSURZ) ,  edge  number (KSURE)  and  surface  type  (KSUR) . 
»»>s:;art  list  of  [KSURZ,  KSURE,  KSUR] 

<<<<<end  list 

In  case  of  a  local  wall  Mach  number  which  is  to  small  to  allow 
the  flow  to  be  turned  parallel  to  the  wall,  the  turn  angle  is 
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multiplied  by  the  constant  DFAC  (  1.  >  DFAC  >  0.). 


*** ♦★OUTPUT  CONTROL***** 

- Printer - 

250  IPRINT  Print  crossflow  plane  if  step  number  is  evenly  divided  by  IPRINT 

1  NSKIP  Print  n  planes  which  are  evenly  divisible  by  NSKIP 

3  MSKIP  Print  m  planes  which  are  evenly  divisible  by  MSKIP 

1  ISKIP  Print  step  size  if  step  number  is  evenly  divided  by  ISKIP 

- PLOTZA/Force-Pressure  Summary - 

0  JSPPR  In  default  mode,  only  edge  1  of  each  zone  will  be  written  to 

PLOTZA.  JSPPR  is  the  number  of  additional  edges  which  are  not 
to  be  written  in  the  default  manner.  For  each  edge  added  or 
deleted  from  this  print  list, include  a  line  which  specifies 
zone  number  (JSPZ) ,  edge  number  (JSPE)  and  print  code  (JSP) . 

JSP  =  1  Will  write  this  zone  edge  tc  PLOTZA,  0  will  not. 

»»>start  list  of  [JSPZ,  JSPE,  JSP] 

«<<<end  list 

0  JSPDIS  1  -  print  grid  point  locations  for  the  edges  specified  by  JSPPR. 

0  -  do  not  print. 

5  KSKIP  Write  surface  properties  on  selected  edges  to  PLOTZA  if 

if  STEP  NUMBER  is  evenly  divisible  by  KSKIP 
100  NPRT  On  summary  sheet,  print  surface  pressures  of  edges  2  or  4  if 
n  <  NPRT. 

1000.  DELZA  Write  to  PLOTZA  if  (ZETA  -  ZETA  at  last  write)  >  DELZA 

- PLOTZC  (Plot  file) 

1000  IPLOT  Write  to  PLOTZC  if  STEP  NUMBER  is  evenly  divisible  by  IPLOT. 

1000.  DELZC  Write  to  PLOTZC  if  (ZETA  -  ZETA  at  last  write)  >  DELZC 

2  IPLOTN  Number  of  target  2  stations  (ZTARGET)  at  which  PLOTZC  will 

be  written.  If  IPLOTN  >  0,  Include  a  list  of  these 
stations  on  the  next  line  (maximum  of  20) .  Stations  must  be 
listed  in  ascending  order. 

»>»start  ZTARGET  list 

21.33  4b. 

<<<<<end  list 


*****AERODYNAMIC  DATA***** 

1.  AREF  Reference  area  used  in  calculating  force  and  moment  coefficients. 

1.  XLREF  Reference  length  used  in  calculating  moment  coefficients  and 

center-of-pressure . 

0  lEFORCE  0,  1  don't,  do  print  force  and  moments  for  individual  edges. 


*****SEPARATION  MODELING***** 

0  IVIS  0,1,2  for  no  modeling,  clipping,  and  forced  separation.  Clipping 

and  forced  separation  can  only  be  used  in  conjunction  with 
cylindrical  coordinates.  If  IVIS  =  2,  specify  the  number  of 
separation  lines  (NSEP) . 

>>»>  NSEP 
«<<< 

For  each  separation  line  specify; 

ISSIDE  -  0  separation  line  located  between  0  and  180  degrees, 
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1  separation  line  located  between  180  and  360  degrees. 
ISEP  -  number  of  points  used  to  define  the  separation  line. 
2SSEP  -  Zeta  value  at  which  separation  is  started 
ZESEP  -  Zeta  value  at  which  separation  is  terminated 
PHICD  -  separation  angle  phic  (degrees) 

PHIAD  -  separation  angle  phiA  (degrees) 

BETACD  -  sparation  andgle  BetaC  (degrees) 

BETAAD  -  separation  angle  BetaA  (degrees) 

»»>  [ISSIDE,  ISEP,  2SSEP,  ZESEP,  PHICD,  PHIAD,  BETACD,  BETAAD] 

««< 

A  list  of  ISEP  pairs  of  points  (ZEPZ,ZEPP)  describing  the  zeta 
and  phi  coordinates  of  the  separation  line.  This  list  should 
start  at  the  most  forward  section  and  move  aft.  Start  the  list 
for  each  separation  line  on  a  separate  line  and  list  in  same 
order  as  used  in  the  above  data. 

»»>  start  list  of  ZEP2,2EPP 
««< 


*****GAS  type***** 

1  IGAS  0-perfect  gas;  1 -equilibrium  air 


**»**UNiTS  •***• 

0  lUNITS  0-english  (p=lbf /ft**2,  rho*slug/ft**3,  t=ran)cine,  xyz*^ft)  ; 

1-metric  (p=n/m**2,  rho-kg/m**3,  t=)celvin,  xyz-m)  . 

All  input  units  must  be  consistent.  Output  units  are  the  same. 


0 


1 

36 

18 


*****INITIAL  DATA***** 

IMOD  0,1  for  new  start  or  restart  respectively.  For  a  restart,  the 
initial  flow  field  is  read  from  TAPE3 .  Additional  information 
is  not  needed  and  the  remaining  entries  in  this  section  should 
be  disregarded.  A  new  start  assumes  a  uniform  flow  field  which 
is  described  in  the  following  section. 


- Mesh  Size - 


IZN 

NA 

Number  of 

Number  of 

zones  . 

cells  in 

XI  direction 

(i .  e . , 

between  edges 

1 

and  3) 

MA 

Number  of 

cells  in 

eta  direction 

(i .  e . 

,  between  edge 

4 

of  zone  1 

and  edge 

2  of  zone 

IZN) 

List  the  number  of  M  planes  in  each  zone,  starting  from  zone  1 
and  ending  at  zone  MA. 

>>>>>Start  list  of  MA  for  each  zone 

<<<<<end  list 


- Flow  Properties-- 


.01 

ZETA 

Zeta  coordinate 

of  initial  data  plane 

1 . 

ALPHA 

Angle  of  attack 

(degrees) 

0. 

BETA 

Angle  of  yaw  (degrees) 

24.5 

XMINF 

Mach  number 

1 . 543 

PINF 

Pressure 

1 . 91356 

1-06  , 

DINF  Density 
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- Problem  Definition - 

PROBLEM  DEFINITION 

i  ICORD  0,1,2  for  cartesian,  cylindrical  or  elliptic  coordinates 

J  lASYM  0  -  no  pitch  plane  symmetry;  1  -  pitch  plane  symmetry. 

1  ISHOCK  Set  to  1  if  edge  3  is  to  be  fitted  by  the  calculation 

(either  as  a  shock  or  sonic  line).-  otherwise  set  to  0. 

ISHOCK  “  1  is  only  valid  for  ICORD  ®1  or  2 .  If  edge  3 
is  to  be  fitted,  the  r  and  phi  (cylindrical  coordinates) 
or  s,  tau  (elliptic  coordinates)  of  the  shock  location  on  each 
M  plane,  (SS,  SSZ)  must  be  defined,  starting  from  M  =  1, 
etc,  . 

>>>>>Start  list  of  [SS,SSZ] 


.01 

5. 

01 

15. 

.  01 

25. 

.01 

35. 

.  01 

45. 

.  01 

55  , 

.  01 

65. 

.  01 

75. 

.  01 

85. 

.01 

95. 

.  01 

105. 

.01 

115. 

.  01 

125. 

.01 

135. 

.  01 

145. 

.01 

155. 

.01 

165. 

.  01 

175. 

>>>>>End  list 


*****MESH  DEFINITION***** 

0  IMESHF  0  for  uniform  mesh  in  xi  direction.  1  for  clustered  mesh.  If 

a  clustered  mesh  is  to  be  used,  the  clustering  function  (FN)  in 
the  xi  direction  at  the  NA  cell  centers  must  be  defined 
starting  at  cell  1,  located  adjacent  to  edge  1  and  moving  to 
edge  3  (usually  body  to  shock) . 

>>>>>start  list  of  FN 
>>>>>end  list 

0  IMESHG  0  for  uniform  mesh  in  eta  direction.  1  for  clustered  mesh.  If 

a  clustered  mesh  is  to  be  used,  the  clustering  function  (GN)  in 
the  eta  direction  at  the  MA(i)  cell  centers  in  each  zone  must 
be  defined,  starting  with  zone  1  and  ending  with  zone  IZN.  Start 
with  cell  1,  located  adjacent  to  edge  4, and  end  at  the  cell 
adjacent  to  edge  2.  The  list  for  each  zone  should  be  started 
on  a  new  card 

>»>>»>start  list  GN 

end  list<<<<<<< 
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6.3.3  Geometry  —  >^^ing  Section 


This  section  of  the  hypersonic  aircraft  is  run  using  a  two  zone  model  with  NA=72,  MAZ(1)=36 
and  MAZ(2)=36  mesh.  The  geometry  is  as  follows: 

ZONE  1,  EDGE  1:  6  =  Cfan^ci  bj  =  tan  0c  .  =  0. 

ZONE  1,  EDGE  2:  top  of  the  wing  surface,  OH^ET  =  ttH. 

ZONE  1,  EDGE  3:  outer  boundary  defined  outside  the  bow  shock  in  the  free  stream. 

ZONE  1,  EDGE  4:  leeward  pilch  plane  symmetry. 

ZONE  2,  EDGE  1;  6  =  C  tan^c,  =  tan  0c  >  =  0. 

ZONE  2,  EDGE  2:  windward  pitch  plane  symmetry. 

ZONE  2,  EDGE  3:  outer  boundary  defined  outside  the  bow  shock  in  the  free  stream. 

ZONE  2,  EDGE  4:  bottom  of  the  wing  surface,  OFFSET  =  t/I. 


The  top  of  the  wing  surface  is  described  in  terms  of  bj,b]^  and  parameters  for  a  flat  fin  surface 
(described  in  Reference  2)  as  follows: 


lo  =  -2.048642 

^LE  =  ^0  -  (C  —  33.75)  tan  1° 

VLE  =  (C  -  33.75)  tan  12° 

Zfu  —  22'o  ^  le 

Vru  =  0 

—  i^LE  —  ^ru) !  {y  LE  —  Vtu) 
bj  —  Xtxi 
bj^  =  tan  1° 

Similarly,  for  the  bottom  of  the  wing  surface. 


o-o  =  -2.048642 
^ LE  —  ^0  —  (C  —  33.75)  tan  1° 

Xm  =  0 

6/  =  xle 

b]^  =  —  tan  1° 


Knowledge  of  bj.  bj^  and  x^,  allows  the  appropriate  a  and  v  functions  to  be  calculated  using 
coding  described  in  Reference  2.  The  wing  edges  are  computed  as  follows: 
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Inner  Edge  : 

o 

II 

3 

Outer  Edge  : 

n  =  0 

if  z  <  48.717 

^  =  yl^LE  +  vh 

otherwise 

The  required  geometry  routines  for  the  wing  section  are 

as  follows; 

SUBROUTINE  BCONST (2ETA, BOF2, BZOFZ) 

THIS  SUBROUTINE  CALCULATES  THE  MAJOR  AXIS  OF  AN  ELLIPSE 
(BOFZ)  AND  ITS  2  DERIVATIVE  (BZOFZ) 

RETURN 

END 

SUBROUTINE  EDGEl ( IN, ZETA, YY, BT, BZT, BYT) 


C . THIS  SUBROUTINE  DEFINES  EDGEl 

C . USER  DEFINES  BT,BZT  AND  BYT  AS  A  FUNCTION  OF  ZETA  (ZETA) 

C  Y  OR  PHI  (YY)  AND  ZONE  NUMBER  (IN) . 

C  INPUT:  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  BT  -  X  OR  R  LOCATION  OF  EDGE  1 

C  B2T  -  Z  DERIVATIVE  OF  BT . 

C  BYT  -  Y  DERIVATIVE  OF  BT . 


BYT=0. 

BZT=0. 

BT=5. 08333 

IF(ZETA.LE.63.0)THEN 
THETAC=4 . 6 
PI-4 . *ATAN(1 . ) 

RAD-PI/180 . 

THETAR=THETAC*RAD 
BT=ZETA*TAN (THETAR) 

BZT-TAN (THETAR) 

BYT-0 . 

ENDIF 

RETURN 

END 

SUBROUTINE  EDGE2 (IN, ZETA, XX, THT, THZT, THXT) 


C . THIS  SUBROUTINE  DEFINES  EDGE2 

C . USER  DEFINES  THT, THZT  AND  THXT  AS  A  FUNCTION  OF  ZETA  (ZETA) 

C  X  OR  R  (XX)  AND  ZONE  NUMBER  (IN) . 

C  INPUT:  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  THT  -  Y  OR  PHI  LOCATION  OF  EDGE  2 

C  THZT  -  Z  DERIVATIVE  OF  THT. 

C  THXT  -  X  DERIVATIVE  OF  THT. 

C 


PI=4.*ATAN(1.) 
RAD-PI/180 . 
IF(IN.EQ.1)THEN 
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0FFSET=PI/2 . 

TAN12-TAN(12 . *RAD) 

TAN1=TAN (1 . *RAD) 

Y0=0. 

XO— 2.048642 
YLE- (ZETA-33. 75) *TA'12 
XLE-XO- (ZETA-33. 75) *TAN1 
RLE-SQRT (XLE*XLE+YLE*YLE) 

XBWI-XLE 
XRU“2.*X0-XBWI 
YRU-0 . 

IF (XX . LE . RLE .AND . ZETA . GE . 48 . 717 ) THEN 
XM= (XLE-XRU) / (YLE-YRU) 

BI=XRU 

BIZ=TAN1 

ELI?E 

XM- (XLE-XO) /  (YLE-YO) 

BI-XO 

BIZ-0. 

ENDIF 

ELSE  IF ( IN. EQ. 2) THEN 
OFFSET-PI 
XM-0. 

BI-O  . 

BIZ-0. 

ENDIF 

XM2P1-XM*XM+1 . 

SINT-(BI/XX+XM*SQRT (XM2P1- (BI/XX) **2) ) /XM2P1 
COST-SQRT (1 . -SINT*SINT) 

Cl=COST+XM»SINT 
THT-OFFSET-ASIN (SINT) 

THXT-BI/ (XX*XX*C1) 

THZT--BIZ/ (XX*C1) 

RETURN 

END 

SUBROUTINE  EDGE3 (IN, ZETA, PHI , CT, CZT, CXT) 


C . USER  DEFINES  CT,CZT,CYT  IF  BOW  SHOCK  IS  NOT  BEING  TRACKED. 

C  Y  OR  PHI  (YY)  AND  ZONE  NUMBER  (IN). 

C  INPUT:  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  CT  -  X  OR  R  LOCATION  OF  EDGE  3 

C  CZT  -  Z  DERIVATIVE  OF  CT. 

C  CYT  -  Y  DERIVATIVE  OF  CT. 


PT=4.*ATAN(1. ) 
RAD-PI/180. 

BT-5. 08333 

IF ( ZETA. LE. 63.0) THEN 
TAN4P6-TAN(4 . 6*RAD) 
BT=ZETA*TAN4P6 
ENDIF 
CFACl-1 . 3 
CT1-BT*CFAC1 
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CZT-TAN4P6*CFAC1 
CXT=0 . 

CT2=0. 

IF (ZETA.GT. 48. 717) THEN 
CFAC2-1 . 1 
TAN12=TAN (12 . *RAD) 

TAN1-TAN(1.*RAD) 

XO— 2.048642 
YLE«(ZETA-33.75) *TAN12 
XLE-XO- (ZETA-33 .75) *TAN1 
RLE«SQRT  {XLE*XLE+yiiE*YLE) 
CT2=RLE*CFAC2 

CZT2'=CFAC2*  (-TAN1*XLE+TAN12*YLE)  /RLE 
ENDIF 

CT-AMAXl (CTl , CT2 ) 

IF (CT2 . GE . CTl ) CZT=CZT2 

RETURN 

END 


SUBROUTINE  EDGE4 (IN, ZETA,XX, SGT, SGZT, SGXT) 

C . THIS  SUBROUTINE  DEFINES  EDGE4 

C . THE  USER  MUST  PROVIDE  A  DEFINITION  OF  SGT, SGZT, SGXT,  AS  A  FUNCTION 

C  OF  ZETA  (ZETA),  X  OR  R  (YY)  AND  ZONE  NUMBER  (IN). 

C  INPUT:  IN  -  ZONE  NUMBER 

C  ZETA  -  ZETA 

C  OUTPUT:  SGT  -  Y  OR  PHI  LOCATION  OF  EDGE  4 

C  SGZT  -  Z  DERIVATIVE  OF  SGT. 

C  SGXT  -  X  DERIVATIVE  OF  SGT. 

C 


PI=4.*ATAN(l.) 

RAD-P1/180. 

IFdN.EQ.DTHEN 
OFFSET=0 . 

XM=0  . 

BI=0. 

BIZ-0. 

ELSE  IF (IN. EQ. 2) THEN 
OFFSET-PI/2 . 

TAN12-TAN (12 . *RAD) 

TANl-TAN (1 . *RAD) 

YO-0. 

XO— 2.048642 
YLE- (ZETA-33 . 75) *TAN12 
XLE-XO- (ZETA-33. 75) *TAN1 
RLE-SQRT (XLE*XLE+yLE*YLE) 

IF (XX . LE . RLE . AND . ZETA . GE . 4  8 . 7 1 7 ) THEN 
XM-0  . 

BI-XLE 
BIZ— TANl 
ELSE 

XM- (XLE-XO) / (YLE-YO) 

BI-XO 
BIZ-0 . 

ENDIF 
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ENDIF 

XM2P1-XM*XM+1 . 

SINT=(BI/XX+XM*SQRT(XM2P1- (BI/XX) **2) ) /XM2P1 
COST=SQRT (1 . -SINT*SINT) 

Cl-COST+XM^SINT 

SGT»OFFSET-ASIN(SINT) 

SGXT=BI/ (XX*XX*C1) 

SGZT— BIZ/  (XX*C1) 

RETURN 

END 

SUBROUTINE  PEDGEdN,  II,  ZETA,  POUTL,  PINL) 


C  THIS  SUBROUTINE  CALCULATES  THE  FIN  EDGE  LOCATION 

C  REQUIRED  INPUT 

C  IN  -  ZONE  NUMBER 

C  II  -  1  FOR  EDGE  4,  2  FOR  EDGE  2 

C  ZETA  -  ZETA 

C 

C . USER  PROVIDES  DEFINITION  OF  POUTL  AND  PINL 


POUTL=0 . 

PINL=0. 

IF (ZETA.GE. 48 . 717) THEN 

IF ( (IN.EQ. 1 . AND. II .EQ. 2) .OR. (IN . EQ . 2 . AND . I I . EQ . 1 ) ) THEN 
PI=4 . *ATAN (1 . ) 

RAD=PI/180 . 

TAN12=TAN (12 . *RAD) 

TAN1-TAN(1.*RAD) 

X0=-2. 048642 
YLE-(ZETA-33.75) •TAN12 
XLE=X0- (ZETA-33.75) *TAN1 
POUTL=SQRT (XLE*XLE+YLE*YLE) 

ENDIF 

ENDIF 

RETURN 

END 


6.3.4  CONVERT  Input  File 

The  flow  field  is  rezoned  from  a  one  zone  36  x  18  mesh  into  two  zones  each  with  a  72  x  36  mesh 
using  the  program  CONVERT.  The  input  file  for  CONVERT  is  as  follows: 


INPUT  file  for  CONVERT.  To  execute,  append  the  geometry  subroutines 
BCONST,  EDGEl,  EDGE2 ,  EDGE3,  EDGE4,  and  PEDGE  to  CONVERT  and  execute 
using  this  file  as  input,  with  ()  replaced  by  the  desired  values  and 
required  lists  inserted  between  lines  containing  >>»>  and  «<<<.  Do 
not  delete  or  add  any  lines  to  this  file. 


VALUE  VARIABLE  DESCRIPTION 

1  lASYMN  symmetry  of  the  final  flow  field;  1  -  symmetric 
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12 

72 


36 

0 

0 


0 


IZN>J 

NAN 

MAN 


36 


ISHOCKN 

IMESHF 


IMESHG 


0  -  asyiruretric 
number  of  new  zones 

Number  of  cells  between  edge  1  and  3  of  new  mesh. 
Number  of  cells  between  edge  4,  zone  1  and  edge2, 
zone  IZNN.  The  number  of  cells  in  each  zone  between 
edges  4  and  2  must  also  be  specified,  starting  at 
zone  1  . 

»»>Start  MAZ(l),  MAZ(2),  ...  MAZ(IZZN) 

««<End  list 

0  -  Edge  3  will  not  be  fitted  with  a  shock 
1  -  Shock  fitted  on  edge  3 

0  -  uniform  mesh  in  xi  direction  of  new  mesh. 

1  -  clustered  mesh.  If  a  clustered  mesh  is  used, 
the  clustering  (FNN)  must  be  specified,  starting 
at  the  cell  center  adjacent  to  edge  1  and  ending 
at  the  cell  center  adjacent  to  edge  3. 

»»>Stait  list  FNN 
««<End  list  of  FNN 

0  -  uniform  mesh  between  edges  2  and  4  of  new  mesh. 

1  -  clustered  mesh.  If  a  clustered  mesh  is  used, 
the  clustering  (GNN)  for  each  zone  must  be 
specified,  starting  with  zone  1  and  ending  at  zone 
IZN.  The  clustering  for  each  zone  should  be  started 
on  a  new  card.  Beginning  this  list  at  the  cell 
center  adjacent  to  edge  4  and  terminate  it  at  the 
cell  center  next  to  edge  2. 

»»>Start  list  GNN 


««<End  list  of  GNN 


6.3.5  Zeus  Input  File  —  Wing  Section 


The  following  ZEUS  input  file  is  used  for  the  wing  section: 


The  following  fi.te  describes  the  required  input  for  ZEUS.  To  execute  ZEUS, 
replace  the  ()  symbols  in  the  value  column  with  the  desired  variable  value. 
Add  information  supporting  specific  options  between  the  lines  containing 
the  »»>  and  «<<<  symbols.  Variable  enclosed  by  []  should  be  placed  by 
themselves  on  a  line.  Do  not  eliminate  any  lines  from  this  file. 


VALUE  VARIABLE  DESCRIPTION 


*****imxegration  control***** 

114.25  ZETAEND  Terminate  calculation  if  ZETA  >  ZETAEND. 


5000. 

KEND 

Terminate  calculation  if 

STEP  NUMBER 

> 

KEND  . 

.9 

FCFL 

Step  size  safety  factor 

(1 .  >  FCFL  > 

0. 

)  Typically  use  .9  . 

1 . 

XKI 

Interior  point  limiting 

consta.nt  (2. 

> 

XKI  >  0.)  Typically  use 

1 

0 

lAPR 

0,1  for  approximate  or  complete  Riemann 

Pioblem.  Typically  use 

0 
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0 

1 


0 


2 


1  3 

2  3 

.9 


200 

1 

6 

1 

2 


1  2 
2  4 

0 

1 


Must  use  0  for  equilibrium  air. 

ISHFLUX  0,1  for  using  complete  Riemann  Problem  or  freestream  properties 
to  compute  the  flux  at  the  cell  edge  next  to  the  bow  shock. 
Typically  use  0. 

ISHANG  0,1  for  using  complete  Riemann  Problem  or  Roe  averaged  variables 
to  compute  shock  angle.  Typically  use  0.  For  equilibrium  air,  a 
value  of  1  is  preferred. 


*****BOUNDARy  CONDITIONS***** 

NXKE  Number  of  edges  at  which  default  limiter  setting  won't  be  used. 
Default  settings  are  2  on  edges  1  or  3  and  0  on  edges  2  or  4 . 
For  each  edge  at  which  the  default  limiter  is  non-standard,  add 
a  line  containing  the  following  information; 

zone  number  (KEZ) ,  edge  number  (KEE)  and  limiter  value  (XKE) 
»»>start  list  of  [KEZ,  KEE,  XKE] 

««<end  list 

NSUR  Number  of  edges  not  featuring  default  surface  types.  Default 
values  are  1  on  edges  1  and  3  and  0  on  edge  2  and  4 .  Program 
will  automatically  account  a  fitted  shock  on  edge  3  or  fin 
surfaces  on  edges  2  or  4 .  For  each  edge  at  which  the  surface 
type  is  not  of  default,  add  a  line  that  specifies  the  zone 
number (KSURZ ) ,  edge  number (KSURE)  and  surface  type  (KSUR) . 
»»>start  list  of  [KSURZ,  KSURE,  KSUR] 

2 

2 

««<end  list 

DFAC  In  case  of  a  local  wall  Mach  number  which  is  to  small  to  allow 
the  flow  to  be  turned  parallel  to  the  wall,  the  turn  angle  is 
multiplied  by  the  constant  DFAC  (  1.  >  DFAC  >  0.). 


*.***OUTPUT  CONTROL***** 


- Printer - 

IPRINT  Print  crossflow  plane  if  step  number  is  evenly  divided  by  IPRINT 
NSKIP  Print  n  planes  which  are  evenly  divisible  by  NSKIP 

MSKIP  Print  m  planes  which  are  evenly  divisible  by  MSKIP 

ISKIP  Print  step  size  if  step  number  is  evenly  divided  by  ISKIP 


- PLOTZA/Force-Pressure  Summary - 

JSPPR  In  default  mode,  only  edge  1  of  each  zone  will  be  written  to 

PLOTZA.  JSPPR  is  the  number  of  additional  edges  which  are  not 
to  be  written  in  the  default  manner.  For  each  edge  added  or 
deleted  from  this  print  list, include  a  line  which  specifies 
zone  number  (JSPZ) ,  edge  number  (JSPE)  and  print  code  (JSP) . 

JSP  =  1  will  write  this  zone  edge  to  PLOTZA,  0  will  not. 

>>>>>start  list  of  [JSPZ, JSPE, JSP ] 

1 

1 

<«<<end  list 

JSPDIS  1  -  print  grid  point  locations  for  the  edges  specified  by  JSPPR. 
0  -  do  not  print . 

KSKIP  Write  surface  properties  on  selected  edges  to  PLOTZA  if 
if  STEP  NUMBER  is  evenly  divisible  by  KSKIP 
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100  NPRT  On  summary  sheet,  print  surface  pressures  of  edges  2  or  4  if 
n  <  NPRT. 

1000.  DELZA  Write  to  PLOTZA  if  (ZETA  -  ZETA  at  last  write)  >  DELZA 


- PLOTZC  (Plot  file) 

1000  IPLOT  Write  to  PLOTZC  if  STEP  NUMBER  is  evenly  divisible  by  IPLOT. 

1000.  DELZC  Write  to  PLOTZC  if  (ZETA  -  ZETA  at  last  write)  >  DELZC 

2  IPLOTN  Number  of  target  Z  stations  (ZTARGET)  at  which  PLOTZC  will 

be  written.  If  IPLOTN  >  0,  Include  a  list  of  these 
stations  on  the  next  line  (maximiom  of  20)  .  Stations  must  be 
listed  in  ascending  order. 

»»>start  ZTARGET  list 

63.917  108.67 

««<end  list 


**** ‘AERODYNAMIC  DATA***** 

1.  AREF  Reference  area  used  in  calculating  force  and  moment  coefficients. 

1.  XLREF  Reference  length  used  in  calculating  moment  coefficients  and 

center-of-pressure . 

0  lEFORCE  0,  1  don't,  do  print  force  and  moments  for  individual  edges. 


*****SEPARATI0N  MODELING***** 

0  IVIS  0,1,2  for.no  modeling,  clipping,  and  forced  separation.  Clipping 

and  forced  separation  can  only  be  used  in  conjunction  with 
cylindrical  coordinates.  If  IVIS  =  2,  specify  the  number  of 
separation  lines  (NSEP) . 

»»>  NSEP 

««< 

For  each  separation  line  specify: 

ISSIDE  -  0  separation  line  located  between  0  and  180  degrees, 

1  separation  line  located  between  180  and  360  degrees. 
ISEP  -  number  of  points  used  to  define  the  separation  line. 

ZSSEP  -  Zeta  value  at  which  separation  is  started 
ZESEP  -  Zeta  value  at  which  separation  is  terminated 
PHICD  -  separation  angle  phiC  (degrees) 

PHIAD  -  separation  angle  phiA  (degrees) 

BETACD  -  sparation  andgle  BetaC  (degrees) 

BETAAD  -  separation  angle  BetaA  (degrees) 

»»>  [ISSIDE,  ISEP,  ZSSEP,  ZESEP,  PHICD,  PHIAD,  BETACD ,  BETAAD ] 

«<« 

A  list  of  ISEP  pairs  of  points  (ZEPZ,ZEPP)  describing  the  zeta 
and  phi  coordinates  of  the  separation  line.  This  list  should 
start  at  the  most  forward  section  and  move  aft.  Start  the  list 
for  each  separation  line  on  a  separate  line  and  list  in  same 
order  as  used  in  the  above  data. 

»»>  start  list  of  ZEPZ,ZEPP 
«<<< 


*****GAS  TYPE***** 

1  IGAS  0-perfect  gas;  1-equilibrium  air 

*****UN1TS  ***** 

0  lUNlTS  0-english  (p=lbf /ft**2,  rho=slug/ft**3,  t  =  ran)cine,  xyz  =  ft)  ; 
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1 


0 

() 

0 


0 

0 

0 

0 

0 

0 


0 

0 

() 


*  ★  * 

0 


0 


it  it 

IMOD 


IZN 

NA 

MA 


ZETA 

ALPHA 

BETA 

XMINF 

PINF 

DINF 


I  CORD 
lASYM 
I SHOCK 


1-metric  (p=n/in**2,  rho»kg/m**3,  t-Jtelvin,  xyz-m)  . 

All  input  units  must  be  consistent .  Output  units  are  the  same. 


♦♦‘INITIAL  DATA***** 

0,1  for  new  start  or  restart  respectively.  For  a  restart,  the 
initial  flow  field  is  read  from  TAPES.  Additional  information 
is  not  needed  and  the  remaining  entries  in  this  section  should 
be  disregarded.  A  new  start  assumes  a  uniform  flow  field  which 
is  described  in  the  following  section. 


- Mesh  Size - 

Number  of  zones. 

Number  of  cells  in  xi  direction  (i.e.,  between  edges  1  and  3) 
Number  of  cells  in  eta  direction  (i.e.,  between  edge  4  of  zone  1 
and  edge  2  of  zone  IZN) 

List  the  number  of  M  planes  in  each  zone,  starting  from  zone  1 
and  ending  at  zone  MA. 

»»>Start  list  of  MA  for  each  zone 
<<<<<end  list 

- Flow  Properties - 

Zeta  coordinate  of  initial  data  plane 
Angle  of  attack  (degrees) 

Angle  of  yaw  (degrees) 

Mach  number 

Pressure 

Density 

- Problem  Definition - 

PROBLEIM  DEFINITION 

0,1,2  for  cartesian,  cylindrical  or  elliptic  coordinates 
0  -  no  pitch  plane  symmetry;  1  -  pitch  plane  symmetry. 

Set  to  1  if  edge  3  is  to  be  fitted  by  the  calculation 
(either  as  a  shock  or  sonic  line);  otherwise  set  to  0. 

ISHOCK  =  1  is  only  valid  for  ICORD  =1  or  2.  If  edge  3 
is  to  be  fitted,  the  r  and  phi  (cylindrical  coordinates) 
or  s,  tau  (elliptic  coordinates)  of  the  shock  location  on  each 
M  plane,  (SS,  SSZ)  must  be  defined,  starting  from  M  =  1, 
etc,  . 

»»>Start  list  of  [SS,SSZ] 

>>>>>End  list 


**MESH  DEFINITION* *’\** 

IMESHF  0  for  uniform  mesh  in  xi  direction.  1  for  clustered  mesh.  If 
a  clustered  mesh  is  to  be  used,  the  clustering  function  (FN)  in 
the  xi  direction  at  the  NA  cell  centers  must  be  defined 
starting  at  cell  1,  located  adjacent  to  edge  1  and  moving  to 
edge  3  (usually  body  to  shock) . 

»>»start  list  of  FN 
>>>»end  list 

IMESHG  0  for  uniform  mesh  in  eta  direction.  1  for  clustered  mesh.  If 
a  clustered  mesh  is  to  be  used,  the  clustering  function  (GN)  in 
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the  eta  direction  at  the  MA{i)  cell  centers  in  each  zone  must 
be  defined,  starting  with  zone  1  and  ending  with  zone  IZN.  Start 
with  cell  1,  located  adjacent  to  edge  4, and  end  at  the  cell 
adjacent  to  edge  2.  The  list  for  each  zone  should  be  started 
on  a  new  card 

»»»»start  list  GN 

end  list«««< 


6.3.6  Results 


Computations  were  performed  for  the  hypersonic  airplane  configuration  of  Figure  9  at  Mach  24.5 
and  1°  incidence.  A  comparison  of  pressure  profiles  between  ZEUS  and  the  Parabolized  Navier- 
Stokes  computations  of  Korte*^  on  the  windward  symmetry  plane  at  z  =  108.67  feet  shows  good 
agreement  except  near  the  surface  where  viscous  effects  are  not  predicted  by  ZEUS.  In  Figure  10, 
pressure  contours  are  also  computed  at  an  axial  station  near  the  end  of  the  body  (z  =  108.67  feet)  for 
equilibrium  air.  Computations  illustrate  good  agreement  except  in  the  vicinity  of  the  wing  tip  where 
ZEUS  does  not  effectively  predict  the  viscous  layer. 


6.4  CONCLUDING  REMARKS 

The  ZEUS  code  has  been  extended  to  high  Mach  numbers  by  incorporating  the  high-temperature 
effects  of  air  in  chemical  equilibrium.  These  high-temperature  effects  are  accounted  for  by  using  the 
curve  fits  of  Reference  10  to  compute  the  thermodymamic  properties  of  equilibrium  air. 


Various  hypersonic  missile-type  bodies  were  computed  including  sharp  cones,  a  spherically  blunt 
nosed  cylinder  and  a  hypersonic  aircraft.  ZEUS  predicted  the  smaller  shock  layer  thickness  and 
lower  temperatures  characteristic  of  equilibrium  chemically  reacting  air.  The  computed  ZEUS  results 
compare  well  with  measured  high-temperature  effects  and  other  computational  efforts  except  in  regions 
where  viscous  effects  are  dominant.  This  enhanced  version  of  ZEUS  is  an  economic  predictive  tool 
for  missiles  and  will  help  provide  accurate  predictions  in  support  of  future  hypersonic  weapons. 
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HGURE  1.  CARTESIAN  AND  CYLINDRICAL  COORDINATE  SYSTEM  FOR 
ZEUS,  A  SPACE-MARCHING  EULER  SOLVER 
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(A.)  Oblique  Shock  Wave 
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nOURE  4.  SURFACE  (A.)  PRESSURES  (B.)  DENSITIES  (C.)  TEMPERATURES  ON 
CONES  AT  AN  ALTITUDE  OF  100,000  FT.  A  COMPARISON  BETWEEN 
CALORICALLY  PERFECT  AND  EQUILIBRIUM  AIR  RESULTS  FOR 
ZEUS  AND  HUDGINS COMPUTATIONS 
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RGURE  5.  COMPUfED  PRESSURE  CONTOURS  FOR  A  SHARP  CONE  WITH  A 
25°  HALF-ANGLE  AT  MACH  20,  AN  INCIDENCE  OF  20°,  AND  AN 
ALTITUDE  OF  100,000  FT.  RESULTS  ARE  COMPARED  BETWEEN  (A.) 
EQUILIBRIUM  AIR  AND  (B.)  CALORICALLY  PERFECT 
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HGURE  6.  DENSITY  CONTOURS  FOR  A  HYPERSONIC  SPHERICALLY  BLUNT  NOSED 

CYLINDER  AT  MACH  20  AND  AN  INCIDENCE  OF  0°  FOR  (A.)  CALORICALLY 
PERFECT  AND  (B.)  EQUILIBRIUM  CHEMICALLY  REACTING  AIR.  RESULTS 
ARE  COMPARED  WITH  COMPUTATIONS  OF  ZEUS  AND  PALMER^^ 
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(A.)  CALORICALLY  PERFECT  AIR  (B.)  EQUILIBRIUM  AIR 


FIGURE  7.  TEMPERATURE  CONTOURS  FOR  A  HYPERSONIC  SPHERICALLY  BLUNT  NOSED 
CYLINDER  AT  MACH  20  AND  AN  INCIDENCE  OF  0°  FOR  (A.)  CALORICALLY 
PERFECT  AND  (B.)  EQUILIBRIUM  CHEMICALLY  REACTING  AIR.  RESULTS 
ARE  COMPARED  WITH  COMPUTATIONS  OF  ZEUS  AND  PALMER^^ 
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Generic  hypersonic  airplane  configuration. 


O.OiOl  0.002  0.002  0.004  0.005 


P 


HGURE  9.  COMPARISON  OF  COMPUTED  PRESSURE  PROHLES  FOR  EQUILIBRIUM  AIR 
ON  THE  WINDWARD  SYMMETRY  PLANE  AT  MACH  24.5,  AN  INCIDENCE  OF 
1°  AND  Z  =  108.67  FT  FOR  ZEUS  AND  THE  PNS  SOLUTIONS  OF  KORTE'”^ 
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APPENDIX  A 
ZEUS  INPUT  FILE 


The  following  file  describes  the  required  input  for  ZEUS.  To  execute  ZEUS, 
replace  the  ()  symbols  in  the  value  column  with  the  desired  variable  value 
Add  information  supporting  specific  options  between  the  lines  containing 
the  >>>>>  and  <<<<<  symbols.  Variables  enclosed  by  [)  should  be  placed  by 
themselves  on  a  line.  Do  not  eliminate  any  lines  from  this  file. 


VALUE  VARIABLE  DESCRIPTION 


*****INTEGRATION  CONTROL***** 

{)  ZETAEND  Terminate  calculation  if  ZETA  >  ZETAEND. 

0  KEND  Terminate  calculation  if  STEP  NUMBER  >  KEND. 

()  FCFL  Step  size  safety  factor  (1.  >  FCFL  >  0.)  Typically  use  .9  . 

0  XKI  Interior  point  limiting  constant  (2.  >  XKI  >  0.)  Typically  use  1. 

()  lAPR  0,1  for  approximate  or  complete  Riemann  Problem.  Typically  use  0. 

Must  use  0  for  equilibrium  air. 

()  ISHFLUX  0,1  for  using  complete  Riemann  Problem  or  freestream  properties 

to  compute  the  flux  at  the  cell  edge  next  to  the  bow  shock. 
Typically  use  0. 

{)  ISHANG  0,1  for  using  complete  Riemann  Problem  or  Roe  avc ;  jed  variables 

to  compute  shock  angle.  Typically  use  0.  For  equil.LDrium  air,  a 
value  of  1  is  preferred. 


••“•BOUNDARY  CONDITIONS***** 

{)  NXKE  Number  of  edges  at  which  default  lim.iter  setting  won't  be  used. 

Default  settings  are  2  on  edges  1  or  3  and  C  on  edges  2  or  4 . 
For  each  edge  at  which  the  default  limiter  is  non-standard,  add 
a  line  containing  the  following  information; 

zone  number  (KEZ),  edge  number  (KEE)  and  iim.iter  value  (XKE) 
>>>>>start  list  of  [KEZ, KEE, XKE 1 
<<<<<end  list 

{)  NSUR  Number  of  edges  not  featuring  default  surface  types.  Default 

values  are  1  on  edges  1  and  3  and  0  on  edge  2  and  4 .  Program 
will  automatically  account  a  fitted  shock  on  edge  3  or  fin 
surfaces  on  edges  2  or  4 .  For  each  edge  at  which  the  surface 
type  is  not  of  default,  add  a  line  that  specifies  the  zone 
number (KSURZ)  ,  edge  number (KSURE)  and  surface  type  (KSUR)  . 
>>>>>start  list  of  [KSURZ, KSURE, KSUB[ 

<<<<<end  list 

0  DFAC  In  case  of  a  local  wall  Mach  number  which  is  to  small  to  allow 

the  flow  to  be  turned  parallel  to  the  wall,  the  turn  angle  is 
multiplied  by  the  constant  DFAC  (  1.  >  DFAC  >  0.). 


*****0UTPUT  CONTROL***** 
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- Printer - 

()  IPRINT  Print  crossflow  plane  if  step  number  is  evenly  divided  by  IPRINT 

()  NSKIP  Print  n  planes  which  are  evenly  divisible  by  NSKIP 

()  MSKIP  Print  m  planes  which  are  evenly  divisible  by  MSKIP 

()  ISKIP  Print  step  size  if  step  number  is  evenly  divided  by  ISKIP 

- PLOTZA/Force-Pressure  Summary - 

()  JSPPR  In  default  mode,  only  edge  1  of  each  zone  will  be  written  to 

PLOTZA.  JSPPR  is  the  number  of  additional  edges  which  are  not 
to  be  written  in  the  default  manner.  For  each  edge  added  or 
deleted  from  this  print  list, include  a  line  which  specifies 
zone  number  (JSPZ) ,  edge  number  (JSPE)  and  print  code  (JSP)  . 

JSP  =  1  will  write  this  zone  edge  to  PLOTZA,  0  will  not. 

>>»>start  list  of  (JSPZ ,  JSPE,  JSP  ’ 

<<<<<end  list 

0  JSPDIS  1  -  print  grid  point  locations  for  the  edges  specified  by  TSPPR. 

0  -  do  not  O’-  ■  nt . 

0  KSKIP  Write  surface  properties  on  selected  edges  to  PLOTZA  if 

if  STEF  NoMB'r-R  is  evenly  divisible  by  KSKIP 

()  NPRT  Or.  summary  sheet,  print  surface  pressures  of  edges  2  or  4  if 

n  <  NPRT. 

()  DELZA  Write  to  PLOTZA  if  (ZETA  -  ZETA  at  last  write)  >  DELZA 

- PLOTZC  (Plot  file) 

0  IPLOT  Write  to  PLOTZC  if  STEP  NUMBER  is  evenly  divisible  by  IPLOT. 

0  DELZC  Write  to  PLOTZC  if  (ZETA  -  ZETA  at  last  write)  >  DELZC 

0  IPLOTN  Number  of  target  Z  stations  (ZTARGET)  at  w.hich  PLOTZC  will 

be  written.  If  IPLOTN  >  0,  Include  a  list  of  these 
stations  on  the  next  line  (maximum  of  2C).  Stations  must  be 
listed  in  ascending  order. 

>>>>>start  ZTARGET  list 
<<<<<end  list 


•*** "AERODYNAMIC  DATA***** 

0  AREF  Reference  area  used  in  calculating  force  and  moment  coefficients. 

0  XLREF  Reference  length  used  in  calculating  moment  coefficients  and 

center-of-pressure . 

0  lEFORCE  0,  1  don't,  do  print  force  and  moments  for  individual  edges. 


.....separation  MODELING***** 

()  IVIS  0,1,2  for  no  modeling,  clipping,  and  forced  separation.  Clipping 

and  forced  separation  can  only  be  used  in  conjunction  with 
cylindrical  coordinates.  If  IVIS  =  2,  specify  the  number  of 
separation  lines  (NSEP) . 

>>>>>  NSEP 
<<<<< 

For  each  separation  line  specify: 

ISSIDE  -  separation  line  located  between  C  and  18C  degrees, 

1  separation  line  located  between  180  and  360  aegrees . 

I  SEP  -  number  of  poi.nts  used  to  define  the  separation  line. 

ZSSEP  -  Zeta  value  at  which  separation  is  started 
ZESEP  -  Zeta  value  at  which  separation  is  terminated 
PKICD  -  separation  angle  phiC  (degrees) 

PHIAD  -  separation  angle  phlA  (degrees) 

BETACD  -  sparaticn  andgle  BetaC  (degrees) 

3ETAAD  -  separation  angle  BetaA  (aegrees) 

[ISSiDr,  -  Sr?,ZSSr..-,Zusr.  r,.'-.n.C..,,n.n-Au,3E.ACi>,£it.  AA..  . 

<<<<< 
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A  list  of  ISEP  pairs  of  points  (ZEPZ,ZEPP)  describing  the  zeta 
and  phi  coordinates  of  the  separation  line.  This  list  should 
start  at  the  most  forward  section  and  move  aft.  Start  the  list 
for  each  separation  line  on  a  separate  line  and  list  in  same 
order  as  used  in  the  above  data. 

»»>  start  list  of  ZEPZ.ZEPP 
<<<<< 


*****GAS  TYPE***** 

()  IGAS  0-perfect  gas;  !-equilibrium.  air 


() 


0 


0 

0 

0 


{) 

0 

0 

0 

0 

0 


0 

0 

0 


*****UNITS  ***** 

lUNITS  0-english  (p=lbf /f*.  *  *2,  rho=slug/f  t*  *3,  t=rankine,  xyz=ft); 

1-metric  (p=n/m**2,  rho=kg/m**3,  t=kelvin,  xyz=m) . 

All  input  units  must  be  consistant.  Output  units  are  the  same. 


*****INITIAL  DATA***** 

IMOD  0,1  for  new  start  or  restart  respectively.  For  a  restart,  the 
initial  flow  field  is  read  from  TAPE3 .  Additional  information 
is  not  needed  and  the  remaining  entries  in  this  section  should 
be  disregarded.  A  new  start  assumes  a  uniform  flow  field  which 
is  described  in  the  following  section. 

- Mesh  Size - 

IZN  Number  of  zones. 

NA  Number  of  cells  in  xi  direction  (i.e.,  between  edges  1  and  3) 

MA  Number  of  cells  in  eta  direction  (i.e.,  between  edge  4  of  zone  1 

and  edge2  of  zone  IZN) 

List  the  number  of  M  planes  in  each  zone,  starting  from,  zone  1 
and  ending  at  zone  MA. 

>>>>>Start  list  of  MA  for  each  zone 
<<<<<end  list 


- Flow  Properties - 

ZETA  Zeta  coordinate  of  initial  data  plane 
ALPHA  Angle  of  attack  (degrees) 

BETA  Angle  of  yaw  (degrees) 

XMINF  Mach  number 
PINF  Pressure 

DINF  Density 

- Problem  Definition - 

PROBLEM  DEFINITION 

ICOHD  0,1,2  for  cartesian,  cylindrical  or  elliptic  coordinates 

lASYM  0  -  no  pitch  plane  sym.metry;  1  -  pitch  plane  symmetry. 

ISHOCK  Set  to  1  if  edge  3  is  to  be  fitted  by  the  calculation 
(either  as  a  shock  or  sonic  line);  otherwise  set  to  0. 

ISHOCK  =  1  is  only  valid  for  ICORD  =1  or  2.  If  edge  3 
is  to  be  fitted,  the  r  and  phi  (cylindrical  coordinates) 
or  s,  tau  (elliptic  coordinates)  of  the  shock  location  on  each 
M  plane,  (SS,  SSZ)  must  be  defined,  starting  from,  M  =  1, 
etc,  . 

>>>>>Start  list  of  [SS,SSZ; 

>>;>.;>>End  list 


*****MESH  DEFINITION***** 

0  IMESHF  C  for  uniform  mesh  i.n  xi  direction.  I  for  clustered  mesh. 
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a  clustered  mesh  is  to  be  used,  the  clustering  function  (FN)  in 
the  xi  direction  at  the  NA  cell  centers  must  be  defined 
starting  at  cell  1,  located  adjacent  to  edge  1  and  moving  to 
edge  3  (usually  body  to  shock). 

>>>>>start  list  of  FN 
>>>>>end  list 

0  IMESHG  0  for  uniform  mesh  in  eta  direction.  1  for  clustered  mesh.  If 

a  clustered  mesh  is  to  be  used,  the  clustering  function  (GN)  in 
the  eta  direction  at  the  MA(i>  cell  centers  in  each  zone  must 
be  defined,  starting  with  zone  1  and  ending  with  zone  IZN.  Start 
with  cell  1,  located  adjacent  to  edge  4, and  end  at  the  cell 
adjacent  to  edge  2.  The  list  for  each  zone  should  be  started 
on  a  new  card 

»>>»>>start  list  GN 

end  list<<<<<<< 
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APPENDIX  B 

SAMPLE  CASE  OUTPUT 


B.1  HYPERSONIC  SHARP  CONE 

B.1.1  ZEUS  Output  —  n 


1 


PROGRAM  2EUS 


••••  INTEGRATION  CONTROL  • 
1 

0.10000 
so  .  00000 
0 

5000 

C. 90000 
l.OOOOOC 
0 


I  CORD 
ZETA 

ZETAEND  - 
KSTART  - 
KEND 

rcn 

XRI 

lAPR 


ISHFLUX  - 
ISHANG  • 


(0-CARTESIAN.  1 -CYLINDRICAL  COORDINATES) 
(INITIAL  ZETA) 

(MAXIMUM  ZETA) 

(INITIAL  STEP  HUMBER) 

(MAXIMUM  STEP  NUMBER) 

(CFL  SAFETY  FACTOR) 

(INTERIOR  POINT  LIMITER) 

(O-APPROXIMATE  RIEMANN  PROBLEM 
1 -COMPLETE  RIEMANN  PROBLEM) 

(FLUX  AT  SHOCK-  0-RlEMANN  PROBLEM 

l-FREESTREAM  PROPERTIES) 
(SHOCK  ANCLE  -  0-RIEMANN  PROBLEM 

1-ROE  AVERAGED  VARIABLES) 


*•••  BOUNDARY 


NXKE  -  :) 

NSUR  -  0 

DFAC 

ISHOCK  -  1 

lASYM  -  1 


CONDITIONS  •••• 

(NUMBER  OF  EDGES  WITH  NON-DEFAULT  LIMITER  SETTINGS) 
(NUMBER  OF  EDGES  WITH  NON-DEFAULT  SURFACE  TYPE) 
0.9000  (REDUCE  TURNS  NHICH  ARE  TOO  LARGE  BY  DFAC) 

(1 -SHOCK.  O-NO  SHOCK  ON  EDGE  3) 

(O-NO  SYMMETRY. :-P:TCH  PLANE  SYMMETRY) 


SURFACE  POINT  LIMITER 

ZONE  EDGE  1  EDGE  2 

1  2.000  O.OOC 


EDGE  3  EDGE  4 
O.OOC  0.000 


SURFACE  TYPE 

ZONE  EDGE  1  EDGE  2 

1  1  0 


EDGE 

0 


••••  OUTPUT  CONTROL  •••* 


PRINTER 
IPRINT  - 
ISKJP 
NSKIP 
MSKIP 


(flie  6) 

25C  (PRINT 

:C  (PRINT 

3  (PRINT 

€  (PRINT 


FLOW  FIELD  IF  STEP  NUMBER  IS  EVENLY  DIVISIBLE  BY  IPRINT) 
STEP  SIZE  IF  STEP  NUMBER  IS  EVENLY  DIVISIBLE  BY  ISKIP) 

N  PLANES  WHICH  ARE  EVENLY  DIVISIBLE  BY  NSKIP) 

M  planes  which  are  evenly  DIVISIP-E  BY  MSKIP I 
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PRESSURE  SUMMARY  (PLOTZA) 

ESRIP  •  SO  (PRIN?  SURFACE  PRESSURE  IF  STEF  NUMBER  IS  EVENLY  DIVISIBLE  BY  KSRIP) 

WPRT  •  100  lOON'T  PRINT  POINTS  IF  N>NPRT1 

DELZA  -  10000.000  (PRINT  SURFACE  PRESSURE  EVERY  2CTA  -  CELZAI 

JSPPR  •  0  (NUMBER  OF  NON- DEFAULT  EDGES) 

JSPDIS  -  0  (PRINT  SURFACE  GRID  POINT  LOCATIONS) 

SURFACE  PRESSURE  OUTPUT  CONTROL:  O-DON'T  PRINT. 1-PRINT. 


ZONE 

EDGE  1 

EDGE  2 

EDGE  3 

EDGE  4 

1 

1 

0 

0 

0 

PLOT  FILE 
IPLOT 

(PLOT2C) 

10000 

WRITE  TO 

PLOTZC  IF 

STEP  NUMBER  IS  EVENLY  DIVISIBLE  BY  IPLOT 

DEL2C  -  lOOOO.OCO  (WRITE  PLOT  TAPE  EVERY  lETA  -  DELZO. 

IPLOTN  -  0  NUMBER  OF  TARGET  PLOT  STATIONS 


•••••AERODYNAMICS  OATA^»*»* 

AREF  •  1.0000  (REFERENCE  A*>EAi 

ALREF  -  l.QQOQ  (REFERENCE  LENGTH) 

lEPORCE  •  0  (0  -  OONT;l  -  DO  PRINT  EDGE  LOADS) 


•••••SEPARATION 
IVIS  •  0 


MODELING'*^** 

(C  -  IKVISCID. 


CLIPPING.  2  -  FORCED  SEPARATION) 


•••••GAS  TYPE** 

IGAS  -  1 


(O-PERFECT,  1-EQUILlBRIUM  AIR) 


••••••UNITS****** 

lUNITS  -  :  (O-ENCLISH.  I-METRIC) 


C 


IMOD 

MACH 

ALPHA 

BETA 

PINF 

DINF 

HO 

SINF 

IZN 

NA 

MA 

KA/ZONE ; 
IMESHF  - 

IMESHG  - 


lAL  DATA  PLANE  •••• 

0  (0  SELF  START::  FROM  TAPE  START) 

20.0000 (FREE -STREAM  MACH  NUMBER) 
20.0000(ANGLE  OF  ATTACK) 

0.0000(YAN  ANGLE) 

11I5.SOOC(FREE-STREAM  PRESSURE) 

7I61E-02 (TREE-STREAM  DENSITY) 

1  ,  72«9E-»12  (STAGNATION  ENTHALPY) 

0 . COOO (FREE-STREAM  ENTROPY) 

1  (rJMBER  OF  ZONES) 

(CELLS  BETWEEN  EDGES 


II 

36 

36 

C 


AND 

(TOTAL  CELLS  BETWEEN  EDGES 

(C-UNIFORM  MESK. 1 -CLUSTERED  I 
BETWEEN  EDGES  1  AND  3) 
(O-UNIFORM.  1-CLUSTERED  MESH 
BETWEEN  EDGES  2  AND  4) 


3) 

2  AND 


4) 


1 


0 

0 


STEP  - 
MACH  NO. 


0  2ETA  -  :  .COOCSCCE-C:  DZETA  -  0 .  OOOCOOCE-^  CO 

2.0000COCE-0:  ALPHA-  2 . OCOCCCCE- C 1  BETA-  C . DCCCOO CE • CC 


ZONE  NC. 


CELL  PLANE  M  - 


C  - 

3 .OOOOOOIE-O:  C2  -  -3.1076 

-36E-C:  C7  - 

C.CCOOCOOE* 

o:  vx  - 

9.99 

C4I233E-C' 

4 . 36:939C6E-C2 

0 .COCOOOOOE- 

CC  ISUR  3 

B  - 

4.6630766E-02  BZ  -  4 . 6630 764E- 0 1  ST  - 

C  .  000000CE-» 

CO  VX  - 

C.OOCCOCOOE*C' 

C.COOCCCCCE+CC 

O.OOOOOOOOE* 

CC  ISUR  1 

N 

S  T 

P 

- 

;; 

V 

w 

t 

M 

11.5 

2.S2B620E-01  4.363323E-C2  I 

.  1I5500E*C3 

.  7|6;CCE-02 

2.C79C97e*C3  -9. 

C7-»533E»C 

5.7:77:3e«:3 

.  •nociE*:3 

2 . 027C47E* 

16.5 

2.64IC9IE-C:  4.363323E-02  I 

.  1I55C0E*C3 

.7I6-.CCE-C2 

2  079:97e*C3  -9. 

C77533E*: 

5  .  7:77;3E*:3 

.  II70CIF*C3 

2  ,  C27047E* 

13.5 

2.225I16E-0:  4.363323E-C2  1 

.  ie55CCE*C3 

.  7|6:0CE-02 

2.079C97e*C 

-9. 

C'’7533E-C 

5.7:77:3e*c3 

.  •l7C0tE*C3 

2  .  C27C47E‘ 
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200  0.663  1.362E*05S  l.S45E-»05S  1.7nE«05S  1.I62E4055  2.05:E-»0SS  2.219E-»05S  2.393E405S  2.S71E405S  2.739E-405S  2.691E405S 

300  1.005  1.3I6E4Q5S  1.557E*0SS  1.721E->C5S  l.tfSE^CSS  2.CS2E405S  2.221E405S  2.394E40SS  2.56iE4-05S  2.732E405S  2.I63E4C5S 

400  1.463  1.367E-f0SS  1.55SE40SS  1.719E405S  1.6d5£->C5S  2  0S2E4OSS  2.221E*0SS  2.394E«0SS  2.S6IE4055  2.733E«05S  2.6I3E4CSS 

500  2.117  1.367E-405S  l.S54E«05S  1.719E-405S  l.tlSE-^OSS  2.053E40SS  2.222£«05S  2.395E405S  2.S69E405S  2.733£«05S  2.i64E40SS 

600  3.057  1.366E405S  1.553E-»'05S  1.719E405S  l.ilSE405S  2.053E405S  2.222E405S  2.394E405S  2.569E40SS  2.733E-4C55  2.I63£-»^C5S 

700  4.411  1.3I6E-405S  1.553E405S  1.719E405S  1.I65E40SS  2.053E40SS  2.222E-»05S  2.396E-»05S  2.569E405S  2.733E-»05S  2.6I4E<405S 

600  6.361  1.3I6E-*05S  1.553E40SS  1.719E405S  1.I64E40SS  2.0S3E«05S  2.222E40SS  2.394E*0SS  2.S6IE405S  2.733£«05s  2.6t3E405S 

900  9.170  1.385E405S  1.553E*05S  1.719E405S  1,66SE40SS  2.0S3E*05S  2.222E*05S  2.395E*05S  2.566E4-05S  2.733E405S  2.I64E+05S 

lOOC  13.220  1.385E-0SS  1.553E-»05S  1.7l9E*05S  1.684E4e5S  2.052E+05S  2.222E*C5S  2.395E40SS  2.568E405S  2.733E4C5S  2.884E*CSS 

1100  19.057  1.385E405S  1.5S3E4-05S  1.719E405S  1.886E405S  2.O53E+0SS  2.222E405S  2.395E405S  2.568E405S  2.733E405S  2.i84E405S 

1200  27.471  1.385E405S  1.553E405S  1.719E*05S  1.884E*05S  2.051E<405S  2.222E40SS  2.395E405S  2.569E4055  2.733E«05S  2.885E<4055 

1300  39.600  1.385E40SS  1.S53E405S  1.720£4>05S  1.885E405S  2.053E405S  2.223£*05S  2.394E405S  2.568E405S  2.733E4-05S  2.883E405S 

1364  49.975  1.385E40SS  1.553E+05S  1.719E40SS  1.885E405S  2.0S3E-0SS  2.222E*05S  2.395E*05S  2.568E405S  2.732E405S  2.884E*05S 


1 

ZONE  «  1 

EDGE  •  1 

MACK  NO.  •  20.000 

$  U  R 

ANGLE  OF  ATTACK  •  20.000  YAW  ANCLE 

FACE  PRESSURES 

OSTEP 

ZETA 

31 

32 

33 

34 

35 

36 

1 

0.104 

3.360E4C5S 

3.474E405S 

3. S67E40SS 

3.638E4C5S 

3.687E«05S 

3.711E*CSS 

100 

0.459 

3.365£*0SS 

3.363E405S 

3.357E405S 

3.381E*0SS 

3.434E*0iS 

3.467E405S 

200 

0.683 

3.023£405$ 

3.134E4055 

3.229E405S 

3.306E4C5S 

3.363E405S 

3.388E405S 

300 

1.005 

3.016E40SS 

3. 131E405S 

3.229E405S 

3.304E405S 

3.359E*05S 

3.383E4C5S 

400 

1.463 

3.018E405S 

3.133E405S 

3.229E4C5S 

3.305E*05S 

3.3S9E405S 

3.384E405S 

500 

2.117 

3.018E405S 

3.133E405S 

3.228E405S 

3.305E405S 

3.361E*05S 

3.384E*055 

600 

3.057 

3.017E40SS 

3 . 133E405S 

3.229E405S 

3.30SE40SS 

3.360E«05S 

3.384£-»05S 

700 

<•.411 

3.017E405S 

3.133E405S 

3.229E40SS 

3.306£*05S 

3.361E405S 

3.384E405S 

800 

6.361 

3.018E405S 

3.133E40SS 

3.229E40SS 

3.305E405S 

3.361E*05S 

3.384E405S 

900 

9.170 

3.017E*05S 

3.133E*C5S 

3.229E4CSS 

3.30SE*05S 

3.360E*05S 

3.385E*C5S 

1000 

13.220 

3.017E405S 

3.133E405S 

3.229E4C5S 

3.30SE405S 

3.358E*05S 

3.384E'»05S 

1100 

19.057 

3.017E-05S 

3.133E-»05S 

3.229E4C5S 

3 , 306E4C5S 

3.360E403S 

3.385E-C5S 

1200 

27.471 

3.018E405S 

3 . 133E405S 

3.229E*05S 

3.30SE-»05S 

3.360E*05S 

3.385E405S 

1300 

39.600 

3.017E405S 

3.132E405S 

3.229E405S 

3.305E4C5S 

3.360E-05S 

3.384E405S 

1364 

1 

0 

49.975 

3.017E405S 

3. 133E4C5S 

A 

FREE 

3.228E405S 

E  R  0  D  Y  N 

S  T  R  E  A 

3.30SE-05S 

A  M  I  C  D 

M  C  0  K  D 

3.360E+05S 

A  T  A 

:  T  I  O  N  S 

3.384E405S 

MACH 

NO.  -  20. 

000  ANGLE  OF 

ATTACK  -  20 

0000  ANGLE 

OF  SIDE  SLIP 

-  0.0000 

REF 

AREA  -  1.000000E400  REF  LENGTH  -  I.OOOOOOE4OO 

OSTEP 

ZETA 

cx 

CY 

cz 

CKX 

CMY 

CMZ 

XCP 

YCP 

1 

0.107 

5.88457E-04 

O.OOOOOE40C 

4  .70158E-04 

O.OOCOOE-OC 

7.41649E-C5 

0 .00000E4CC 

1 .26C33E-C: 

0.O0C00E*CC 

100 

0.460 

7.01070E-02 

0.O0O00E40C 

S.50821E-02 

O.OOO0OE4O0 

2.75I95E-02 

O.OOOOOE-00 

3.92535E-01 

O.OOOOOE-OC 

200 

0.685 

1.63137E-01 

O.OOOOOE40C 

1 .26780E-01 

O.OOOOOE*OC 

9.30832E-02 

0 . 00C0CE*0C 

5.70584E-C: 

0.  OOOOE-CO 

300 

1.007 

3.59929E-01 

O.OOOOOE400 

2.79272E-0: 

O.OOOOOE-»OC 

2.98114E-0: 

C.OOCOCE  OC 

•.28259E-01 

O.OCOOCE-OC 

40C 

1.465 

7.68080E-C1 

O.OOOOOE40C 

5.96013E-C1 

O.OOOOCE-OC 

9. 192:8E*0: 

C.00CCCE*CC 

1 . 19677E*0C 

C.OOOOOE-CC 

50C 

2 .121 

1.61481E400 

0.00000E«0C 

1  .2534CE*0C 

C.00CC0E-»0C 

2.78758E*C0 

C.00000E-»0C 

1  .  72626E*CC 

C  .  OOCCOE-OC 

600 

3.063 

3.37413E400 

0.00000E40C 

2.6:921E40C 

C-0C000E*0C 

8 . 3984;E'»0C 

0.00000E*CC 

2.48906E«0C 

O.OOOOOE^OC 

700 

4.419 

7.02909E40C 

O.OOOOOE400 

5.45666E«CC 

C  .  0CCC0E*0C 

2 . 52216E*C: 

D.00C00E*CC 

3.588;7E*00 

C  .  OOOOOE-OC 

800 

6.372 

1.46236E401 

O.OOOOOE4OC 

:  .  1350IE«C: 

C  .  OCOCCE^OO 

7.56375E«C: 

0. OOOOOE^CO 

5  . 17228E*CC 

O.OOCCOE-CC 

900 

9.187 

3.04037E401 

O.OOOOOE-OC 

2.3598:E*C1 

O.OCOOCE-OC 

2.26676E*C2 

O.OOCCOE'CC 

7.45556E40C 

O.OOOOCE-CC 

1000 

13.244 

6.31926E401 

O.OOOOOE-*00 

4 . 90456E«C: 

O.OOOCCE-CC 

6  79121E*02 

0  .  OOCOOE^CC 

1  .  07468E*0: 

O.OOOOOE-CC 

1100 

19.092 

1.31327E-02 

0.00000E40C 

1  .r'.922f.»02 

0 . 0CC0CE*0C 

2.0344:E*C3 

O.CCC0CE*CC 

1  .S49:2E«01 

O.OOOOCE-CC 

1200 

27.521 

2.72901E402 

C.OOOOOE-CC 

2 . : :79:i*c2 

O.CCCCCE»CC 

6.C9393E403 

0.0000CE»0C 

2.23302E40: 

0  .  OOCCCE-CO 

1300 

39.672 

5 . 67081E402 

0 . oooooE*c: 

4  .  4CC8IE*C2 

C.C0C00E*CS 

1.82b35E‘04 

O.OOOOOE-OC 

3.21885E*01 

C.OOCOCE-CC 

1364 

50 . OCO 

9. 00777E402 

0 .  OOOOCE-*OC 

6.99054E«C2 

C  .CD00CE*CC 

3.65426E*04 

O.OCOOGE*0.: 

4  .OS679E*01 

O.OOOOCE-OC 

NORMAL  STOP 
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B.2  HYPERSONIC  BLUNT  NOSE  CYLINDER 


B.2.1  BLUNT2  Output 


IHypcrtonlc  Blunt  Nose  Cylinder,  H*2C.  Alph«>0.  Alt«20.0CC  m 

EXTERNAL  PARAMETERS-KENSE-  OOlOOOOOOO 


TEST(  1) 
TEST!  2) 
TEST<  3) 
TE$T<  4) 
TEST(  5) 
TEST (11) 
TEST(12) 


l.OOOOE-04 

4.0000E-03 

5.0000E-07 

S,0000E-04 

5.0000E-04 

1.0000E-0< 

1 .OOOOE-04 


IRG  -  1  AIR 

RIDEAL  EIDEAL  GIDEAL  RREAL  PO  HO  CO  PREF  TREF  SRREF 

2.t70SE+02  l.OOOOE  +  00  1.3990E+00  2.32S5E*03  l.OOOOE+00  2.0245E^01  2.3222E-»02  3.1450E+01 

♦  2.t70SE+02  1.0133E+05 

DEL  -  1  .OOOOOOOE-fOO 

IMAX  -2000  IMIN  -2000  IHOO  -  2i  IPRINT-2000 

IRSPCE-  0  ISTART-  >1 

JMAX  -  50  JPRINT-2000 

KMAX  -  25  HKAXI  -  2S 

NMAX  -  II 

RN  -  S.OOOOOOOE-01 

STAB  -  1.0000000E-»00 

XO  -  5.0000000E>01  XSTAG  -  0 .  OOOOOOOE-^OO 

YMAX  -  1.2000000E+C2  YSTAR*  4 . 5000000£-»01 


EEUS  RELATED  INPUT: 

NA  •  36  MA  •  34  lASYK  •  1 

ISTART  -  5.0000000E-01  20  -  5 . OOOOOOOE-01  ALPHA  •  2  .OOOOOoOE+01  BETA  •  0 . 0000000E4  00 

IBL  •  0  NSTREM  •  9  FTCOKV  -  1 . OOOOOOOE^OO  lUNITS  -  1 


FREE  STREAK  CONDITIONS 

ALTITUDE  -  20000. 

1959  ARDC  TABLES 

?  S/R  H  0  M  RHO  T  2  G  OQ 

S.S29IIE403  2.53209E401  2.14947E405  5.I73I0E*03  I.OOOOOE^Ol  I.IBUlE-OI  2.16662E402  l.OOOOOE^OO  1.3I4I7E400  1 . 533I5E404 


91  ENTRIES  IN  SHOCK  TABLE 


UWl/OKMIN)  -  l.OOOOOOE-Oi 


UMl/OKMAX)  -  1 . 100000E400 


PSTAG-  2.94I5E+04  TSTAG-  7.3519E*03  RHOSTG-  l.0549E*00  GSTAG-  1.1845E*00  HTOT-  1.744IE407 

PSH-  2.I032E  +  04  TSH-  7.3:i9E'»03  RHOSH-  I.OIOIE  +  OC  CSH-  1.1I39E*00  HSH-  1.7334E*07 

5RSH-  3.9222E401  AKSH-  2.IS24E-0I 

DSH  -  7.0793211E-01 

IHypertonlc  Blunt  Nose  Cylinder.  K-20.  Alpha-0.  Alt-20,000  re 
DEFAULT  SPHERICAL  GEOMETRY  USED 
RN  -  5.00000E-01 

IHypersonlc  Blunt  Nose  Cylinder.  M-20.  Alpha-0,  Alt-20,000  re 


K 

Y 

RB 

RBTH 

1 

O.OOOOOOOE-^OO 

5 . OOOOOOOE-01 

O.OOOOOOOE*OC 

2 

8.724444SC-02 

5 . OOOOOOCE-O: 

0 .OOOOOOOZ*OQ 

3 

1 . 7453294E-01 

5.0000000E-01 

O.OOCOOOOE*OC 

4 

2.4’79940E-01 

5.0000000E-01 

o.oooooooE«oo 

5 

3.4904587E.C: 

5.0000000E-01 

O.OOOOOOOE40C 

6 

4 .3633235E-01 

5.000C0COF>Ol 

-1 .53554l5E-0e 

7 

5.2359I79E-01 

5.0000000E-01 

O.OCCCOOOE+OO 

1 

6. 10I6529E-01 

5.0000000E-01 

0.0000000E40G 

9 

e.9«:3174E-01 

5.0000000E-C1 

0.0000000E«0C 

10 

7.i539i:9E-C: 

5 . OOOOCOOE-O: 

2.56927l9E-0t 

11 

I.7246469E-C: 

5,OOOCCOOE-0: 

0 . 0000000E400 

12 

9.5993114E-01 

5 . OOOOOOOE-01 

C.OOOCOOCE*CO 

13 

1 .0471976E400 

5.0000000E-0: 

o.oooooooE*oo 

14 

1 . 1344440£«00 

5.000000CE-Ci 

0 . OOOOQOOE*CO 

15 

1 .2217306E-»00 

5 . OOOOOOOE-OI 

O.OOOOOOOE*OC 

14 

1 .30I9970E40C 

4.9999997E-01 

0.00G0000E*CC 

17 

1 .394263SE400 

5.0000000E-01 

0.0000000E«0C 

B-6 
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le  1 .4I3S2»9E>00  5.0000000E-01  0 . OOOOOOOB^OO 

19  1 .57Q1964E*00 '  S . OOOOOOOE-01  0 . OOOOOOOE+OO 

20  1 .4SIO€29£-*^00  5  .  OOOOOOOE-01  0  .  OOOOOOOE-»00 

21  1 .74S3294E<»00  S .  OOOOQOOE-01  O.OOOOOOOE-^00 

22  l.i32S9SIE+00  S . OOOOOOQE-01  0 .OOOOOOOE+OO 

23  1.919l423E-f00  3 .  OOOOOOOE-Ol  0 . 0000000E-»00 

24  2.0071297E400  S . OOOOOOOC-01  Q .OOOOOOOt^QO 

25  2.0943952E400  5 . QOOOQOOE>01  0 . OOOOOOOE«00 


IRypereonlc  Blunt 

.  Note  Cylinder, 

.  H-20,  Alph*-0, 

.  Alt-20. 000  m 

STEP 

TIME 

DELTA 

PSTAG 

WRMS 

1 

5.5728725E-07 

6.5766737E-02 

2.9429010e406 

5.0127559E+00 

2 

1.1144186£'06 

6.5759756E-02 

2.9360023E-f06 

9.4505558E400 

3 

1.6714025E-06 

6.5748170E-02 

2.9278905E406 

1 .4038545E401 

4 

2.2283821E-06 

6.S731928E-02 

2.9186713E>06 

1.8675589E«0l 

5 

2.78S3721B>06 

6.S711141E-02 

2.9084403B-f06 

2.3286003e40l 

(mi/OKHIX)  -  4.7524f0E-01 

imi/QKHM)  - 

1.003999B4^00 

IHypereonlc  Blunt 

.  Nose  Cylinder, 

,  M-20,  AIphe-0, 

,  Alt-2Q. OCC  n 

STEP 

TIME 

DELTA 

PSTAG 

«RMS 

201 

1 .0931705E-04 

4.0710162E-02 

2.9427413E406 

2.6345813E402 

202 

1 .0978320E-04 

4.0722951E-02 

2.9403565E406 

2.6139301C402 

203 

1.1024I78E-04 

4.0736869B*02 

2.9376098E-^06 

2.5923465E402 

204 

1.1071380E-04 

4.0751837E-02 

2.9345803E-«>06 

2  569925SE402 

205 

1.1117827E>04 

4.07677I9E-02 

2.9313635£^06 

2.S467589E402 

KMAX  - 

25  NMAX  -  It 

-  3.239257£>-01 

UNl/OKMAX)  - 

9.966639E-01 

STEP 

TIME 

DELTA 

PSTAG 

ITRMS 

401 

1.9852728£‘04 

3.7458688E-02 

2.9282870e<»06 

4.5209958E-402 

402 

1.9895344E-04 

3.7450656E'02 

2.9283540E406 

4.4930344E402 

403 

1.9937946E-04 

3.74428:iE>C2 

2.9284715E-»06 

4.4651422E4C2 

404 

1 .9980536E-04 

3.7435144E-02 

2.9286308E-»06 

4.4373187E402 

405 

2.0023114E-04 

3.7427694E-02 

2.9268176Ef06 

4.4095642E40P 

MMAX  - 

25  NMAX  -  18 

UMI/QKMIN)  -  3.4347I4E-01 

UN1/Q1(HAX1  - 

9. 983664E-01 

STEP 

TIME 

DELTA 

PSTAG 

MRMS 

601 

2.8191469E>04 

3.6079410E-02 

2.9361170E+06 

9.7029724E401 

602 

2.8232345E-04 

3.6075525E-02 

2.9360975E406 

9.6875244E401 

603 

2.8273216E-04 

3.6071673E-02 

2.9360833E406 

9.6731621E401 

604 

2.8314081E-04 

3.6067856£'02 

2.9360750E4Q6 

9.659e328£40l 

MMAX  - 

25  KMAX  -  16 

UNl/QKMZN)  -  3.616154E>01 

IH*1/Q1(MAX)  - 

9.991269e'01 

STEP 

TIME 

DELTA 

PSTAG 

MRMS 

801 

3.6272925E-04 

3.5485916E-02 

2  SSOSTSSE^OS 

8.2728256E*0l 

802 

3. 6312939E-04 

3.5483882E-02 

2. 9385738E*C6 

6.27048e7E40: 

•  03 

3. 635294 9£<04 

3. 5481848E-02 

2.9385738E406 

8.2679306C*01 

804 

3.6392955£>04 

3.S479810E-02 

2.93e5738E406 

8.26514S1E401 

•  05 

3.6432958E<04 

3.5477776E-02 

2.9385738E406 

8.262n93£40l 

KKAX  -  25  KKAJC  «  19 

Wl/QKMIN)  -  3.57SI45E-01  UK1/Q1(MAX)  -  9.995127E-01 

STEP  TIME  DELTA  PSTAG  NRMS 

1001  4.422S4I9E-04  3 . S176776E-02  2 . 9396670E*06  4 . 5741 1 tOE^Ol 

1002  4 .426S250E-04  3  .  S17S651E-02  2 . 9396725E406  4 .5658ee2£-»01 

1003  4.4304I09£>04  3. 5174560E-02  2. 9396B38£«04  4 . 5577564E401 

1004  4 .4344366E-04  3 . 5173472E-02  2 . 9396923E406  4 . 5497143E401 

1005  4.4343921E-04  3 . 51 72361E-02  2 . 939697IE406  4 .5417595E401 

MKAX  •  25  NKAX  *  16 

UWl/QKMIN)  -  3.5S2336E-C1  UWl/QKMAXl  -  9.997421E-01 

STEP  TIME  DELTA  PSTAG  MRMS 

1201  5.2111130E-04  3 . 5010364E-02  2 . 940S473E*06  2 . 74 54748E401 

1202  5.2150449E-04  3 . 5009763E-02  2 . 9405500E406  2 . 7373650E4 01 

1203  5.2169769E-04  3 . 5009202E'02  2 . 9405500E406  2 .7293016C401 

1204  5.22290I6E-04  3 . 5006620E- 02  2 . 9405526E«06  2 . 721 2e39E«01 

1205  5.22684C7E-04  3 . S008039E-02  2 . 9405585E*06  2 . 7133 1 1 0E*01 

1246  5.3879584E-04  3 .4995255E~02  2 . 9406680E*06  2  .421 1 7t4£fCl 

1247  5.3918662E-04  3 . 4 98474 9E-02  2 . 9406680E406  2 . 4 147606E«01 

1246  5.3958141E-04  3 . 4 984238E- 02  2 . 9406708E«06  2 . 4083683£*0 1 

1249  5.3997413E-04  3 . 4963728E-02  2 . 94  C  67  3SE>  0  6  2 .4020014E-»0; 

1250  5.4036686E-04  3 . 4963221E-02  2 . 9406735E406  2 . 3956598E«01 

MMAX  -  25  NMAX  -  18 

UMl/OKMIN)  -  3.547I83E-01  UWl/QKMAXI  -  9.997761E-01 

COKVERGENCE  CRITERIA  SATISFIED  AT  STEP  1258 
lHyp«r»orUc  Blunt  No»«  Cylinder.  H*20.  Alpha^O.  Alt-20.000  m 

STEP  TIME  DELTA  '•STAG  MRMS 

1251  5.4075959E-04  3 . 498271 1E-C2  2 . 9406763E+06  2 . 3893423E*C 1 

1252  5.4115232E-04  3 . 4 982204E-C2  2 . 9406790E«06  2.3630505E40I 

1253  S.4154504E-04  3 . 4981694E-02  2 . 9406790E406  2 . 3767822E40 1 

1254  5.4193777E-04  3 . 498 1 1 I7E-02  2 . 9406790E*06  2 . 37053^3E*01 

1255  5 .4233044E-04  3 . 4 98067 7E-02  2 . 9406820E406  2 .3643064E«01 

1256  5.4272311E-04  3 . 4  98  0 1  70E- 02  2 . 9406848E-»06  2 . 3581018E*01 

1257  5.4311578E-04  3 . 4979660E-02  2 . 9406875E*O6  2 . 351 91 59E*01 

1258  5.4350845E-04  3 . 49791 53E-02  2 . 9406933E*06  2 .3457474E«D1 

MMAX  -  25  KMAX  -  18 

UWl/OKMIN)  -  3.S47164E-01  UM1/01(MAX)  -  9.997822E-01 

IMypersorslc  Blurt  Ho*e  Cylinder.  M-2C,  Alph4-C,  AJC-20.000  tr. 

STEP  **•  DT  -  3.9266E-07  AT  W  •  1 8  M  -  7 

TIME  -  5.4351E-04 

I  -1258 

M  y  RST  RSTH  PRESSURE  CHECK  XB  YB  XS  YS  SIGMA 
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2.5941E403  *3.5538E*03 

3  7367E401 

3.5057E40C 

7.9189E-01 

-2.5703E-03 

I  - 

1258 

N  -  12 

H 

P 

RHO 

T 

U  V 

S/R 

MACH 

R 

DK 

1 

2.893SE*06 

1 .0385E4-00 

7.3357E+03 

-2 .987SE402  O.OOO0E4OO 

3.9218E401 

1.6445E-01 

5.2263E-01 

-7.7492E-04 

3 

2.8043E-*06 

1 .0188E-»00 

7.2788E+03 

-2.8960E402  •7.S843E402 

3.9120E401 

4.4980E-01 

L.2316E-01 

-7.6270E-04 

5 

2.5505E<»^06 

9.6040E-01 

7.1164E+03 

-2.S948E402  -1.4I86E403 

3.884ae401 

8.5303£>01 

5.24S2E-01 

-7.7415E-04 

7 

2.I790E4-06 

8.7156E-01 

6.B265E>03 

-1.9974E432  >2.1«02E403 

3.8427e401 

1.2(33E400 

5.2795E>01 

I.1182E  05 

9 

1.7S0€E+06 

7.5803E-01 

6.4521E403 

-9.4250E401  -2.7449E403 

3.7929E401 

1.8658E400 

S.3312E>01 

-3.7701E-03 

11 

1.3337E*0« 

6.3871E-01 

5.98572403 

7.4852E402  -3.2300E403 

3.7415E401 

2.0575E400 

5.4120E-01 

-9.2729E-03 

13 

9.7204E-^05 

5.1322E-01 

5.51C2E4Q3 

3.1809E402  -3.5871E403 

3.7012E401 

2.4111E400 

5.5338E-01 

-1 .3239E-C2 

15 

6.9337E405 

4.0760E-01 

5.0299E403 

6.3BltE402  -3.8333E403 

3.««€5E401 

2.7421E400 

5.7153E-01 

-1 .6039E-02 

17 

4.9392E-I-05 

3.1355E-01 

4 .7034E'»03 

1.0425E403  "3.9969E4G3 

3.8569E401 

3.0261E400 

5.9816E-01 

-2.5350E-04 

19 

3.5390E+05 

2.42S0E~01 

4.4065E403 

1.4920E4C3  -4.0296E403 

3.e44££401 

3.2688E400 

6.3673E-01 

4.5977E-03 

21 

2.5599E*05 

1.9352E-01 

4.0B08E403 

1.9865E403  -3.9«46t*03 

3  - 6:25E401 

3.5320E400 

6  9248E>01 

-•.8475E-03 

23 

1 .8764E-»05 

1 . 5147E-01 

3.8735E4C3 

2.5371E-»03  -3.8222E4C3 

3.5992E401 

3.7680E400 

7.7472E-01 

-1.9‘l5E-03 

25 

1 .3991E+05 

1 .2227E-01 

3.8449E403 

3.0949E403  -3.5721E403 

3.5668E401 

4.D254E400 

9.013SE-01 

-5.9186E-03 

1  - 

1258 

N  -  15 

M 

P 

RHO 

T 

U  V 

S/R 

MACH 

R 

OH 

1 

2.856SE406 

1.0269E400 

7.32S6E403 

'4 . 0274E402  C.OCCCE^OO 

3.  9218E*01 

2.2131E-01 

5.2881E-01 

-5.7846E>04 

3 

2.7757E+06 

1.0127E+00 

7.2596E4C3 

-3  8924F.C2  -8.4:72L-C2 

3.908eE4Cl 

5.1494E-01 

5.2947E-01 

-5.9706E-04 

5 

2.5465£-»06 

9. 7107E-01 

7.0611E-»Q3 

-3.4,549L*02  -1.6493E403 

3.8704E4C1 

9.5743E-01 

5.3159E-01 

>8.5163E-04 

7 

2.2104E-»06 

9.0622E-01 

6. 7247E*03 

-2.6005E402  -2.3855E403 

3.8134E-»01 

1 .4150E400 

S.3557E-01 

'6.4486E-04 

9 

1 .8209E<>06 

8.2604L-C1 

6.2450E4Q3 

-•1.1452E+02  -3.0212E403 

3.7430E-»01 

1.87S9E400 

5.4215E-0i 

-4.54B6£>03 

11 

1 .4361E-»06 

7.3726E-01 

5.6735E403 

1.0967E402  -3.5387E*C3 

3.6698E401 

2.3361E400 

5 . 5243E-01 

-1.14S2E*C2 

13 

1 .0943£<^06 

6.4275E-C: 

5.0715E403 

4.2374E402  -3.9181E403 

3.SO65E401 

2.  ?769E-»00 

5.6793F  O: 

-8.2090E-03 

15 

8.2020E«05 

5.4821E-Q: 

4.S670E-»03 

8.28e5E402  -4.1806E4C3 

3 . 5513E401 

3.1961E*0C 

5.9104E-01 

-1  .5545E-03 

17 

6.1567E*05 

4 . 5984E-01 

4.1911E-»03 

1 .3090E+03  -4 .3C13E403 

3.5035E-01 

3.5532E400 

6.2493E-01 

-5.4294E-04 

19 

4.6403E-*^05 

3.8S40E-01 

3.8677E403 

1.8491E403  -4.3141E403 

3.4534E401 

3.8968E400 

6.7402E-0: 

-4.0585E-03 

21 

3.49e3E-»'C5 

3.2132E-01 

3 . 5817E403 

2.4218E-»03  -4.2172E*03 

3.4037E401 

4.2256E400 

7.4497E-01 

>€.2888E'03 

23 

2.6439E-»C5 

2 .677SE-C1 

3.3154E-*03 

3.0232E*C3  -4.0086E403 

3.3561E401 

4.5577E400 

8.4964E-01 

-4 .3620E-03 

25 

1 .9887E«05 

2.2716E-01 

2.9998E403 

3.6303E*03  -3.6902E403 

3.2972E40i 

4 . 9568E400 

1 .0108E40C 

-4 .2408E-03 

I  -1258  N  -  18 


K  P  RHO  T  •  U  V  S/R  MACH  R  DH 

1  2.80m*06  1.0100E*00  7 . 3 1 0€E-»03  -  5  . 1 82  “JE*  02  0.OOOOE*OO  3.9219E*01  2.I597E-01  5.3*98E-01  -3.9011E-04 

3  2.7316E*C€  1.0C20E*00  7.2347E*03  -4.9969E-02  -9.2620E»C2  3.9047E-0:  S.SSOIE-Ol  3.3579E-01  -3.8892E-C4 

S  2.5323E-06  9.8003E-01  6.9953E+03  -4.4002E*02  -1.81S5E*03  3.8549E*0l  l.C694E*00  5.3836E-01  -4.3017E-04 

7  2.2407E-^06  9.4S33E-01  6.6073£*03  -3.2508E*02  -2.6224E*03  3.7i04E*0J  1.58nE*00  5.4319E-0:  -7.8297E-04 

9  1.9032E*06  9.0989E-0:  6.0172E*03  -1.3488E*02  -3.3228E*C3  3.E585E-01  2.1*74E*00  5.5ai8E-Ci  -3.9637E-04 

II  :.S871£«06  8.8496E>0I  5.2634E403  1.S160E402  -3.888CE«03  3.S8ME«0l  2.68eSE*0C  S.8367E-C;  '1.IS04E-03 

13  1.2687E-.06  8.6562E-0:  4.5359E  +  03  5.4399£<»02  -4.3059E*03  3.4789E*C1  3.2852E*0C  5.8249E-01  -6.8939E-04 

15  1.0201E*C6  8.1I23E-01  4.C2€3E*03  1.0310E-»03  -4.5S57E*C3  3.3787r»0l  3.80B9E*0C  8.3CS5£-C1  -9.3299E-04 

17  •.2228E+05  7.S790E-01  3.62C3E+03  1.5928£*03  -4.6531E*03  3.2948E*01  4.2720E+0C  6.5170E-01  -1.4128E-03 

19  6.6379E*05  6.9963E-01  3.2391E*03  2.2042E-»C3  -4.4C77E+03  3.2191E*0l  4.7114E*0C  7.1131t-01  -1.9228E-03 

21  5.32S5£*0S  $.4705E-0:  2.I490E*C3  2.8333E+03  -4.4339£»C3  3.1499E+01  5.1739E+00  7.9747E-01  -2.0899E-03 

23  4.233IE*05  6.02I9£-01  2.4439£«03  3.4S14£-»03  <4.I324£-^03  3.0850£«0l  5.688CE*00  9.2458E-0I  -.3.46S9E-D3 

25  3.3099E+05  S.6701E-01  2.0336E-*’03  4.03S2E*03  -3.71S9E*C3  3.01t9E*01  8.3132E*0C  1.12C3E-0C  -5.1491E-C3 

12EUS  INITIAL  DATA  GENERATED  BY  PIT 

I2N  -  1  NA  -  36  MA  •  36 

PINT  -  5.52988E*C3  DINE  -  8.89148E-02  XMINF  -  2C.0C  ALPHA  -20.00  BETA  -  C.OO 

KSTART  -  0  2START  -  5.000COE-01  lASYH  -  1  ICORD  -  1  ISHOCK  -  1 

XFX,  XFY,  XF2  -  7.56896E-»04  O.OOOOOE+00  2.33032E*05  XKX.  XKY,  XM2  -  C.OOOOOE+00  1.24163E'»05  O.OOOOOE-CO 

PLANE  M  -  1  PHI  -  2 . 50C 


N 

R 

P 

RHO 

u 

V 

N 

5.05893E-C: 

4  .  50897E-*04 

3  .  n47lE-C 

1 .27880E*C2 

-6.0€:03E*C: 

4.19273E*C3 

3 

5.29464E-CI 

6. 19923E-04 

4 . 17315E-C 

5.95243E*02 

-3.787e8E401 

4 . 11372E*03 

5 

5.53036E-01 

7.76170E*04 

5.07866t-0 

2  9.17b66E*02 

-2.24753e*01 

4.0055tE*03 

7 

5.76607E-C1 

9 . 20749E«04 

S.89468E-0 

2  :.lth09t*G3 

-1 .02045E'»C: 

3.t9545£*03 

9 

t.00179£-Cl 

1 .05993E*05 

6. 75C26E-0 

2  1.382€3E*03 

-6. 5809CE-01 

3.82352E*03 

11 

6.23750E-0: 

1 . 19354E*05 

7.60630E-0 

2  1.57321E*03 

7 . 79345E*C0 

3-75556E*03 

13 

6.47322E-01 

1.32654E*C5 

8.66395E-C 

2  1.75t42E403 

1 .6:820E*C: 

3.7317eE.C3 

15 

6.70893E-01 

1 .46549E*C5 

9.9474:e-0 

2  1.89733E403 

2.27«78E401 

3.752S3E*C3 

17 

6.94464E'01 

1 .61006E40S 

1  .  11496E-C 

2.037nE403 

2.86028E-»01 

3.73170e*C3 

19 

7. lt036E-C: 

l.76311E*05 

1  .24C23E-C 

2.:785SE403 

3.434:6E*01 

3.69344E*C3 

21 

7.41607E-C: 

1.92797E405 

1  .  39335E-C 

2.3255fE-'C3 

4  .OJ628E-»02 

3. 67$79E*03 

23 

7.65179E-C1 

2 . 10943E*C5 

1  .  5958IE-C 

2.47622E403 

4.60887E*C1 

3.70605E*03 

25 

".88750E-C: 

2.315:8E-*CS 

1  .85743E-C 

2.62524E403 

5. 17797E*01 

3.76957E*C3 

27 

8.12322E-0I 

2 . S5723E«C5 

2.  19139E-C 

2.77733E403 

5.7393;E4C1 

3.8388I£«03 

29 

8 . 35893E-01 

2  .  B4824E-»C5 

2 . 6:990E-C 

2 . 93297E*03 

6.3044CE*C: 

3.89883E«C3 

31 

8.5946SE-01 

3.20557E*05 

3.  19316E-C 

3.0B639E403 

6 . 86052E401 

3 .95446E*03 

33 

t.83036E-01 

3.65237E405 

3.98443E-C 

3.23089E-»03 

7.39098E-»C1 

4.00467E«03 

35 

9.06607E-C: 

4 .22:23E*C5 

5.11269E-C 

3.36CCIE-C3 

7.879C9E*C1 

4.05064E*C3 

36 

9. i8393E-0: 

4 , 56959E*C5 

5  .  837C5E-C 

3.4:795E-»C3 

8.  :'.49CE*0: 

4.064:4E*C3 

P 

LANE  K  -  6 

PH:  -  27.500 

N 

p 

RHC 

- 

V 

N 

5 . 05443E-C; 

4 . 5C897E«C4 

3. : 

4  .  2v6:3E*C2 

-5.59446E*C2 

4  :3632E*C3 

3 

5  27213E-C: 

6. 19923E*04 

4  .  173: 5E-Ci 

•'  .  129F-C2 

-3.39:93t*C2 

4.06663E«C3 

B-9 


NAVSWC  TR  91-86 


5 

5.4t9l3E-01 

7. 76170E+04 

5  07966E-C2 

1.05281E403 

-l.*7350E*C2 

3.967I3E403 

? 

S.70754E-0: 

9.20749E404 

5.99464E-C2 

1 .27470E403 

-6.4***7£«5: 

3.86652E403 

9 

5.92524E'01 

1 . 05993£<»05 

6.75026E-02 

1.43tt4E4C3 

3.12828E4C’. 

3.*0267£«03 

11 

6.14295E-01 

l.i9353£-»05 

7.60629L-G2 

1 .60396E403 

1.18882E-4C2 

3.74C64E403 

13 

6.36055E-01 

I.32654E405 

•.66395E-02 

1.75075E403 

1 .99872E402 

3.72911E403 

15 

5.S7i3SE-01 

1 .45S49C+05 

9.94741E-02 

1.86776E403 

2.6S549E402 

3.75t02E403 

17 

5.79505E-01 

1 .61006£«05 

1 . 11496E-01 

1.99553E403 

^.292^^t*Q2 

3.7397$E403 

19 

7.01375E-01 

1.7€311E-*05 

1 .24023E~01 

2.1233tE403 

3.905e2E402 

3.70505E403 

21 

7.23147E-01 

1.92797E+05 

1 .3933SE-01 

2.25462E403 

4.51951E4C2 

3.59339E403 

23 

7.44917E-01 

2.10942E«05 

l.S95liE-01 

2.3*5*6£403 

5.12511E402 

3.73007E403 

25 

7.55597E-01 

2.31519£<»05 

1.95743E-01 

2.51762E403 

S.71610E402 

3.799**E403 

27 

7.IS4SIE-01 

2.55723E+05 

2. 19139E-01 

2.55474E403 

5.30915E402 

3.t7375e403 

29 

•.10229E'01 

2.l4t24E405 

2.61990E-01 

2.79529E403 

6.90618E4C2 

3.93911E403 

31 

•.31999£~01 

3.205S7E405 

3.19315E-01 

2.93426e403 

7.49S87E402 

3.99957E403 

33 

•.S3759E>01 

3.5S237E+05 

3.9*443£-01 

3.0€S57£«03 

•.0€136£«C2 

4.053I7E403 

35 

•.75539E-01 

4.221i3E-fOS 

S.11269E-C1 

3.1t249£403 

*.St063E-»C2 

4 .10360£-^03 

35 

•.•5424£>0: 

4.S59S9E405 

S.i3704E-01 

3.23641E403 

*.*312*£«S2 

4.11893E403 

PLANE  M  -  12 

PHI  -  57.500 

N 

R 

P 

RHO 

U 

V 

N 

1 

5.04225E-0 

6.34664E404 

4.1070*E-02 

1.04S*9E4C3 

-5.54818E4C2 

3.90466E403 

3 

5.21125£>01 

•.1€2I4B*04 

5.21969E-02 

1.19143E*03 

-2.72905E*C2 

3.877S8E403 

S 

S.3902CB-01 

9.t5473E*04 

<.15427£-C2 

1.3095te4-03 

>5.*9*00E4-01 

3.80134E4D3 

7 

5.54927E-01 

i.l44QSE405 

7.0270*E-02 

1 .41774E*03 

1 .22095E«C2 

3.72462E-03 

9 

::  71I27E-01 

1 .29591E*05 

7.96*93E-02 

1 .49515E403 

2.67083E402 

3.694I1E403 

11 

S.9l72t£*01 

1.444t4E405 

8 

:  .  58199E403 

4.04385£«C2 

3.6847SE403 

13 

5.05629E-01 

2  .  59635E*0.‘i 

’  7.*>;;5E-01 

1 . 64603E«03 

5.17393E*02 

3. 70446E403 

15 

6.22529E-01 

1  .  7502«i’/.4(J:> 

1  i3O*9E-0: 

1.72294E*C3 

6.30603E402 

3.68213E403 

17 

6.39430E'01 

1.90641E4Q:> 

\  .23906E-C1 

1.81116E4Q3 

?.47163E*02 

3.63451E403 

19 

5. 5C330E-C1 

2.0764t£405 

1 ,39327E-0. 

1  .  88523E403 

8.47*20£4C2 

3.6S202E403 

21 

5.73231E-01 

2.26294E40S 

1.59374E-01 

1 . 96009e403 

9.45574E*02 

3.71328E403 

23 

€.90132E>01 

2.47109E+05 

l.*3111E-01 

2.03619£<»03 

1.04083E-»03 

3.78562E403 

25 

7.O7O32E-01 

2.7044SE4Q5 

2 . 103S8E-01 

2.12140E403 

1.140S0E403 

3.t4431E4D3 

27 

7.23933E-01 

2.97344E405 

2.44598E-01 

2.21350£-»C3 

1 .24370E*C3 

3.91111E403 

29 

7.4nt33£-01 

3  29214E*05 

2.i9741E-01 

2.30849E*03 

1 .34771E-03 

3.99205E*03 

31 

7.57734E-01 

3.6tl33E405 

3.49240E-C1 

2.40493E-C3 

1-45292E*C3 

4 . 06*44E<-C3 

33 

7.74634E-0; 

4 . 163:0£4C5 

4.2I576E-0: 

2 .4973iE*:3 

1.5SS53E*03 

4.13355E-»C3 

35 

7.91535E-01 

4.76444E4C5 

5.3807SE-01 

2 . 581*0E-»03 

1.65203E403 

4.19O76E403 

PLANE  M  •  11 

PN:  -  *7.500 

N 

R 

p 

RHC 

U 

V 

N 

1 

S.0302IE-Q1 

1 .23005E405 

7.08036E-02 

1.2«77eE-C3 

-1 .97149E*0C 

3.5I388E*C3 

3 

5.1S139e-0l 

1.43S6?£*0S 

<.31044E>C2 

1.27S54E-*C3 

2.07542E-»C2 

3.610nE*03 

5 

S.27249E-0: 

1 .63061E405 

9.3063SE-02 

1.28304£*03 

3.79843E*C2 

3.55344E403 

7 

5.393S9E-C1 

l.*1443e405 

1.03159E-C: 

1 .29183E403 

5.34lt9E*C2 

3.52405E4C3 

9 

S.5145tE-0l 

1 .99460E*05 

1 . 14567E-C: 

1 .29398E*03 

6.S9440E-02 

3.54041E4C3 

11 

S.€3S7tE-01 

2.17451E*05 

1.2617IE-01 

1.30940E*03 

7.*7416E*02 

3.54576E*C3 

13 

5.7S<«9E>01 

2.35613E405 

1  .38530E-01 

1 .32026£«03 

8.96784E402 

3.S6114E*03 

15 

5.I7799E-C1 

2 . 54830E*05 

1 . 54077E-01 

1.33812E4C3 

1.00956E403 

3.61124E«03 

17 

5.9990IE-01 

2.74I14E«05 

1.72S82E-01 

1.34998E«03 

1.10528E403 

3.68427E*03 

19 

6.12019E-01 

2.96206E<»05 

1  .**<50E-01 

1.37526E*03 

1 .21367E403 

3.68C75E4C3 

21 

6.24129E-0: 

3. lf532E4CS 

2.0S7I4E-C: 

:  .40528£‘C3 

1 .32622E*C3 

3.66€:6£.C3 

23 

6.36239E-01 

3.43313E*05 

2.30763E-0: 

1 .43249E-*03 

1 .4324fE-‘03 

3.716?9E*C3 

25 

6.49349E-01 

3 . 70979E*0S 

2 .64277E-01 

1 .46235£*C3 

1.53979E-C3 

3.e:il9E*03 

27 

6.60459E-01 

4.03127E+C5 

3 . 05099E-01 

1 .49569E403 

:.6497SE*C3 

3.91289E*C3 

29 

6. 72569E-01 

4.40569E405 

3  .  537*3E-C: 

1 .53265E-»03 

1.75422E*03 

4.0C332E-C3 

31 

6.t4679E-0l 

4.*4*61E-*05 

4.147C1E'G: 

2  .  S7193E*C3 

1.88136£*:3 

4 .09C26E*C3 

33 

6.967I9E-0: 

5.3i671E*C5 

4 .94:72E*Ci 

1  .  6ii35E*C3 

1 .9966CL‘C3 

4.17519E4C3 

35 

7.09f99E-0: 

6.03931E405 

5.97105E-C1 

1 . 64937E*C3 

2.109467.»03 

4.251I1E*03 

36 

7 . 14954E-01 

6.43207E405 

6.63077E-.:: 

:  .66419£*03 

2.16I15E*03 

4.29:93E*03 

PLANE  N  -  24 

pk:  -  117. 50C 

N 

R 

p 

RHO 

U 

V 

N 

1 

5.02176E-01 

2.28381E405 

1 .  19312E-01 

9.28784E*C2 

5. :9228E4Cr 

3.34506E4C3 

3 

5.10878E-01 

2.50513E405 

1  .32477E-01 

8 .662*7e.C2 

6.26r76E.C2 

3.40472E*C3 

5 

5.195I0E-01 

2.72160E*05 

1.42a05E-01 

8.24651E«'02 

7.28943E*C2 

3 . 369i7E*03 

7 

5.282I2E-01 

2.93069E405 

1 . 54380E*C1 

7.84388E402 

8.22790E*C2 

3.38072E-.C3 

9 

5.369I4E<01 

3.13836E4C;5 

1 .6663iE-C: 

7. 5C668E402 

9.063C7E.G2 

3.40409E*03 

11 

5.456*€E-0' 

3.34619r*r« 

1  '»9596E-0: 

7.ri9i:E-C2 

9.97i:2E*C2 

3.42983E*C3 

*3 

S .  343*i£ 

3.55826E405 

\:  .95192E-C1 

6.93737E-»02 

i  .07126E*C3 

3,49630E*C3 

15 

5.63090E-01 

3.7I247E40S 

2 . 10722E-C1 

6. 71875E^C2 

1  .  15772E*03 

3.52805E*C3 

17 

5.71792E-01 

4.01063E*05 

2 .2t749L-Cl 

6.49778E*02 

1 .23*95E*C3 

3 . 58100E«03 

20 

5.84845E-0I 

4.38t41E405 

2 . 67750E-C: 

6. 17791F-C2 

1.357C4E*C3 

3.74900E-C3 

21 

5.89196E-01 

4.52474E405 

2.eC733E-Ci 

6.0940CE*02 

1 .39597E*03 

3.7|C01E*C3 

23 

5.97898E-01 

4.I1672E405 

3.04456E-01 

5.94946E*C2 

1 .48750E*03 

3.79640E-03 

25 

€.0<600E-01 

5.13473E405 

3.35090E-01 

5 . 80409E402 

1.57746E4C3 

3.I4300E4C3 

27 

6. 15302E-01 

5.49507E*05 

3 . 76422E-0: 

5 . 65184E*D? 

1 .66594E*C3 

3.92522E*03 

29 

6.24004E-01 

5.90424E405 

4 .26995E-01 

5.  5i83:E-»02 

1.7597CE‘03 

4  .C2243E-»C3 

31 

6.32706E-01 

6.38042L405 

4.8939CE-C: 

5 . 4C513E*02 

1.85747E403 

4 . 12672E*C3 

33 

6.4140IE'01 

6.94979E405 

5 . 688C;E-C: 

5.294C7E*C2 

1 .956S6E*C3 

4.23I42E«C3 

35 

6. 5On0E-Ol 

7.62267E*05 

6  69329E-01 

5  .  17479E*C2 

2  C5:89E*C3 

4 .34954E*03 

36 

6. S4461E-01 

•.01766E*05 

7 .26968E-01 

5.I06S8E«02 

2.09835E*C3 

4.39239E4C3 

PLANE  M  «  30  PHI  -  147, 5CC 


B-10 


NAVSWC  TR  91-86 


N 

R 

P 

RHO 

U 

V 

N 

1 

5.01703E-01 

4 .0S408E405 

1.95198E-01 

3.C4170E402 

S.13634E4C2 

3.C7553E4C3 

3 

5.09516E-01 

4.30444E-^05 

2.08815E-01 

2.40892E4^02 

5-57370E402 

3.16573E403 

5 

5.1S329E>01 

4.52607E-»05 

2.19100E-01 

1.85903E402 

6.Q6221E402 

3.15583E403 

7 

5.22142e-01 

4.74431E*05 

2.31443E-01 

1.32243E«02 

6.50613E402 

3.19875E403 

9 

5.2i954E-01 

4 .96483E405 

2.44829E-01 

8.42503E401 

6.93089E402 

3.25552E403 

11 

5.35767E-01 

5.19221E-»05 

2.58831E-01 

3.64799E401 

7.3891SE402 

3.30857E403 

13 

5.425IOE*'01 

5.42419E-^05 

2.73770E-01 

-7.84019E400 

7.80442E402 

3.36065E403 

15 

5.49393E>01 

5.66474E405 

2.92048e-0l 

-5.54704E401 

8.27221B402 

3.44004E403 

17 

5.56206E-01 

5.91772E+05 

3.11717E-01 

>9.67004E401 

8.69075E402 

3.51383E403 

19 

5.(3019E-01 

6.18997E+05 

3.32967B-01 

-1 .41673E402 

9.19302E402 

3.S7623B403 

21 

S.69I31B-01 

6.48039E405 

3.63394E>01 

<1.88979E402 

9.65733E402 

3.69532E403 

23 

S.76644E>01 

6.79971E+05 

3.95218E-01 

-2.33293E402 

1.01254E403 

3.7B615E403 

25 

5.93457e>01 

7.14208E-f05 

4.22519E-01 

-2.75813E402 

1.06301E403 

3.82090E403 

27 

5.90270E>01 

7.S1810E405 

4.S9481E-01 

-3.21601E402 

1.11676E403 

3.89004E403 

29 

5.97042E-01 

7.94100E405 

5.13409E-01 

-3.68408E«02 

1.17146E403 

4 .00I91E403 

31 

6.03995e-01 

8.42786E405 

5.76693E-01 

-4.15e74E*02 

1.22720E403 

4.12177E403 

33 

6.10704E-01 

8.99287E-f05 

6.50744E-01 

-4.63<21E4C2 

I.28481E403 

4.23323E403 

35 

6.17521E-01 

9.64052E405 

7.40231E-0I 

-*.n99QE-*C2 

1.34240E403 

4.34964E403 

36 

6,20927E-01 

1.00107E+06 

?.9119#£--C1 

-5.36t50£4C2 

1.37190E403 

4.40335E403 

FLAKE  M  -  36 

PHI  -  177. 50C 

M 

R 

P 

RHO 

U 

V 

M 

1 

5.01536E-01 

4.08408E405 

1.95198E-01 

-1.31155E40: 

4 .99661E401 

3.13250E*03 

3 

S.076I2E-01 

4.30444E+05 

2.08ei5£-01 

«7.377H£.o: 

5.28534E401 

3.22216E«03 

S 

5.13428E-01 

4.52607E-4-05 

2.19100E-0: 

-1 .32389E402 

5.63776E401 

3.21569E403 

7 

5.19974E>01 

4.74431E405 

2.31443E-C1 

-•1.88779E402 

5.95523E401 

3.26092E403 

9 

5.26120E-01 

4.96463E«05 

2.44829E-01 

-2.41482E*02 

6.26842E401 

3.32019E*03 

11 

5.32266E-01 

5.19221E405 

2.5a831E-01 

-2.96377E»02 

6.615S8E401 

3.3766SE403 

13 

5.34412E-01 

5.42419£-f05 

2.73770E-01 

-3.47417£»02 

6.93123E401 

3.43185E403 

15 

5.445S8E-01 

S.66474E40S 

2.92048E-01 

-4.01517E402 

7.28219E401 

3.S1493E403 

17 

5.S0704E-01 

5.91772E40S 

3.11717E-01 

-4 . S2259E402 

7.61256E401 

3.59184E403 

19 

5.56850E-01 

6. 18997E-*05 

3.32967E-01 

-S.09706E402 

8.01403E401 

3.65901E403 

21 

5.62996£'01 

6.48039£*05 

3.63394E-01 

-5.59717E402 

8.34615E401 

3.78200E403 

23 

5.69141E-01 

6.79971E405 

3.95219E-Q1 

-6.13079E402 

8.70884E401 

3.I7700E403 

25 

5.75287E-01 

7.142O0E*O5 

4.22519E-01 

>6.71988E402 

9. 12928E401 

3.91732E403 

27 

5.81433E-01 

7. 51ll0E-»05 

4.59481£<01 

-7,31551E402 

9.  56097E401 

3.99232E*D3 

29 

5.87579E-01 

7.94100E*05 

5.13409E-01 

-7.91341E402 

9.99577E401 

4.11621E403 

31 

5.93725E-01 

8.42786£<»05 

5.76695E-01 

-8. 53208E402 

1.04475E402 

4.23420E403 

33 

S.99871E-01 

8.99286E^05 

6.50744E-01 

-9. 16009E402 

1.09079E402 

4.35143E403 

35 

6.06017e-01 

9.640S2E'»05 

7.40230E-01 

>9. 78731E402 

1.13688E402 

4.47356E*03 

36 

6.09090E-01 

1.00107E+06 

7.9119eE-0l 

-1 .01152E403 

1.16073E402 

4.53025E403 

1 SHOCK 

DATA  FOR  ZEUS 

M 

C 

cz 

CPHl 

PHI 

1 

9.24286E-01 

9.61111E-01 

-1 .41036E-02 

2 . SOOOOE^OO 

3 

9.17666C>01 

9.48112E>01 

-6. 85406£*02 

1 .25000E40I 

5 

9.02486E-01 

9.18274E-01 

-1.15430E-01 

2.25000E401 

7 

8.79361E-01 

e.72594E*01 

-1 . S0059E-01 

3.2S000E*C: 

9 

8.50172E-01 

8.14522E-01 

-1.7027eE-01 

4.25000E401 

11 

8.20553E-01 

7.56569E-01 

-1 .79262E-01 

5.25000E401 

13 

7.88659E-01 

6.93760E-01 

-1 .76687E-01 

6.25000E*01 

15 

7, 58e57E-01 

6.33946E-01 

-1 .66993E-01 

7.25000E*0: 

17 

7.31285E-01 

5.79489E-0: 

-1 .52921E-01 

8.25000E*01 

19 

7.05739E-01 

5.28376E-01 

-1 .35594E-01 

9.25000E*0: 

21 

6.83649E-01 

4.82485E-0: 

-1.17210E-01 

1.0250CE402 

23 

6.65103E>01 

4 .42779E-01 

-9.9027eE-C2 

1 . 12500E*02 

25 

6.49050E-01 

4.07447E-0: 

-8.11792E-02 

1.22500E402 

27 

6.36596E-01 

3.79871E-01 

-6.48930E-02 

1.32500E*02 

29 

6.26445E-01 

3.56429E'01 

-4 .94732E-02 

1.42500E402 

31 

6. 19289E-01 

3.39007E-01 

-3.5283iE-02 

1.52500E402 

33 

6.14146E-01 

3.26190E-01 

-2.19254E-02 

1 .62500E*C2 

35 

6. 11218E-01 

3.18767E-01 

-9.25615E-03 

1.72500E*02 

36 

6.10626E-01 

3.17255E-01 

-3.07567E-03 

1 . 77S00E*C2 

USER 

TIME  -  305, 

.0160 

SYSTEM  TIME  -  1.236340 

FORTRAN  STOP 


B.2.2  ZEUS  Output 


1  PR0CRAM2EUS 


B-n 


NAVSWC  TR  91-86 


••••  IKTEGRATION  CONTROL 


I  CORD 

- 

1 

ZETA 

- 

o.soooo 

ZETAEND 

- 

10.00000 

KSTART 

- 

0 

KEND 

• 

500C 

FCFL 

> 

0.90000 

XKI 

• 

l.OOOOCO 

lAPR 

■ 

0 

1 SHF LUX 

- 

0 

ISHANG 

• 

1 

(0-CANT£SlAN.  I -CYLIKt>RXCAl  COORDINATES* 
(INITIAL  ZETA) 

(MAXIMUM  £ETA» 

(INITIAL  STEP  NUMBER* 

(MAXIMUM  STEP  NUMBER* 

(CFL  SAFETY  FACTOR} 

(INTERIOR  POINT  LIMITER) 

(O-APPROXIMATE  RIEMANN  PROBLEM 
INCOMPLETE  RIEMANN  PROBLEM) 

(FLUX  AT  SHOCK-  0-RIEMANN  PROBLEM 

l-FREESTREAM  PROPERTIES) 
(SHOCK  ANGLE  -  0-RIEMANN  PROBLEM 

1-ROE  AVERAGED  VARIABLES) 


BOUNDARY  CONDITIONS  •••• 


KXKE 

NSUR 

DFAC 

ISHOCK  - 
lASYM 


(NUMBER  OF  EDGES  NITH  NON-DEFAULT  LIMITER  SETTINGS) 
(NUMBER  OF  EDGES  NITH  NON-DEFAULT  SURFACE  TYPE) 
(REDUCE  TURNS  NHICH  ARE  TOO  LARGE  BY  DFAC) 

(1-SHOCK.  0-NC  SHOCK  ON  EDGE  3) 

(O-NO  SYMMETRY. \-PnCH  PLANE  SYMMETRY) 


SURFACE  POINT  LIN.ITER 


ZONF.  EDGE 

1  l.OCC 


EDGE 

O.OOC 


EDGE  3 
2.00C 


EDGE 

O.OOC 


SURFACE  TYPE 

ZONE  EDGE  1 

1  1 


EDGE  2 
0 


EDGE  3 
1 


EDGE  4 
0 


•••  OUTPUT  CONTROL  •• 


PRINTER  (fil« 
IPRIN?  •  2SC 

ISKJP  -  10 

NSKIP  •  S 

MSKIP  -  < 


(PRINT  FLOW  FIELD  IF  STEP  NUMBER  IS  EVENLY  DIVISIBLE  BY  IPRIKT) 
(PRINT  STEP  SIZE  IF  STEP  NUMBER  IS  EVENLY  DIVISIBLE  BY  ISKIP) 
(PRINT  N  PLANES  NHICH  ARE  EVENLY  DIVISIBLE  BY  NSKIP) 

(PRINT  M  PLANES  NHICH  ARE  EVENLY  DIVISIBLE  BY  MSKIP) 


PRESSURE  SUMMARY  (PLOTZA) 

KSKIP  «  40  (PRINT  SURFACE  PRESSURE  IF  STEP  NUMBER  IS  EVENLY  DIVISIBLE  BY  KSKIP) 

NPRT  •  ICC  (DON'T  PRINT  POINTS  IF  N>KPRT) 

DELZA  iOOOO.COC  (PRINT  SURFACE  PRESSURE  EVERY  ZETA  -  CELZA* 

JSPPR  “  0  'NUMBER  OF  NON-DEFAULT  EDCES) 

JSPDIS  -  0  (PRINT  SURFACE  GRID  POINT  LOCATIONS) 

SURFACE  PRESSURE  OUTPUT  CONTROL-  O-DON'T  PR  I  NT. 1 -PR INT. 

ZONE  EDGE  1  EDGE  2  EDGE  3  EDGE  4 

1  :  0  C  0 


PLOT  FI 

I.PLCT 

DELZC 

IPLOTN 

ZPLOT 


(PLOTZC) 

IQOCC  WRITE  TC  PLC72C  IF  STE.»  NUMBER  IS  E'/ENIY  DIVISIBLE  BY  IPLCT 

lOOCO.OCC  (NPITE  PLOT  TAPE  EVERY  ZETA  -  DElZC). 

2  NUMBER  OF  TARGET  PLOT  STATIONS 

O.SCCLC  lO.CCCOC. 


•••••AERODYNAMICS  D 
AREF  -  I.OOCC 

alref  -  l.OCOC 

lEFORCE  -  C 


ATA***** 

(REFERENCE  AREA) 

(REFERENCE  LENGTH) 

<C  -  DONT.;  -  DO  PRINT  EDGE  LOADS) 


•••••SEPARATION  MODELI NC» • • • • 

jVIS  -  C  (C  -  INVISCID.  1  '  CLIPPING.  2  -  FORCED  SEPARATION) 


CAS  TYPE 


B-12 


NAVSWC  TR  91-86 


1GA5  -  1  (O-PERFEG*!-  '  EQUILIBRIUM  AIR) 


••••••UNITS*-**** 

XUNITS  -  1  <0-ENGLlSH.  1-KETRlC) 


0 


••••  INITIAL  DATA  PLANE  •••• 


IMOD 

MACH 

ALPHA 

BETA 

PINT 

DIKF 

HO 

SINF 

IZN 

NA 

MR 

HA/ZOKE: 
IMESHF  - 

IMESHG  • 


1  <0  SELF  START: 1  PROM  TAPE  START) 

20.0000 (FREE'STREAM  MACH  NUMBER) 
20.0000<ANGLE  OF  ATTACK) 

O.OOOOIYAM  ANCLE) 

SS29.I779(FREE>STREAM  PRESSURE) 
•.•915E-02<rREE-8TR£AM  DENSITY) 
}.5162E<»12(STAGNATIOK  ENTHALPY) 
O.OOOOiFREE-STREAM  ENTROPY) 


I 

3$ 

36 

3C 

0 


(NUMBER  OF  ZONES) 

(CELLS  BETKEEN  EDGES  1  AND  3) 
(TOTAL  CELLS  BETWEEN  EDGES  2  AND 

(O-UNIFORM  MESK. 1-CLUSTERED  MESH 
BETWEEN  EDGES  1  AND  3} 

<0-UNZFORM. 1-CLUSTERED  MESH 
BETWEEN  EDGES  2  AND  4) 


4) 


1  STEP  -  0  ZETA  •  5.0000000E-01  D2ETA  -  0 . OOOOOCOE+00 

MACH  NO.  -  2.0000000E401  ALPHA  -  2 . OOOOOOOE-Ol  BETA  •  0 . 0000000£*00 

0 

0  ZONE  NO.  1 


CELL  PLANE  M  -  1 

C  •  9.242t60lE-0I  C2  • 

B  -  5.0000000E-01  B2  - 


9.6U1101E-01  CT  -  -1 .4103631E-02  VX  - 
O.OOOOOOOC-00  BT  >  0 . 0000000E-»00  VX  • 


9.9S3a262IE-01  5 . •a638149E-02 

O.OOOOOOOOE«00  O.OOOOOOOOEtOO 


0.00000000£«00  ISUR  3 
O.OOOOOOOOE400  ISUR  1 


N 

S 

T 

P 

D 

U 

V 

W 

E 

M 

36.3 

9.1I3931E-01 

4.363323E-02 

4 . 5€9594£«05 

S.I37047E-01 

3.420235E*03 

-6.810439£^01 

4.06413IE^03 

I.944247E403 

S.566885E400 

31.5 

I.59464SE-01 

4.363323E-02 

3.205S6#£>05 

3.193162E-01 

3. 096449E-03 

-6-60»672E^Cl 

3.954461E^03 

9.694987E403 

4 .497383E400 

23.3 

7.SI7302E-01 

4.363323E-02 

2.31S177E*05 

1  .•57429E-01 

2.625004e*03 

-6.278110E*01 

3.769571E403 

1.027318e404 

3.700651E400 

19.3 

7.1I0339C-01 

4.363323E-02 

1.763107E*05 

1.240225E-Q1 

2.177979E+03 

-6.071832E«01 

3.693437E^03 

1.063360E404 

3.276337E400 

13.3 

6.47321SE-01 

4.363323E-02 

1 .326544E-05 

I.6639S0E-02 

1.7554S3E*03 

-6.044731E*01 

3.731761E*03 

1 .D854i9£404 

3.085242E400 

7.5 

3.766072e-01 

4.363323E-02 

9.207491E404 

5 . t946€3E-02 

2  .  1<3502£-»03 

-6. 188703B+01 

3.»95451E^03 

2 . 2029D3E404 

3.029335E400 

1.3 

3.05S929C-01 

4.363323C-02 

4.50<974E*04 

3. 1147i5E-02 

1 .251147E*02 

-6.613063E*0l 

4.192732E^03 

1 .  n4449E404 

3.218029E40C 

CELL 

PLANE  H  • 

12 

C  - 

8.0421042E-01 

C2  -  7.2419li2E-01  CT 

-  -1.7874149E 

-01  VX  -  3 

,49807t58E-01 

9.62836862E-C1 

O.OOOOOOOOE^OO  ISUR  3 

B  - 

5.0000000E-01 

B2  *  0 . OOOOOOOE-00  BT 

-  O.OOOOOOOE 

♦00  VX  -  0 

.OOOOOOOOE*O0 

O.OOOOOOOOE*OC 

O.OOOOOOOOE^OO  ISUR  1 

N 

S 

T 

P 

D 

U 

V 

w 

E 

M 

36.3 

7.999I53E-01 

1 .003564E40C 

5.131466E405 

6.04063eE-01 

2.840194E^03 

-1 .2973t6E*03 

4.212109E*03 

9.046637E403 

5.323775E*00 

31.5 

7.57733SE-01 

1 .003564E*0C 

3.6ai330E«05 

3.492402E-01 

2.517548E+03 

-1 .247642E*03 

4.06I437E*03 

9.747701E*03 

4 .32I0S6E400 

25.5 

7.070321E-01 

1 . 003564E-fOO 

2.704453£«05 

2 . 103579E-01 

2.10171SE*03 

-1 .  178383E^03 

3.844306E«03 

1.030159E404 

3 . 586752E4C0 

19.5 

6.563303E-01 

I.003564E-00 

2.076475E>05 

1 .393267E-01 

1.727974E*03 

-1 . 134454E+03 

3.652021E^03 

1 .064846E404 

3. 157021E*00 

13.5 

6.0S62t6E-01 

1 .003564E-00 

1.5963S2E*05 

1 .019252E-01 

1 .320608F«03 

-1 .110360E^03 

3.704460E^03 

1 . 0B3014E404 

3.029101E400 

7.5 

5.549269E-01 

1 .003564E>00 

1 .144051E^05 

7.027075E-02 

I.647255E*02 

-1.130109E^03 

3.724623E^03 

1 . 103635E404 

2.9165S4E4C0 

1.5 

5.042251E-01 

1 .003564E*00 

6.346638E*04 

4  . 107Oi0E-02 

•.559204E*01 

-1.185570E*03 

3.904660E^03 

1 . 11465SE404 

3.D53542E400 

CELL 

PLANE  M  • 

24 

C  - 

6.S6636t9E-01 

C2  -  4.242109ie-01  CT 

-  -I,9929541E 

-02  VX  -  -5 

.83229065E-01 

8.23772132E-01 

0 . OOOOOOOOE 

♦00  ISUR  3 

B  - 

5.0000000E-01 

B2  -  0 . 0000000E*0C  BT 

-  O.OOOOOOOE 

♦00  VX  -  0 

.ooooooooE*oo 

0  OODOOOOOE400 

D.OOOOOOOOE 

♦CC  ISUR  1 

N 

s 

T 

P 

D 

u 

V 

W 

E 

M 

36.5 

6.544614E-01 

2.050762E*00 

•.017663E*05 

7.2696I2E'01 

1 .62S465£«03 

-1 .421870E*03 

4 . 3923I7E*D3 

9.486033E403 

4 .264193E400 

31 . 5 

6.327063E-01 

2.050762E*00 

6.3I0416E405 

4.I93I99E-01 

1 .398C19EtC3 

-1 .337128E^03 

4.1267nE^03 

9.9913:2E403 

3.595S01E*00 

25.5 

6.066001E-01 

2.050762E-00 

5.134731E+05 

3.350902E-01 

1 . 131219E+03 

-1 .243218E«03 

3 . •43C02E^03 

1 . 037733E4C4 

3.2C3972E*00 

19.3 

5.t04939E-01 

2.050762E-*00 

4.257402E«05 

2  .  540733E-01 

•  .791S66E-»02 

-1 . 165272E*03 

3.695B23E^03 

1.06Q195E*04 

2.856104E40C 

13.5 

3.543I7IE-01 

2.050762E*00 

3.55I261C*05 

1.951924e-01 

6.298855E>02 

-1 .110004E^03 

3.496295E^03 

I .086:38E404 

2.584120E400 

7.3 

3.2I2I17E-01 

2.050762E*00 

2.930691E«05 

1.543I02E-C; 

3.676338E*02 

-1 .0756I3E*03 

3.380723E*03 

1 . 104831E404 

2.431878E400 

1.3 

5.021755E-01 

2.050762E-»00 

2.2I3I14£«05 

1 . 193U4E-01 

3.169659E^0: 

-1.063594E^03 

3.345056E«03 

1 . 1 :7764E*04 

2.3i47SlE400 

CELL 

PLANE  M  - 

36 

C  - 

6. 106262IE-01 

C2  -  3. 1725460E-01  CT 

-  -3.075672IE 

-03  VX  -  -9 

.99267936E-01 

3.I5I71679E-02 

O.OOCOOOOOE400  ISUR  3 

B  - 

5.oooooooe-oi 

B2  -  O.OOOOOOOE*OC  BT 

•  O.OOOOOOOE 

*00  VX  -  0 

.OOOOOOOOE^OC 

O.OOOCCOOOE^OC 

0. OOOOOOOOE 

♦00  ISUR  1 

N 

s 

T 

P 

D 

u 

V 

W 

E 

H 

36.5 

6.090999E-0: 

3. 097960E*00 

1 .0C:C66£*06 

7.9:i97iE-L; 

1.015618t^03 

-7.1I4073E^CI 

4.53C247E*C3 

9.74I788E*03 

3.726900E*00 

31.5 

5.937250E-01 

3.097960E*OC 

•  .42‘>|79E+05 

5.766953E-C1 

e.569534E^02 

-6.71S912E^01 

4 .2342D2E«03 

1 . 012967E404 

3.234094E40C 

25.5 

5.752I73E-0: 

3.097960E*0C 

7 . 142076E*05 

4 .2251I5E-C1 

6.753306E^C2 

-6.1I9426E«C1 

3 . 9n32iE^03 

1 .053936E*04 

2 . 778456E400 

19.5 

5.568496E-01 

3.097960E*0C 

6. 149971E*05 

3.329672E-0: 

5.127:66E*02 

-5.783103E*0: 

3.659011E^03 

1 .06I942E*C4 

2 . 53t745E4C0 

13.5 

5.3f4119C-01 

3.09796CC*0C 

5.4241#9E405 

2. 737695E-C: 

3.5C:2COE*C2 

-5.40922:E*01 

3. 43:i50E*D> 

1 . D90889E*04 

2.302670E*C0 

7.5 

5 . 199742E-01 

3.097960E-00 

4 . 7443C9E-C5 

2  .  314432E>C: 

1 , 9:i973E*C2 

-5.126116E*CI 

3.260924E»C3 

1 . 106457E*04 

2 . 144908E40C 

l.S 

5.0:5365E-0; 

3.097960E*00 

4 . OI4C7aE«C5 

1  .  951976E-01 

1 . 528248£^C: 

-4.934645E*C; 

3 . 1325C:E*C3 

I . 11939l£«04 

2.035921E*00 

B-13 
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SECANT 

METHOD 

JN  DECODE 

.DID  NOT  CONVERGE  FOR 

M. 

N 

1  34 

SECANT 

METHOD 

IN  DECODE 

.  DID  NOT  CONVER 

SE  FOR 

M. 

N 

2  34 

OSTEP 

- 

20 

EETA  - 

6.10449910£>01 

B2ETA 

- 

6.26I37742E-03 

MCFL.MCFL 

- 

1 

35 

SECANT 

METHOD 

IN  DECODE 

.DID  NOT  CONVERGE  FOR 

M, 

N 

14  35 

SECANT 

METHOD 

IN  DECODE  DID  NOT  CONVERGE  FOR 

M, 

N 

17  36 

OSTEP 

- 

30 

EETA  - 

6.771507I6E-01 

D2ETA 

- 

7.00617116E-03 

NCFL.MCFL 

• 

1 

3d 

SECANT 

METHOD 

IN  DECODE 

;  DID  NOT  CONVERGE  FOR 

M. 

N 

15  34 

SECANT 

METHOD 

ZN  DECODE  DIO  NOT  CONVERGE  FOR 

M. 

N 

14  36 

OSTEP 

• 

40 

EETA  - 

7.51731277E-01 

D2ETA 

- 

7.77146500E>C3 

NCPL.MCFL 

• 

1 

36 

SECANT 

mntoo 

ZN  DECODE  DID  NOT  CONVERGE  FOR 

M. 

N 

14  34 

SECANT 

METHOD 

IN  DECODE  DID  NOT  C^fVERGE  FOR 

M. 

N 

20  35 

OSTEP 

- 

SO 

EETA  - 

•.33479047E>01 

DZETA 

- 

4.54412373E>03 

NCFL.MCFL 

- 

1 

35 

SECANT 

METHOD 

IN  DECODE 

.  DID  NOT  CONVERGE  FOR 

M. 

N 

29  36 

SECANT 

METHOD 

IN  DECODE 

.  DID  NOT  CONVERGE  FOR 

M. 

N 

32  35 

OSTEP 

- 

40 

EETA  - 

9.2299tl30E-01 

DZETA 

- 

9.22467t39E-03 

NCFL.MCFL 

• 

1 

36 

SECANT 

METHOD 

IN  DECODE 

,  DID  NOT  CONVERGE  FOR 

M. 

N 

27  34 

SECANT 

METHOD 

ZN  DECODE  DID  NOT  CONVERGE  FOR 

M. 

N 

24  34 

OSTEP 

- 

70 

EETA  - 

1.01756394E400 

DZETA 

- 

9.61779151E-03 

NCFL.MCFL 

- 

1 

36 

OSTEP 

• 

•  0 

EETA  - 

1.115i3364E-^00 

DZETA 

- 

1.00026354E-02 

NCFL.MCFL 

- 

1 

35 

OSTEP 

- 

90 

EETA  - 

1.21771550E+00 

DZETA 

- 

1.03135379E-02 

NCFL.MCFL 

- 

1 

36 

OSTEP 

- 

100 

EETA  - 

1.321t6139£40C 

DZETA 

- 

:.04742:53E-0? 

MCFL.MCFL 

- 

1 

36 

OSTEP 

- 

110 

EETA  - 

1 .4271194 9E4QC 

DZETA 

- 

1 .05639524E-02 

NCFL.MCFL 

• 

1 

35 

OSTEP 

• 

120 

EETA  - 

1 . 53316963E-f00 

DZETA 

- 

1 .06327394£<02 

NCFL.MCFL 

- 

1 

35 

OSTEP 

- 

130 

EETA  • 

1.639703i7£400 

DZETA 

- 

1 .06570404E-02 

NCFL.MCFL 

- 

1 

36 

OSTEP 

- 

140 

EETA  - 

1.74617064E40Q 

DZETA 

- 

1.06357164£>02 

NCFL.MCFL 

1 

35 

OSTEP 

- 

150 

EETA  - 

l.l522403IE-*'00 

DZETA 

- 

1.05427354E-02 

NCFL.MCFL 

- 

1 

35 

OSTEP 

- 

160 

EETA  - 

1.95765221E4>00 

DZETA 

- 

1 . 049676442-02 

NCFL.MCFL 

- 

1 

35 

OSTEP 

• 

170 

EETA  - 

2.06214452E400 

DZETA 

- 

:.04126129E-D2 

NCFL.MCFL 

- 

1 

35 

OSTEP 

- 

110 

EETA  - 

2.16567969E40C 

DZETA 

- 

1.0297b9S5E-02 

NCFL.MCFL 

• 

1 

35 

OSTEP 

- 

190 

EETA  • 

2.26793265E-^00 

DZETA 

- 

1 .01717506E-02 

NCFL.MCFL 

- 

1 

35 

OSTEP 

- 

200 

EETA  - 

2 .36t94994E-»00 

DZETA 

- 

1 .00452143E-02 

NCFL.MCFL 

• 

1 

36 

OSTEP 

• 

210 

EETA  - 

2.46I55021E-»QO 

DZETA 

- 

9.44014699E-03 

NCFL.MCFL 

• 

1 

36 

OSTEP 

- 

220 

EETA  - 

2.56644564E-»00 

DZETA 

- 

9. 71429426E-03 

NCFL.MCFL 

1 

36 

OSTEP 

- 

230 

EETA  - 

2.66279626E400 

DZETA 

- 

9. 56043440E-03 

NCFL.MCFL 

- 

1 

35 

OSTEP 

- 

240 

EETA  - 

2.75750709E400 

DZETA 

- 

9.3974499SE-03 

NCFL.MCFL 

- 

1 

36 

OSTEP 

- 

250 

EETA  • 

2.I5059977E«00 

DZETA 

- 

9.2446;i49£-03 

NCFL.MCFL 

• 

1 

36 

OSTEP 

- 

260 

EETA  - 

2.94215294E40Q 

DZETA 

- 

9.0I44324SE-03 

NCFL.MCFL 

• 

1 

36 

OSTEP 

- 

270 

EETA  - 

3.03217936E400 

DZETA 

- 

4.93734307E-03 

NCFL.MCFL 

- 

1 

36 

OSTEP 

- 

210 

EETA  - 

3.12Q64349E400 

DZETA 

- 

4.41141145E-03 

NCFL.MCFL 

- 

1 

36 

OSTEP 

- 

290 

EETA  - 

3.20t43l63E400 

DZETA 

- 

4.72C4S569E-03 

NCFL.MCFL 

• 

1 

36 

OSTEP 

- 

300 

EETA  - 

3.29529333E*C0 

DZETA 

- 

4.66061724E-03 

NCFL.MCFL 

• 

1 

36 

OSTEP 

• 

310 

EETA  - 

.3.3li713«7E-C0 

DZETA 

• 

8.62495790E-C3 

NCFL.MCFL 

- 

1 

36 

OSTEP 

- 

320 

2ETA  - 

3.467|99l0E-*00 

DZETA 

- 

4.60946672E-03 

NCFL.MCFL 

1 

36 

OSTEP 

• 

330 

EETA  • 

3. 5539264 7E400 

DZETA 

• 

4  .  S9935209E-03 

NCFL.MCFL 

• 

1 

36 

OSTEP 

- 

340 

EETA  - 

3.63992476e400 

DZETA 

- 

4 . S99256:7E-03 

NCFL.MCFL 

• 

1 

36 

OSTEP 

• 

350 

EETA  - 

3.7254536IE40C 

DZETA 

- 

4.6005954:e-C3 

NCFL.MCFL 

• 

1 

36 

OSTEP 

• 

360 

EETA  - 

3.4119471ie«00 

DZETA 

- 

4.63548:42E-C3 

NCFL.MCFL 

- 

1 

36 

OSTEP 

• 

370 

EETA  - 

3.49I69714E*0C 

DZETA 

• 

4.  72740336E-03 

NCFL.MCFL 

• 

1 

36 

OSTEP 

- 

340 

EETA  - 

3.9I677015E*DO 

DZETA 

- 

4.44944134E-03 

NCFL.MCFL 

- 

1 

36 

SECANT 

METHOD 

ZN  DECODE 

.  DIO  NOT  CONVERGE  FOR 

M. 

N 

36  7 

OSTEP 

• 

390 

EETA  - 

4.0767C307E*00 

DZETA 

• 

9.04970740E-03 

NCFL.MCFL 

• 

1 

36 

SECANT 

METHOD 

IN  DECODE 

DIO  NOT  CONVERGE  FOR 

H. 

N 

36  7 

OSTEP 

- 

400 

EETA  - 

4.16t9171IE400 

DZETA 

• 

9.32439606E-03 

NCFL.MCFL 

.. 

1 

36 

OSTEP 

- 

410 

EETA  - 

4 .2636«237e40C 

DZETA 

• 

9.54475247E-C3 

NCFL.MCFL 

• 

1 

36 

SECANT 

METHOD 

IN  DECODE 

;  DID  NOT  CONVERGE  FOR 

M.N 

35  6 

SECANT 

METHOD 

IN  DECODE 

DID  NOT  CONVERGE  FOR 

M. 

H 

36  6 

OSTEP 

- 

420 

EETA  - 

4 .360921I6E40C 

DZETA 

- 

9. 42640963E-O3 

NCFL.MCFL 

• 

1 

36 

OSTEP 

- 

430 

EETA  - 

4 .46059604E40C 

DZETA 

- 

1.00939022E-C2 

NCFL.MCFL 

• 

1 

36 

OSTEP 

- 

440 

EETA  - 

4.5633320IE400 

DZETA 

• 

1  .  04226097E-02 

NCFL.MCFL 

• 

; 

36 

OSTEP 

- 

450 

EETA  - 

4.66950130C40C 

DZETA 

- 

1.07S633C3E-02 

NCFL.MCFL 

• 

1 

36 

SECANT 

METHOD 

IN  DECODE 

DID  NOT  CONVERGE  FOR 

M. 

N 

35  5 

SECANT 

METHOD 

IN  DECODE 

DID  NOT  CONVERGE  FOR 

M. 

y 

35  5 

SECANT 

METHOD 

IN  DECODE 

DID  NOT  CONVERGE  FOR 

H. 

N 

36  5 

OSTEP 

- 

460 

EETA  - 

4 . 77955532E400 

DZETA 

- 

l.il9;73ltE-02 

NCFL.MCFL 

• 

1 

35 

OSTEP 

- 

470 

EETA  - 

4.I9393473E-*0C 

DZETA 

- 

1 . 164C79C4E-02 

NCFL.MCFL 

1 

35 

OSTEP 

- 

490 

EETA  - 

5.C:244700E-*CC 

DZETA 

- 

: .21093737E-Or 

NCFL.MCFL 

• 

1 

35 

SECANT 

METHOD 

IN  DECODE 

,  DID  NOT  CONVERGE  FOR 

N. 

N 

35  4 

SECANT 

METHOD 

IN  DECODE 

DIB  NOT  CONVERGE  FOR 

N. 

N 

35  4 

SECANT 

METHOD 

ZN  DECODE 

DID  NOT  CONVER 

GE  FOR 

H, 

N 

35  4 

SECANT 

METHOD 

IN  DECODE 

DID  NOT  CTNVERCE  FOR 

M. 

N 

36  4 

SECANT 

METHOD 

IN  DECODE 

DID  HOT  CONVERGE  FOR 

M. 

K 

36  4 

OSTEP 

- 

490 

EETA  - 

5.13669:57E40C 

DZETA 

- 

1.26054632E-C2 

NCFL.MCFL 

• 

1 

35 

OSTEP 

- 

500 

EETA  - 

5.265472i9E«0C 

DZETA 

- 

I . 31032II6E-02 

NCFL.MCFL 

- 

1 

35 

1  STEP  -  500  2ETA  -  5 . 26  54 '72  9E*0C  D2ETA  •  1.3:032m-02 

MACH  NO.  -  2.00COCOOE-*C1  ALPHA-  2 .  OOOOOOOE-01  BETA-  0 .  C  OOS  LCCE  *  C  L 

0 

0  ZONE  NO.  1 

0 

LEFT  EDGE  PLANE  GEOMETRY 


N 

X 

Y 

5G 

SGZ 

SG5 

ve: 

VE2 

VE3 

:s; 

36.5 

3.944509e-»00 

0.000000E«00 

0.0000O0E-»0C 

0 . OOOOOOE«OC 

C.OOOOOCE-CO 

0.000000E*00 

I .27I564E-03 

O.OOOOOOE-00 

0 

31.5 

3.462403E«00 

0.000000£*00 

O.OOOOOOE«OC 

C . OOOOOOE-CC 

C . OOOOOOE-OC 

0. OOOOOOE-CC 

1 .271566E-C3 

O.OOOOOOE-OC 

c 

25.5 

2.479957E*0C 

O.000000E*0C 

0 . OOOOOOE-CC 

C.OCOCOOE-CC 

C . OCCOOOE-OC 

O.OOOCCOE-CC 

1  271563E-C3 

0.0C0000E«0C 

c 

19.5 

2.297110E*00 

C.0C0000E*0C 

O.OOOOOOE*CC 

O.OOOCCCE-CC 

C.OOOOCCE^OC 

C-OOCOCOE-OC 

: .271566E-C3 

O.OCOOOCE-OC 

c 

13.5 

1 . 714264E*0C 

O.OOOOOOE*OC 

O.OOOOOOE«OC 

O.OOOOOCE-CC 

C.OOCOCOE-^CC 

O.CCOOOOE-CO 

I  .27;566E-03 

0. OOOOOOE-CC 

c 

7.5 

1 . 131417E40C 

0.O00000E*0C 

0. OOOOOOE-CC 

O.OCOOCOE-Of 

C . OCOOOOE-00 

C.OOOOOOE-OC 

I .271566E-03 

0 . OOOOOOE-CC 

c 

1  .  5 

S.445706E-C1 

C.COOOOOE-00 

O.OOOOOOE-OC 

C  OOOCCCE-CC 

0 . OCCCOCE-CC 

c.oocoocE-c: 

1  .271565E-03 

C , CCOOOCE-CC 

c 

LEFT 

’  EDGE  PLANE 

CEF  -  3.9970794E»00  BEF  - 

5  .  COCCCOOE-C: 

N 

S 

T 

P 

c 

u 

V 

N 

E 

36.5 

3.944S09E«00 

O.OOOOOOE«00 

5.273677E-C4 

3  .  C19022E-  Cl 

2 . 569:73E-C3 

C.OCCCCOE-OC 

5.22450CE-03 

' , 726:62E-C3 

1 .2r 

B-14 
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31. S 

3.442t03E400 

0.000000E-<-00 

2.202334E4^04 

6.423993E-02 

2 

25.3 

2.479957E+00 

O.OOOOOOE-fOO 

2.50UI2E«04 

3.115077E-02 

2 

19.5 

2.297110E<»00 

O.OOOOOOE^OO 

2.7914l3E-»04 

2.495513E-02 

1 

13.5 

1.714244E-»-00 

0.000000E400 

3.215629E^04 

2.529749E-02 

1 

7.5 

1 . 131417E400 

O.OOOOOOE+00 

2 .793257E+04 

2.115690E-Q2 

4 

1.5 

5.4I5704E-01 

O.OOOOQOE-»00 

1 .920391E«04 

1 . 51t229E-02 

-7 

CELL 

C  - 

PLANE  M  -  1 

3.9970794E4^00  CZ  • 

5.46307SEE-01 

CT  - 

-2.143t050E-01 

B  • 

S.OOOOOOOE'Ol  BZ  - 

O.OOOOOOOE^'OO 

BT  - 

o.ocoooooe^oo 

N 

S 

T 

P 

D 

37.0 

3.997079E-»00 

4.343323E<02 

5.529t74E«03 

I.I914I5E<Q2 

2 

34.5 
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